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FOREWORD 


The value of microwave scai•^erometers and radiomehers as remoi'e sea wind 
sensors has been independenhiy demonsfrai’ed by a number of inves^igai•ors. However, 
near-simulfaneous observaHons by a composi^e radipme^er and scafferomefer 
(RADSCAT) insf•rumen^ have been proposed as a me^hod of making improved wind 
es^^mates, The improvement' is derived from fhe complemenfary as well as 
supplementary features of the sensors. To demonstrate this potential a joint effort 
among New York University, General Electric Space Division, the University of 
Kansas and NASA Langley Research Center was undertaken through the Advanced 
Application Flight Experiment program of NASA. This document is submitted in 
support of these efforts . 

Specifically, this document was prepared by the Remote Sensing Laboratory 
of the University of Kansas Center For Research, Inc. under NASA Contract NAS 
I-'I0048. The principal investigator under the contract is Dr. R. K. Moore and 
project engineer Is Dr. A. K. Fung, Several individuals rendered valuable 
assistance in the development of the computer programs. The control and integration 
of the routines were partially achieved by Glen E. Elliott of the University of 
Kansas. John L. Mitchell of LTV operated the programs repeatedly on the CDC“ 
6600 while they were under final scrutiny by the author. 



ABSTRACT 


The theory^ design and operaf'ion of fhe compui'er programs which 
aufomafe fhe reducHon of join}- radiome^er and scaH'eromefer observations conducted 
by the AAFE RADSCAT instrument is presented. The programs reduce scatterometer 
measurements to the normalized scattering coefficient; whereas the radiometer 
measurements are converted into antenna temperatures. The programs are both 
investigator and user oriented. Supplementary parameters are provided to aid In 
the interpretation of the observations. A hierarchy of diagnostics is available to 
evaluate the operation of the instrument, the conduct of the experiments and the 
quality of the records. 

General description of the programs and their data products are also presented. 
This document therefore serves as a user's guide to the programs and is therefore 
intended to serve both the experimenter and the program operator. 
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I. INTRODUCTION 


The AAFE RADSCAT* Data Reduction Package, as prepared by the University 
of Kansas, is a family of computer programs** which together form a basis for 
automated reduction and presentation of scatterometer and radiometer measurements 
taken by the RADSCAT instrument. The primary member in the family validates the 
raw RADSCAT data records and uses instrument and aircraft parameters to yield 
emission and scattering characterisrics of the scene represented by the data. The 
scatterometer data is reduced to the normalized scattering coefficient whereas the 
radiometer data is reduced to the antenna temperature. Other programs in the 
family are employed to maintain a file of instrument characteristics, list the output 
(reduced data), and prepare duplicate output tapes for use on other computing machines. 
Two engineering routines are also provided to compute certain antenna parameters 
essential to the reduction of the microwave data; these routines are not employed 
regularly. The programs were written in FORTRAN IV, developed on a HW 635 
system and then modified and demonstrated on a CDC-6600. 

The program outputs provide a comprehensive basis from which to interpret 

r 

the results of RADSCAT experiments. The basic output parameters have been 
augmented with others to reconstruct the influence of the aircraft attitude on the 
experiment, to validate the measurements, and to reflect the status of the RADSCAT 
instrument. Full evaluation of the data is accomplished by examining contents of 

t 

an activity report as well as the output records. A fuller data base from which to 
analyze RADSCAT data is only provided when these outputs are augmented with seq/ 
ground truth and certain aircraft parameters available from other sources. The 
programs are therefore largely experimenter oriented. 

Attempts have been made to generalize the reduction algorithms so as to 
easily adapt the programs for changes in the character or operation of the instrument. 

In doing so the user is expected to provide many more of the calibration and 


A composite rad iometer and scatterometer sensor. 

The family of programs does not include a merge and edit routine which prepares 
the raw input data. This routine was developed by the Flight Instrumentation 
Division of NASA Langley Research Center. 



da^a paramefers than vvouid otherwise by required. For this reason a historical file 
is maintained on the instrument parameters. It Is also provided to trace the develop- 
ment, modifications and aging of the instrument. 

No attempts have been made to develop efficient coding; rather extensive 
efforts have been made to develop programs which will, to some degree, document 
themselves. Seif documentation is provided by a liberal insertion of comment cards. 
This approach had been instituted with the full realization that the programs are 
likely to undergo changes as experience or changes in the instrument dictate them. 

It is hopeful that this approach together with this document will ease the insertion 
of program modifications. However, the authors can assume no responsibility for 
changes made without their knowledge. It is strongly recommended that changes to 
the reduction algorithms be made only with the approval of an engineer who is 
thoroughly familiar with the operation of the instrument, the theory behind the 
algorithms and the needs of the investigators. 

Some strides have been made to design production type programs. Externa! 
directives give the users the ability to select the files to be processed or the service 
function to be executed. Unrecognized' directives will halt the execution of the 
program. Data entries are qiso validated and when too' many bad files are encountered 
the job is aborted. Activity reports and termination statements are provided so that 
the user can trace the progress of the job . 

The details on the use, theory and design of these programs are presented in 
the remainder of this document. To place these programs in context with RADSCAT 
experiments the technical background is briefly described In Section II. The , 
experimenter and programmer will find Section III helpful in understanding the 
general features and data products of the programs. He may wish to refer to Section 
IV and the Appendices for more specific information. The program operator will 
find Section IV essential to the operation of the programs. Again he may wish to 
consult the appendices for specific help. The cognizant engineer will want to 
become familiar with all aspects of this document to develop a thorough understand- 
ing of the program content. To a large extent this document assumes some knowledge 
of the operation of the RADSCAT instrument*. The programmer and cognizant 

' 

A description of the instrument may be found in "Field Service Procedure Handbook," 
vol . 1, prepared by General Electric Space Division, Philadelphia, Pennsylvania 
under NASA Contract NAS 1-10161, 
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engineer wHI also wani" to consult the document which describes the merge 
and edit routine which prepares the raw data records. This document Is In 
preparation at NASA Langley Research Center. Those wishing to learn the 
details on the theory and design of the various routines in the family of programs 
will find the appendices helpful. 
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IL TECHNICAL BACKGROUND 


A. InfroducHon 

This section is primarily written for the investigator or programmer not 
familiar with scatterometric or radiometric measurements. The motivation for joint 
measurements is presented and the elementary theory behind the measurements is 
discussed. The instrument and Its relation to the C-I30 aircraft is also briefly 
described. It is hoped that this section will place the AAFE RADSCAT Data Reduction 
Programs in context of experiments. 

B. Motivation for Joint Measurements 

In recent years a strong interest has developed in interpreting the microwave 
signatures of natural surfaces. In these efforts a number of potential earth resources 
and scientific applications have been demonstrated by measuring the emission or 
scattering properties of various scenes. The scatterometer is employed to measure the 

scattering characteristic whereas the radiometer is used to measure the emission 

« 

properties. To date the investigations have been largely conducted with one instrument 
or the other. However, several Investigators have proposed that joint scattering and 
emission measurements would form a strongerTnterpretational basis, because each 
sensor measures a different aspect of the bi“static scattering characteristic of the 
surface. From an entirely different Viewpoint the sensors are complementary In nature. 
It is known that the radiometer is very sensitive to clouds and rain whereas the 
scatterometer is relatively insensitive to clouds and rain. Over the ocean this 
property has potential in correcting the scotterometer measurements for small 
attenuations caused by clouds. 


C . Description of the RADSCAT Sensor 

The AAFE RADSCAT instrument is a composite sensor consisting of a pencil" 
beam microwave scatterometer and radiometer. It was designed to conduct measure- 
ments from an aircraft and is capable of operating at either of two frequencies, 

9.3 GHz or|I3.9 GHz. In three of Its four’modes of operation it makes scatterometer 
and radiometer measurements on a near-simultaneous basis. In the fourth mode it 
makes radiorneter rneasurements only. Observations can be conducted at any one 
of six adjustable angles covering the sector between nadir and 60 degrees. 

4 



The AAFE RADSCAT ins^rumen^ was primarily designed fo perform experimenis 
over the ocean to verify the remote anemotnetric capability of the composite sensor. 
Although primarily designed for oceanographic work, the extension to observations 
of agrarian or urban scenes is straightforward. The instrument is capable of conduct” 
Ing .experiments at one of four altitudes, viz, 2000, 5000, 10,000 and 20,000 feet.* 
The wide range of altitudes is provided to study atmospheric effects such as humidity, 
temperature, clouds and rains. To avoid difficulties introduced by operating from 
behind a radome, the RADSCAT instrument mounts on the tail gate of the NAS//MSC 
C-130 aircraft. The Litton Navigator LTN”51 aboard the C~130 aircraft is extremely 
useful to the interpretation and reduction of the RADSCAT observations. It provides 
aircraft attitude and location information from which the point of observation on 
the surface may be established. The attitude parameters are especially useful when 
measurements are performed from a roll maneuver, since these parameters enable 
one to compute the incident angle. The radar altimeter information provided by 
the aircraft is essential to the inversion of the scatterometer measurements. 


D . Estimating the Scattering and Emission Properties from Measurements 

The scatterometer is a radar which is capable of accurately measuring the 
backscatter properties of a rough surface. The backscotter ability of the surface is 
denoted by the normalized scattering coefficient which is given in terms of the 
scattering cross-section per unit area. In general it is dependent upon the point 

t 

of observation, the incident angle, frequency, polarization, etc. It is well known 
that the radar return is predicted by (see Figure 1), 


P. = 


(4 'T)' 


P. 


iJ- 

-^Tr 


3 ( 0 ,^) 0 - 


d A 


0 ) 


where = normalized scattering coefficient 

= received power 
= transmitted power 
X = wavelength 
g = antenna gain function 

I 

R = radar range tolelemental area dA 


The radiometer by itself can operate at any altitude. 
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For narrow (or equivalenfiy narrow) beam s c at fero meters the above integral 
expression may be inverted rather simply for <r°. As a consequence the problem of 
accurately determining a° reduces to the problem of accurately measuring the 
transmitted and receiver powers, the pattern function, look angle, and radar range. 
The area of illumination may be shown to be related to the pattern function, look 
angle, and radar range. The scatterometer provides outputs which are measures of 
P and P.. The power transfer function of scatterometer completes the association 
between the outputs to yield P^ and P^^ . A member of the family of programs 
actually employs the scatterometer transfer function and an algebraic approximation 
of equation (1.) to estimate o'®. 

The radiometer, on the other hand, is a very sensitive receiver capable of 
measuring very low level emissions which emanate naturally from objects. In the 
microwave region the natural emission is governed by the Jeans-Rayleigh law which 
states that the power emitted is proportional to the physical temperature T of the body. 
In actuality, however, different objects at different temperatures emit different 
powers depending until their relative ability to absorb energy. Better absorbers 
emit more efficiently than poor ones. The relative ability to emit is described in 
terms of a parameter called emissivity e . The emitted power per unit bandwidth 

* A 

in direction 0 is given by 


P(0) 


4tt k € (Q) T 


( 2 ) 


where k is Boltzmann’s constant. 

We may now suppose that the surface emits in accord with this relation. Then a 
radiometer whose antenna gain function is given by g(9,<l’ ) will measure (ignoring 
atmospheric contributions) a total power per unit bandwidth given by 


’’a = 7T JJ 9(e,4>) k Tg<e) dil 


Actually only the ratio P/P. is needed. 

** r r 

For one polarization. 

This relation is only representative of the nature of the problem. See Claassen, J, P. 
and A. K. Fung, "An Efficient Method for Inverting Antenna Temperatures for the 
Apparent Temperature Distributiloh" University of Kansas Center For Research, Inc., 
TR 186“8, January 1973, for a precise relationship. 
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where 


Jf 9 dil (4) 

4*rr 

• and 

Tg = € (9) T 

1 ' 

The anfenna power may be relafed ^o antenna temperature by For 

narrow beam efficient antenrias the antenna temperature is an approximate measure 

of the brightness temperature Tg In the look direction. The radiometer therefore 

must accurately measure the antenna power. The RADSCAT data reduction programs 

use the calibration data and radiometer transfer function to estimate T for each 

a 

polarization . 
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Ill, GENERAL DESCRIPTION 


A. InIroducHon 

This secHon is primarily wrif■^en lo familiarize fhe invesfigafor or programmer 
wifh the family of RADSCAT data reduction programs. The functions of the particular 
programs and the relationships between the programs are presented. Input and output 
data products are described to various degrees. 

B. Summary of Programs 

i * 

The AAFE RADSCAT Data Reduction Package consists of a family of six 
programs: 

1. Conversion Program (CONVERT) 

2. Output Program (OUTPUT) 

3. ProgVam for the Storage and Retrieval of Instrument Characteristics 
(ICHAR) 

4. Output Translation and Duplication Program (TDUPE) 

5. Equivalent Beamwidth Program (WIDTH) 

6. Antenna Gain and Efficiency Program (GAIN) 

The actual reduction of data is accomplished in CONVERT. The input data for 
CONVERT is prepared by an edit and merge routine developed by NASA Langley 
Research Center, The merged tape consists of RADSCAT and aircraft data. To 
assist in the computations, program ICHAR prepares a historical file of instrument 
characteristics. CONVERT selects a particular set of characteristics from the file 
and applies them to the input data. The computed results are stored on an output tape. 
Program OUTPUT will list the contents of the output tape. To prepare duplicate 
output tapes for use on other computing machines program TDUPE may be employed. 

WIDTH and GAIN are special routines which compute several antenna 
parameters. The parameters are employed as elements of the instrument characteristics. 
Both programs require antenna pattern information to compute the parameters. Their 
use will seldom be required unless the antenna pattern changes. 

Of the six programs the Initial four programs are production types and the 
lost two are engineering application programs. 
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C, Conversion Program 

The reducHon oF RADSCAT dai'o is performed enHrely in CONVERT. By 
means of exfern al direcfives, physical files of a specified fype ore processed with 
a designofed sef of insfrumenf characferisfics wifhdrown from a separafe file source. 
Exfernal direcfives permif fhe user fo inifiafe processing of inpuf files by fype or 
sequenfially regardless of fype. The files are organized by fiighf runs and are 
labeled in accord wifh a procedure insfifufed by the Flight Operations Division of 
NASA/MSC. 

Directive options are given by MISFLT, FLTLIN, SPECIF, or ALL, When 
the first directive is appended with mission and flight numbers, all files from that 
mission and fiighf are' processed in the order that they appear on the tape. Mission, 
flight and line numbers must appear wifh the second directive option. Files of 
the designated line type will thus be processed. Consecutive directives which 
exhaust the line types will order the output files by line type. This feature is 
helpful to the interpreter who often associates an experiment condition with a 
line type. The third directive is employed when a specific file designated by 
mission, flight, line and run is to be processed. This directive is appropriate for 
a quick or selective look at the experiment results. The last directive (ALL) will 
initiate processing of all files regardless of type. Combinations of directives may be 
employed; however, indiscriminant use of the processing types is cautioned against 
since duplicate processing of files may occur. For each directive the file search is 
effectively conducted from the beginning of the tape. 

Each of the above directives must include a designator which selects a set 
of instrument characteristics from a historical file of characteristics or which permits 
them to be read from cards. It Is likely that the most current characteristics will 
be chosen; however, any set on the file may be selected. Alternatively, characteristics 
may be entered from the card reader. However, operation under this latter option is 
restricted to one data file per directive. Entry through the card reader is mainly 
reserved for testing a set of instrument characteristics before storage on the historical 
file. Also as a part of the same directive the number of files to be processed of the 
indicated type and the line density of the activity report must be designated on the 
control card. The latter parameter can be used to limit the number of data records 
on which reports are given. 
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During processing CONVERT genera^es a conversion report. The report 
routinely issues statements on the colibration parameters and critical receiver 
temperatures as they are computed from the input file. It will also routinely state 
that 0 record had been reduced. The reports are appended with record numbers. 

Other special statements are generated whenever critical situations are encountered 
in the records. When these special situations occur, the record is normally processed 
regardless; however, that record is flagged by an interpretation a I code which becomes 
a part of the output record. In some cases, such as invalid data, the record is by- 
passed and an appropriate error message is generated. When an excessive number of 
bad records is encountered', the file is bypassed and an appropriate message printed. 
These messages are primarily intended to reflect the status of the processing, the 
condition of the instrument, the conduct of the experiment, and the quality of the data. 
A sample listing of the conversion report is illustrated in Figure 2. 

As processing progresses output records are accumulated in blocks (arrays) and 
the blocks are written to the output tape, A description of the contents of the output 
tape appears in the following subsection, 

D . Output Program 

Program OUTPUT will list the processed records produced in the conversion 
program. External directives similar to those In CONVERT permit listing files of 

designated types or ail files regardless of type. Again general directives are given 

* 

by MISFIT, FLTLIN, SPECIF or ALL. The appropriate mission, flight, line and run 
designators must accompany the directives. 

A sample listing produced by OUTPUT is shown in Figure 3, A double page 
of output consists of the file labe!,column headings, 80 records and a key. The sel- 
of instrument characteristics which produced the results is Identified by set number 
and date in the last entry of the file label.' Each measurement record occupies one 
line on the listing. A record number is appended to each line and should agree with 
the number assigned in the conversion report. The time of measurement in mean 
Greenwich time is presented In the second column. It is to be redd in hours, minutes, 
seconds and 1/IOth seconds. The measurement can therefore be associated with data 
from olher sources. The mode from which the instrument conducted the measurement 
is reflected in column 3. The incident and cross-tracks angles appear in the next two 
columns. They were computed from aircraft attitude parameters and the RADSCAT 

view angle with respect to the aircraft platform. The cross-track angle is that angle 
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V 

V 

• 595 

5.39826*00 

7.3 

243.3 

01000 

lOOOOO 

9958 

l.l 

bO 

151928.0 

F.A. 

13.1 

-1.4 

V 

V 

.595 

4.66996*00 

6.7 

244.2 

O.OOQ 

LOOOOO 

9961 

1.2 


KEY TO COOES- 

0 KC MAO 1000 POLAR REVERSAL 

1 FCSSItiLY OOTSJOE KA^OE CATE lOOOO EXCESSIVE oEPUL 

10 CUSIOE DYJsANiC range lOOOOO KEC TEKP AuACKHAL 

103 CXCESSJVC OUPPLEA tOOOOOO kARK LOAD TEPP ESlIKATEU 

^ CIPUIKA IJlAS Op plains CAN CCCUK. 



FIGURE 3. 


SAMPLE LISTING OF OUTPUT RECORDS 
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be^ween I’he "negat'ive" flighf- direci-ion and the azimuthai look direction. The 

positive cross“track angle is defined to the right of the aircraft when looking aft. 

It is induced by aircraft roll and drift angle. The incident and cross-track angles 

are precisely computed for any aircraft attitude. The angles enable the investigator 

to locate the illumination area. The scatterometer transmit polarization and the 

scatterometer and radiometer receive polarization are presented in the next two 

columns, respectively. The scatterometer integration time is indicated in the 

following column. The radiometer integration time is not listed since it is constant 

(128 milliseconds) regardless of mode or angle. The normalized scattering coefficient 

2 2 

(m /m ) is shown in decimal as well as dB units. The radiometer antenna temperature 
Is' shown In column 11. 

A depolarization factor appears in column 12. The depolarization factor when 
expressed as a percentage Indicates the percent power transmitted at the opposite surface 
polarization. Whenever DEPOT exceeds 50 per cent a polarization reversal flag is 
generated under the S/R validation code (column 13) and the complementary depolari- 
zation factor (1-DEPOL) is posted in column 12. If a reversal Is indicated and the 
(complementary) depolarization factor is small, the polarization designators in 
columns 6 and 7 may indeed be regarded in the opposite state. Polarization 
reversals will frequently occur when observations ore conducted from an aircraft 
roll condition. When DEPOT exceeds 2 per cent a flag in the S/R code is raised. 

Other flag conditions can occur under the S/R code. The flags are used as 
a quick indicator of an abnormal condition. The interpretations of flags are given in 
a key at the bottom of every double page. The range gate flag is specifically set 
whenever the aircraft is not within + 300 feet of the range gate setting. On some 
occasions this condition may be purposely generated to check the scatterometer 
zeros. A dynamic range flag is set whenever the scatterometer return Is above or 
below the dynamic capability of the scatterometer receiver. If any of the receiver 
tempera tures;are not within the proper operating range a flag is set.' Whenever the 
receiver temperature data is missing for long periods of time (data drop out), the 
receiver temperatures are estimated and a flag is generated. The other keys are 
self-explanatory. For additional assistance In diagnosing the flagged condition 
the conversion report may be consulted. 

The latter two columns denote aircraft radar altitude and flight direction. 

Other aircraft parameters are available from the listings of the aircraft mission tape. 
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Nof lisl'ed bui' presenf on fhe fape is l^he ambien^ aif temperat'ure a^ alfifude. .This 
la^^er parameter is available should fhe inveshlgai-or wish fo use i^ ho eshimahe hhe 
apparenh temperahure disfribuhion* . 

E. Program ICHAR 

Program ICHAR maintains a historical tape file of RADSCAT instrument 
characteristics. The historical file will (eventually) consist of many sets of instrument 
parameters which have been entered chronologically and labeled numerically. The 
historical file is maintained to chart the changes in the instrument parameters os may 
occur through updating, instrument modifications or simply aging. The Conversion 
Program withdraws an appropriate set of parameters from the historical file and applies 
the withdrawn parameters to reduce the measurements, 

ICHAR performs several useful functions. It will initially establish the file. 

Once the file has been established it will permit the user to append new sets of 
characteristics. To assist in the maintenance of the file, ICHAR may be used to 
print and punch old sets of characteristics. The punched sets may be modified, for 
example, and submitted to ICHAR as a new set to be appended to file, 

F, Program TDUPE 

Program TDUPE provides a master duplicate tape of the output tape formed by 
the Conversion Program for use on other computing machines. The bit conversion 
problem between machines is avoided by scaling and integerizing all the numerical 
data before transferring the output data to the duplicate tape. This technique has 
been successfully applied to tapes prepared by a CDC-6600 and read on a Honeywell ■ 
635. The CDC 6600 is a 60 bit machine which performs, "ones complement" arithmetic 
whereas the Honeywell is a 36 bit machine which performs "twos complement" arithmetic 
The scale factors may be withdrawn from the listing of TDUPE shown in Figure D~2 . 

Each output file consists of a header label (a flight run), the Instrument characteristics 
applied to the data, and the output data records. 

_ ^ ^ ^ . 

See for example, Claassen, J. P. and A. K. Fung, "An Efficient Technique for 

Estimating the Apparent Temperature Distribution," University of Kansas Center for 

Research, Inc., Technical Report 186-8 , January 1973. 
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Additional inlegerized tapes ma/ be prepared from fhe masi'er dupltcal-e ^ape 
by using syslem uhilify programs. 

I 

G . Special Programs 

In addition to the above production type programs two engineering routines 
WIDTH and GAIN are included in the reduction package. Program WIDTH computes 
the equivalent beamwidth when it is presented with the mainbeam antenna pattern. 

The equivalent beamwidth parameter is essential to the inversion of the scatterometer 
data. The program GAIN computes the beam solid angle and directivity. The initial 
factor is important to the inversion of the radiometer measurements whereas the second 
is required in the inversion of the scatterometer measurements. 

Functional representations of the antenna pattern are required by both programs. 
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IV, OPERATING INSTRUCTIONS 


A, IniroducHon 

General operating requirements, preparation of data input files, construction 
of the job deck, descriptions of the data products, flags and aborts, etc., are treated 
in this section. It is primarily written for the user having little or no knowledge of the 
theory and design of the programs. The reader will find that Section II and Section 
III will be helpful in understanding the contents of this section. Supplementary as 
well as detailed information on these programs may be found in the Appendices. 
Sufficient information is intended to permit the user to operate the primary programs 
CONVERT, OUTPUT, ICHAR and TDUPE. The remaining programs WIDTH and GAIN 
are engineering routines and should be handled by an engineer (or equivalent) familiar 
with the RADSCAT antenna. However, It is also strongly advised that a competent 
engineer, thoroughly familiar with the RADSCAT instrument,be involved in the 
preparation of the data and in the operation of the primary programs. 

System control cards are not specified. It is presumed that this section together 
with sufficient knowledge of the operating system will dictate the necessary system 
control cards. 

B. Convers ion Program (CONVERT) 

1 . General 

This program reduces 50 word input records to 20 word output records. Various 
computational algorithms described in detail in fhe Appendices are applied to the 
input records to yield data products of interest to the investigator. The input records 
are members of a file having an indefinite number of records . The files are elements 
of a multi“file tape. Each file is labeled with the entries shown in Table 1 . The 
files are footed with an "end of file" mark. It is anticipated that files (flight runs) 
from several missions can be stored on a single tape. Tapes of this type are designated 
Raw Data Tapes. 

To Initiate processing of any of the files or all of them, processing directives 
are employed. The directives permit seeking files of certain types and must appear 
as data cards In the job deck. Each processing directive must be accompanied with 
a file referral argument to select a set of instrument parameters from the Instrument 
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Characterishics Tape* on which se^s of insfrumenf paramefers ore indexed numerically 
{and presumably sfored chronologically), A sef of these parameters is employed by 
CONVERT to reduce the row data records. Alternatively, the instrument parameters 
can be entered from the card reader should the file index be set at zero. Files 
processed under a particular directive are transferred to the Output Tape. Various 
processing directiyes may be applied successively. In each instance the file search 
is effectively initiated from the beginning of the tape. Processing continues until 
ail directives have been exhausted. 

Other types of directives perform service functions. The directive POST 
permits newly processed files to be appended to an old Output Tape. This directive, 
when used, must precede all processing directives. The PRINT and PUNCH commands 
will list and punch, respectively, a set or several consecutive sets of instruments 
parameters from the Instrument Characteristics Tape, When the PUNCH command is 
employed, a listing of the set is also generated. 

The composition of a CONVERT job deck is illustrated in Figure 4, 

2. Preparation of the Row Data Tape 

As mentioned above each file is identified by a header label . The first 
seven items in the lobe! (Table 1) are integer variables. The integer variables must 
be right justified. Note that the date requires three words. The sixth and seventh 
items are appended to the label to assist in selecting the proper subset of instrument 
characteristics from the set index designated on the processing directive card. See 
Section IV D for a description of the instrument characteristic files. The frequency is 
specified in GHz, The feed type is declared in the seventh word and is an alpha" 
numeric word. 

Presently two feed types are employed by the RADSCAT Instrument. One feed 
uses a mechanical waveguide switch and the other a ferrite circulator. The 
abbreviations WG SW and CIRC are recommended as designators. To establish a 
convention, they should be left justified. The label must be created with an unformatted 
WRITE statement. 

Beneath the header label are a sequence of raw data records each a mixture of 
real and Integer words. The records are blocked by groups of 10 records. The blocks 

— ^ 

This tape is physically separate from the Raw Data Tape. 
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SYSTEM CARDS 

lNSTRU(ViENT CHARACTERISTICS 



IX POST 



figure 4. 


^OR PROGRAM CONVERT 


19 




Table 1 

. File Label Entries Illustrated 


Position 

Name 

Entry 

1 

Mission 

207 

2 

Flight 

n 

3 

Month 

7 

4 

Day 

21 

5 

Year 

1972 

6 

Line 

4 

7 

Run 

9 

8 

Frequency 

13.9 

9 

Feed 

WG SW 


mus^ be created by unformatted WRITE statements. When insufficient records exist 
to fill the last block on a file,^ the remaining records must be filled with 99999.0 
(floating point). These entries are used to flag the end of the records (file). In 
addition on EOF is appended to the foot of the file. As many files as practical 
may be created on the tape. The end of the files must be flagged with the label 
END OF REEL which fills the beginning of a nine word array and in effect is another 
file label . 

The files are corhposed of three record types: 

1. Normal Calibration 

2. Baseline Cal ibrafion 

3. Measurements 

Calibration records may appear anywhere in the file and as frequently as desired. In 
a properly constructed file both types of calibrations should appear in the first 
records. However, should one (or both) type(s) of calibration be missing from the 
top of the file, the program will search for calibration records deeper in the file. 
Records bypassed in the search are transferred to a scratch unit for subsequent 
processing. Once found, the calibration information is extracted and applied to the 
bypassed records. Records following the calibration records are also treated with 
the same parameters until new calibration records are encountered. When the cali- 
bration information is totally missing or incomplete, the file is bypassed and ah 
appropriate error message is printed. 

The entries in the records are shown in Table'll. Type and units of the 
variables are indicated. Parenthetical designators behind some of the entry names 
describe the computed results stored In those positions during processing. The Integer 
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TABLE II RAW DATA VECTOR 


NAME ; DATA 
POSITION 

ENTRY 

TYPE 

UNITS 

I 

TIME 

R 

HR MIN SEC 

2 

FREQUENCY INDICATOR 

I 


3 

MODE INDICATOR 

I 

— 

4 

DATA TYPE INDICATOR 

I 

— 

5 

ANGLE DESIGNATOR 

I 

— 

6 

TRANS POL INDICATOR 

I 


7 

REC POL INDICATOR 

I 

— 

8 

SCAT 1 (OR CONVERTED SCAT) 

R 

VOLTS 

9 

SCAT 2 

11 

ti 

10 

SCAT 3 

ir 

II . 

11 

SCAT 4 

It 

II 

12 

RAD POL INDICATOR 

I 

— 

13 

RAD (OR CONVERTED RAD) 

R 

VOLTS 

14 

SCAT'ANGLE (OR INCID ANGLE) 

R 

II 

15 

RAD ANGLE (OR X-TRACK ANGLE) 

R 

ti 

16 

RANGE GATE INDICATOR 

I 

— 

17 

T^, 9.3 DICKE SWITCH 
TRIPLEXER ' 

R 

VOLTS 

18 

R 

11 

19 

Ty T/R CIRCULATOR 

R 

M 

20 

T^, 13.9 DICKE SWITCH 

R 

II 

21 

T^, 9.3 TDA 

R 

11 

22 

T^, 9.3 WARM LOAD 

R 

11 

23 

Tyr 9.3 LIMITER 

R 

It 

24 ■ 

Tg, BRACKET 

T^, COMMON T/R GUIDE 

R 

11 

25 

R 

11 

26 

T^q, 11.7 DICKE SWITCH 

R 

n 

27 

T^^, 13.9 WARM LOAD 

R 

11 

28 

polar switch OR CIRC. 

R 

ri 

29 

T^3, scat calib switch 

R 

11 

30 

T^4, 9.3 IMAGE REJECTION FILTER 

R 

1) 
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TABLE II RAW DATA VECTOR (continued) 


NAME : DATA 




POSITION 

ENTRY 

type 

UNITS 

31 

T^^. baseplate 

R 

VOLTS 

32 

OMT 

R 

n 

33 

HOT LOAD 

R 

VOLTS 

34 

V 

R 

VOLTS 


age 



35 

^supply 

R 

It 

36 

FLAG (T^^ MISSING) 

I 

”” 

37 

EMPTY 



38 

M 



39 

It 



40 

11 



4.1 

RADAR ALTITUDE 

R 

' FEET 

42 

GRD SPEED 

R 

FT/SEC 

43 

PITCH 

R 

DEGREES 

44 

DRIFT 

R 

DEGREES 

45 

ROLL 

R 

DEGREES 

46 

AMBIENT AIR TEMP 

R 

°C 

47 

FLIGHT DIRECTION 

R 

DEGREES 

48 

empty' 




49 

50 
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indical’or variables -are -defined in Table III, These indicafors are employed fo cue 
certain program actions or to index variables and must be right justified. Ifems 2 
through 13 and 16 are extracted from the RAD5CAT microwave data channel. Items 
17 through 32 are retrieved from the RADSCAT ATM channel . Item 33 is presently 
entered by hand but may eventually appear on the ATM channel. Items 34 and 35 
are presently ignored but may also eventually appear on the ATM channel . Item 36 
Is presently entered by hand and is used to flag missing temperature data (DATA(17) - 
DATA(32)). Items 41 to 45 and 47 are provided by the Litton Navigator LTN-51 
aboard the C~130. Item 46 is provided by an external sensor aboard the aircraft. 

The measurement records must be completed as shown in Table II. However, 
the calibration records may be abbreviated to include only the data type indicator 
(DATA(4)) and the calibration entries in the appropriate SCAT ancj/or in the RAD 
position. The normal calibration records must be ordered by SCAT channels (one 
through four) with the normal RAD calibration occurring in the fourth record (contain- 
ing the SCAT channel 4 calibration). The routine which extracts calibration data 
will orient itself on the calibration record for SCAT channel 1 In an attempt to meet 
contingencies where calibration record drop-outs occur. Any number of normal 
calibration sequences (calibrations on channels one through four in each sequence) 
may occur successively. All will be averaged. Typically the instrument produces 
three sequences of normal calibrations. The baseline calibration record simply 
requires the appropriate RAD baseline calibration in the DATA(13) position (as 
well as the type indicator in DATA(4)) and' may occur anywhere in the file except 
within a sequence of normal calibrations. It is essential that the baseline and normal 
calibrations be grouped by three or more consecutive records to assure that a sta- 
tistically significant average will be established. The program wilt average con- 
secutive baseline records, even if they are interspersed with sets of normal calibrations. 
The norma! calibrations will also be averaged in this circumstance. 

A routine prepared by Langley Research Center actually merges data 
from several sources to meet the above specifications. 

3, Preparation ond Use of the Control Cards 

There are four types of processing directives and three service directives. 

The directives consist of a control word and some arguments. All control words are 
left justified; the arguments are right justified. Sample directives and their formats 
are defined in Table IV. 
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TABLE m. DEFINITION OF INTEGER INDICATOR VARIABLES 



0 

1 

2 

1 

3 

4 

! 1 

‘ 5 1 

1 

! 6 

Frequency Indicoror 

i 

9.3 GHz 

11. 4* GHz 

13.9 GHz 


j 

! 

Dctc Type • 

Me CSV re merit 

Normol 

Cclibrotion 

Base 

Cclibral-ion 

1 

1 

1 




1 

Angle Designator j 


X 

X 

X 

X 

X j 

X 


Horizonte! 

Vortical 

[ 



1 

1 


Range Gate 

2000 ft. 

5000 ft. 

10,000 ft. 

20,000 ft. 



1 

Mode Indicc'or'** 

R.O. 

t 

S.S. 

F.A. - 

A. A. 





’not operational 


R. O'. = RADIOMETER ONLY 

S. S. = SHORT SCAT 
F.A. = FIXED ANGLE 

A .A . = ALTERNATING ANGLE 


























TABLE IV. CONTROL CARD OPTIONS ILLUSTRATED 


CONTROL WORD 

MISSION 

FLIGHT 

LINE 

RUN 

ARGUMENTS 
LINE PRINT MODULUS 

ICHAR SET NO. 

FILES TO BE PROCESSED 

A6, 4X 

i 

1 

15 

15 

15 

15 

15 

15 

15 

SPECIF 

207 

12 

2 

3 

1 

2 

1 


208- 




\ 



FLTLIN 

203 

1 

4 

-- 

5 

3 

81 

MISFLT 

209 

4 

. 

— 

10 

3 

10 

ALL 

— 

— 

— 

— 

10 

2 

999 

POST 

— 

— 

— 

— 

“ 

— 

— 

PRINT 

— 

— 

— 

— 

— 

2 

2 

PUNCH 

— 

— 

— 

— 


1 

1 



The conhroJ word for ihe processing direcHves designafes J-he cJassiflcoHon of 
f'he file search To be conducted. The search is conducfed by comparing fhe argument's^ 
mission, flight, line and run on the directive with the file labels on the Raw Data 
Tape, The depth of comparison in the label will depend on the type of processing 
directive. ALL will process all files on ihe tape, MISFLT will process all lines and 
runs under the designated mission and flight, FLTLIN will process all runs from the 
specified mission, flight, and line. SPECIF will process only a particular file 
labeled by the specified mission, flight, line and run. When consecutive processing 
directives are employed, the file search is effectively conducted from the beginning 
of the' tape . 

Arguments other than the mission, flight, line and run must also be provided. 
The sixth word LNPMOD controls the density of records on which reports will be 
mode in the activity report. The use of this parameter will limit the amount of 
print out acquired during the reduction of records. It will not affect the density of 
output records . 

The set of instrument characteristics (ICU) to be applied to the specified files 
is designated In the seventh word. The entry specifies the file position on a 
master tape of characteristics. If a zero is entered, the program will expect 
instrument parameters to be read from cards appended behind the program control 
card. When characteristics are read from the card reader, processing is limited to 
one physical file regardless of the directive. This control option is primarily reserved 
to test a set of instrument characteristics. See Section IV D for a description of the 
format required for the instrument characteristics. If 999 is specified,the last set 
of characteristics on the master (historical) tape will be used. 

The last (eighth) entry NFP on the control cord specifies the number of files 
to be processed under the directive. To process all files under the directive simply 
make NFP sufficiently large. 

The special service command POST may be used to position a previously 
produced output tape to a point behind the last file on the tape. The use of this 
command will, therefore, permit new files to be appended to an old output tope. This 
command, when used, must appear before any processing directives. There are no 
arguments with this directive. 

The two service directives PRINT and PUNCH are provided to list and list 
and punch a designated set of Instrument characteristics, respectively. The desired 
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sef number is enfered In columns 36lhrough 40 of fhe direcHve card. If n“I sefs 
behind the designated set are to be listed or punched also, NFP should be set to n; 
otherwise it is set equal to one. NFP appears in column 41 through 45 of the directive 
card, 

4. Peripheral and Memory Requirement 

The input and output unit requirements are illustrated in Figure 5. Under a 
processing directive units 01, 02, 04, 06, 08 are required; under the service directives 
PRINT or PUNCH units 04, 06, and 43 are required. Unit 08 is a scratch unit, either 
disc or tape, perferably disc. It is used to temporarily store bypassed records when 
it is necessary to search for calibration data. 

Approximately 17000 decimal words of memory were required to load this 
program on a CDC-6600 to include buffer space and system routines. 

5. Flags and Aborts 

The activity report produced by the program conveys various messages to ^ 
indicate the status of the processing. It is normal for the program to echo check 
the directives, to list the progress on the processing, and to indicate a normal 
termination. The Conversion Report will also make regular statements reflecting the 
progress of the processing and the quality of the measurements. A sample listing of 
the Conversion Report Is shown in Figure 2 of Section III. In some instances however, 
improper use of the program or construction of the data will result in an error message 
or an abort. 

When on abort condition is encountered, an error message is generated, an 
abort code is given, an end of file mark is appended to the output tape on unit 02, 
an END OF REEL statement is appended, and control is yielded to the system. 

Table V defines the various abort conditions. When invalid data records are 
encountered, an error message will identify the Improper element and the entire record 
will be dumped in octal. The message will appear within the Conversion Report. 

When more than 10 bad records are encountered in a file, the remaining records In 
the file are bypassed and processing is directed to the next internal activity. When 
more than three bad files are encountered the processing is aborted. 
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FIGURE 5. 


PERIPHERAL REQUIREMENTS FOR THE PROGRAM CONVERT 





TabI 

e V. Definition of Abort Codes 

UD 

= Unrecognized Directive 

IC 

= Instrument Characteristic Missing 

MF 

= Missing Raw Data File 

FF 

= Unable to Find Instrument Characteristics 
with the Feed and Frequency Designated 

ER 

= Too Many File Errors 

NFP 

= NFP less than 1 

6. Description 

of the Data Products 


Normally two data products are qenerated by CONVERT^ an activity report 
and an output tape. When the PUNCH option is employed a set or sets of punched 
Instrument parameters can also occur. The character of the punched output is 
Identical to that described in Section III D, The activity report reflects the program 
progress and largely contains items In the Conversion Report. See Figure 2 of Section 
III for an Illustration of the Conversion Report. The Output Tape, however, will be 
described. 

The output file structure is similar to that of the input raw data files. The 
file consists of the file label and blocks of output records. However, occurring 
between the label and output record blocks Is the set of instrument parameters that 
had been used to process the input records. The entries in the file label are the same 
as that for the input files shown in Table I. The label and instrument parameters are 
created with separate unformatted WRITE statements. The entries in the instrument 
characteristics are described in Section III D, The output records occurring beneath 
the instrument parameters are blocked in 20 records groups. Each block Is created 
with a binary WRITE statement. The elements In the output records consist of mixed 
integer and floating point entries. The entries are defined in Table VI, When 
insufficient records exist to fill ihe the last block in the file, the remaining entries 

i 

are filled with integer zeros. An "end of file" mark is created after the last block 
In the file. 

When all, directives have been exhausted or the program aborted an "end of 
reel" statement is appended behind the last output file. 
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TABLE VI OUTPUT DATA VECTOR 



POSITION 

ENTRY 

TYPE 

UNITS 

1 . 

TIME 

R 

HR, MIN, SEC 

2. 

MODE 

I 

— 

3. 

INCIDENCE ANGLE 

R 

DEGREES 

4. 

CROSS TRACK ANGLE 

R 

DEGREES 

5. 

TRANSMIT POLARIZATION 

I 

— 

6. 

RECEIVE POLARIZATION 

I 

— 

7. 

SCAT INTEGRATION TIME 

'R 

SECONDS 

8. 

NORMALIZE SCATTERING COEF. 

R 

mVm^ 

9. 

NORMALIZE SCATTERING COEF. 

R 

dB 

10. 

ANTENNA TEMPERATURE 

R 

DEGREES KELVIN 

11. 

DEPOLARIZATION FACTOR 

R 

— 

12. 

DATA VALIDATION CODE 

I 

— 

13. 

ALTITUDE 


FEET 

14. 

FLIGHT DIRECTION 


DEGREES wri 
TRUE NORTH 

15. 

AMBIENT AIR TEMP 


DEGREES KELVIN 

16-20. 

EMPTY 
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C. Qul^puf Program 


1 . General 

The Output Program simply lists the output records prepared by program 
CONVERT, The output listings are organized by files which are selected by externa! 
directives similar to those described for the Conversion Program, A sample of a list- 
ing is illustrated in Figure 3 of Section III. The file heading together with instrument 
characteristic identifiers compose the page heading. Beneath the heading are listed 
80 record images; each image consists of the first fourteen entries in an output record 
as defined in Table VI. The fifteenth entry is not printed but is available on tape 
should the investigator choose to use it to further reduce the radiometer data. Each 
line of print-out is appended with the record number. The record number corresponds 
to those used in the Conversion Report. The 80 records are intended to fill two 
consecutive pages of print-out where the page size is 8-1/2" by 15". Appended to 
the foot of each double page is a key defining the 5/R validation code in column 13. 

The instrument characteristics although present on the output tape are not 
listed. The set which has produced the data is Identified in the heading on the 
right by the set number and date. No option is provided to print the characteristics. 

The job deck for the program OUTPUT Is organized identically to the conversion 
program (however, no Instrument characteristics must appear in the data deck). File 
types are selected from the output by proper use of the directives: 

1 . ALL 

2. MISFLT 

3. FLTLIN 

4. SPECIF 

The meaning and use of these directives have been presented in Sections III B.l and 
III B.2. The format of these control cards is identical to those described above. In 
the argument field of the control card the appropriate file identifiers must be present. 
No provision is made to limit the number of files to be processed, nor can the record 
density be controlled through a line print modulus. It is, of course, meaningless to 
specify the Instrument characteristics ICU. 

2 , Program Messages 

In addition to listing the output file, OUTPUT will keep the user posted on 

Its internal activity as it searches for the specified files. The directive Is initially 

echo checked, bypassed files are identified and the user is notified when the reel Is 

encountered. Unrecognized directives will terminate the processing. 
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3 . Peripheral and Memory Requrrement 

Peripheral requirements are illustrated in Figure 6. The program required 
approximately 7000 decimal words of memory on a CDC 6600 to include buffer space 
and system routines. 

D. Program ICHAR i 
' ^ ♦ Generol 

The primary purpose of program ICHAR is to maintain a file of Instrument 
characteristics. The file is maintained not only to serve as a source of parameters 
to reduce raw data records but also to serve as a history of the changes undergone by 
the RADSCAT instrument. The file can consist of many sets of instrument parameters 
which are appended to the file numerically. The user may also identify the set by a 
date which is entered with the parameters. 

Each set of instrument characteristics can consist of four subsets. The four 
subsets permit characterizing the instrument for a combination of two feeds and two 
frequencies. Each subset is composed of a set of radiometer parameters and a set of 
scatterometer parameters. The subsets are identified by feed and frequency. These 
identifiers are used by the Conversion Program to select the appropriate subset of 
characteristics . 

This file maintenance program is operated by means of external control cards. 
Four directive (processing) options are offered: 

1 . INITIALIZE 

2. APPEND 

3. PRINT 

4. PUNCH 

They must be left justified and occupy columns 1 through 10 of the control card. 
INITIALIZE will create the historical file by storing the first set of characteristics 
on the tape. The set must appear behind the control card. APPEND will permit 
other sets to be added to the existing historical file. The characteristics set to be 
appended must appear behind the APPEND command card. When several sets are 
added at one operation, each set must be led by an APPEND card and intervening 
PRINT or PUNCH commands are not permitted. When sets ore appended to the file, 
the old sets are tronsferred from the old tape to a new tape and the new sets are then 
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FIGURE 6. PERIPHERAL REQUIREMENTS FOR THE PROGRAM OUTPUT 


m 


added to the new tape. Before appending a new set of characteristics, it is 
advisable to have tested them in the Conversion Program by entering them from the 
card reader. 

The PRINT and PUNCH comrriands perform the same service as they do in the 
Conversion Program. However, one command must be used for each set to be printed 
or punched. The set to be printed or punched must be specified on the control card’ 
in columns 11 through 15. The set number must be right justified . A listing 
accompanies the punched output. The PRINT or PUNCH may proceed or follov/ a set 
of APPEND commands but must not intermingle with them. These commands are 
provided so that new sets may be formed easily from old sets of characteristics. 

A sample deck is illustrated in Figure 7. 

2 . Preparation of the Instrument Characteristics 

The instrument parameters cards must be prepared in groups of four subsets. 

When insufficient data exists to complete one of the subsets, a complete subset 
may be repeated an appropriate number of -times to fill the requirement. A sample 
subset is shown in Figure 8.' The lines actually represent card images (except for the 
word SET). The entries in the upper half are radiometer parameters whereas the entries 
in the lower half are scatterometer parameters. This subset would be used for 
measurements taken at 13.9 GHz with a waveguide switch (WG SW) as the polarization 
switch. The radiometer and scatterometer parameters are Identified with the set 
number of which it Is a member and with a date of entry (origin). The first line of 
radiometer or scatterometer parameters must be punched with the format 415, lOX, 
E10.3, AlO, E10.3 whereas the remaining lines must-be punched according to the 
format 7E10.3. Excluding the set number and date, the radiometer parameters array 
consist of 75 words and the scatterometer parameter array consists of 85 words. Not 
all words are necessarily filled. 

The entries In the radiometer and scatterometer arrays are defined in Tables 
VII and VIII, respectively. Many of the entries can be clarified by referring to a 
description of the reduction algorithms in the Appendices. These algorithms have 
been developed from the transfer and calibration characteristics of the RADSCAT 
Instrument. As a consequence the entries may change (and perhaps the program 
algorithms) when the instrument Is modified. Many of the entries must be based on 
special calibrations and measurements. Explanations are required for some of the 
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FIGURE 7. JOB DECK FOR PROGRAM ICHAR 
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FIGURE 8. SAMPLE OF INSTRUMENT CHARACTERISTICS 
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TABLE VII RADIOMETER PARAMETERS 


VARIABLE NAME ; RPAR 


POSITION ENTRY 


TYPE UNITS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


FREQUENCY 

FEED (WG SW OR CIRC) 

CROSS FEED CONSTANT 

II II 

HOT LOAD REFERENCE CONST 
WARM LOAD REFERENCE CONST 
HOT LOAD BIAS FACTOR 
WARM LOAD BIAS FACTOR 
SWITCH COMPENSATION FACTOR 

(1 II It 

OMT 

II U II 

CUTLER FEED " - " , 

II It 11 

INPUT GUIDE " 

RESIDUAL 

II 

MEASUREMENT PERIOD 

CALIBRATION PERIOD 

UPPER TEMP LIMIT 

LOWER TEMP LIMIT 

THERMISTOR CONVERSION FACTOR 


"HV 

’VH 


a 


IH 

^IV 

^ 2 H 

^ 2 V 

"3H 

^ 3 V 

''4H 

^4V 

PSH 


" 5 V 


T 

T 


a 


C 

UP 

LOW 

'1 


a 


ct„ 


a 


4 

^4 


R 

A 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 

R 


GHZ 

UNITLESS 




UNIT LESS 
UNITLESS 


SEC 

SEC 

°K 

°cAolt 

°C 

°cAolt 


°CAOLT 

°C/VOLT 
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TABLE VII RADIOMETER PARAMETERS (continued) 

* 

POSITION ENTRY TYPE UNITS 


31 

THERMISTOR 

CONVERSION FACTOR 

“5 

F6.2 

°CAOLT 

32 

11 

II 

It 

^5 

11 

°C 

33 

II 

M 

M 

a 

6 

II 

°C/VOLT 

34 

II 

t 

II 

11 

^6 


°C ’ 

35 


M 

f 

II 

a 

7 

11 

“C/VOLT 

36 

II 

11 

II 


II 

“C 

37 

II 

II 

21 


It 

“C//OLT 

38 

11 

1! 

It 


11 

°C ‘ 

39 

11 

tl 

It 

ao 

II 

°c/von 

40 

II 

n 

If 

^9 

U 

■ °C ' 

41 

II 

It 

11 ' 

^'lO 

11 

°C/VOLT 

42 

II 

11 

]| 

^10 

II 

°c 

43 

II 

It 

II 

. “n 

It 

°cA^olt 

44 

11 

II 

II 

Bll 

11 

°c 

45 

II 

11 

tl 

“12 

11 

°C/VOLT 

46 

II 

It 

II 

B12 

II 

°C 

47 

II 

II 

tl 

“l3 

II 

°C/VOLT 

48 

II 

tl 

II 

313 

H 

. °c 

49 

II 

II 

II 

“l4 

11 

“CAOLT 

50 

11 

11 

11 

B]4 

II 

°C 

51 

tl 

11 

n 

“l5 

It 

“CAOLT 

52 

tl 

tl 

11 

^15 

1] 

°C 

53 

It 

tl 

II 

“16 

II 

“cAolt 

54 


tl 

11 

®16 

II 

°c 

55 

II 

II 

It 

“17 

II 

°cAolt 

56 

II 

II 

11 

^7 

- 11 

°c 


57 EMPTY 


75 EMPTY 
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TABLE VIII SCATTEROMETER PARAMETERS 


POSITION 

ENTRY 


TYPE 

UNITS 

1 

FREQ 


R 

■ HZ 

2 

FEED (WG SW OR CIRC) 


A 

— 

3 

CALIBRATION PERIOD 

C 

R 

SEC 

4 

INT. TIME ^1 

1 

R 

II 

5 

INT. TIME *5^2 

2 

R 

II 

6 

INT. TIME ^3 

3 

R 

11 

7 

INT. TIME H 

4 

R 

II 

8 

INT, TIME ^5 

5 

R 

1! 

9 

INT. TIME #6 

6 

R 

II 

10 

'(SHORT SCAT INT TIME 

7 

R 

II 

11 

TRANSFER FUNCTION , 

^HT 

R 

— 

12 

tl II 

®VT 

R 

— 

13 

11 n 

^VHT 

R 

— 

14 

n 11 

®HVT 

R 

— 

15 

II II 

*^HR 

R . 

— 

16 

11 n 

®VR 

R 

— 

17 

II II 

^C 

R 

— 

18 

CAL ATTN CHAN 1 


R 

_ 

19 

tl 11 II 2 

Aj 

R 

— 

20 

II II n 2 , 

^3 

R 

— 

21 

" " '■ 4 


R 

— 

22 

FILTER LOWER LIMIT (NEG) 


R 

HZ 

23 

FILTER UPPER LIMIT (POS) 


R 

HZ 

24 

FILTER RESOLUTION 


R 

HZ 

25 

ANT GAIN H 

Tv 

R 

— 

26 

ANT GAIN V 

R 

— 

.27 

CHAN ^1 SATURATION LEVEL 


R 

VOLTS 

28 

29 

EMPTY 

MINIMUM SQUARE LAW LEVEL FOR 


R 

VOLTS 

30 

ALL CHANNELS 

ANGLE CONVERSION FACTOR 


R 

degAolt 
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TABLE VIII SCATTEROMETER PARAMETERS (continued) 


POSITION 

ENTRY 


TYPE 

UNITS 

31 

ANGLE CORRECTION FACTOR* 

R 

DEG. 

32 

EQUIVALENT BEAMWIDTH H 9eq 

R 

RADIANS 

33 

EQUIVALENT BEAMWIDTH V 9eq 

R 

II 

34 

EMPTY 



_ 

35 

ir 



— 

36 

n 



_ 

37 

n 



— 

38 

RELATIVE FILTER GAIN, 

f= 3000 

R 

_ 

39 

(1 

f= 2750 

R 


40 

ir 

f= 2500 

R 

— 

41 

n 

f= 2250 

R 

— 

42 

It 

f= 2000 

R 

— 

43 

It 

f= 1750 

R 


44 

II 

f= 1500 

R 


45 

II 

f= 1250 

R 

— 

46 

rt 

f= 1000 

R 

— 

47 

ti 

f= 750 

R 

— 

48 

ir 

f= 500 

R' 


49 

ir 

f= 250 

R 

— 

50 

t 

II 

f= 0 

R 

— 

51 

ir 

f= -250 

R 

— 

52 

n 

f= -500 

R 

■ — 

53 

II 

f= -750 

R 

— 

54 

II 

f=-1000 

R 

— 

55 

ir 

f=-1250 

R 

— 

56 

II 

f=-1500 

R 

— 

57 

II 

f=-1750 

R 

— 

58 

II 

f=-2000 

R 

— 

59 

It 

f=-2250 

R 

— 

60 

II 

f=-2500 

R 

— 


•A- 

Includes relative pitch between aircraft frame and antenna gimbal. 
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TABLE VIII SCATTEROMETER PARAMETERS (continued) 


POSITION 


ENTRY 

• 

TYPE 

UNITS 

61 

RELATIVE FILTER GAIN, 

f=-2750 

R 

— 

62 

i < 


f=-3000 

■r 

— 

63 


If 

f=-3250 

R 

— 

64 


11 

f=-3500 

R 

— 

65 


U 

f=-3750 

R 

— 

66 


It 

o 

o 

0 

1 

R 

— 

67 


n 

f=“4250 

R 

— 

68 


ri 

f=-4500 

R 

— 

69 


II 

f=-4750 

R 

— 

70 


n 

f=”5000 

R 

— 

71 


n 

f=-5250 

R 

— 

72 


ti 

f=-5500 

R 

— 

73 


li 

f=-5750 

R 

— 

74 


ir 

f=“6000 

R 

— 

75 


ti 

f=-6250 

R 

— 

76 


ij 

f=-6500 

R 

— 

77 


II 

f=-6750 

R 


78 


M 

f=-7000 

■ R 

_ 

79 


ir 

f=-7250 

R 

— 

80 


It 

f=-7500 

R 

— 

■81 


ti 

f=-7750 

R 

~ 

82 


II 

o 

o 

o 

00 

1 

R 

— 

83 

'EMPTY 





84 

ir 





85 

IS 
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entries. The ^hermisfor conversion factors in RPAR(23) through RPAR(56) convert the 
thermistor voltages in DATA{17) through DATA(33) to temperature in centigrade 
(see Table II). A sampled representation of the scatterometer doppler filter frequency 
characteristic appears in SPAR(34) through SPAR(85). The sample frequency interval 
is designated In SPAR(24) and the frequency domain of the doppler filter is specified 
in SPAR(22) and SPAR(23), 

3. Peripheral and Memory Requirements 

The peripheral requirements are illustrated in Figure 9. To load the program 
required approximately 10,000 decimal words of memory on a CDC-6600 to include 
buffer space and system routines. 

E . Program TDUPE 

Program TDUPE scales and integerizes output files prepared by program 
CONVERT and transfers the result to another tape, • The resulting tape serves as a 
master tape which may be duplicated and sent to the investigator for use on another 
type of computing machine. 

The scaling and transferring of an output tape is implemented by an external- con- 
trol card appended to the source or object deck. The card simply states the files by mission 
and flight which are to be transferred. The control cord requires a 2110 format. The 
input files must be located on tape unit 01 . The specified file types are transferred 
to a tape on unit 02. An activity report, stating the status on the transfers, is 
generated. The program requires approximately 8000 decimal words on a CDC-6600. 

The peripheral requirements are illustrated in Figure 10. 

F .«• Special Programs 

Several engineering application programs are provided to compute anc{/or 
verify certain antenna parameters employed in the instrument parameters. Program 
WIDTH computes the equivalent pencil beamwidth, an important parameter to the 
inversion of the scatterometer data (SPAR(32) and SPAR(33)). The program GAIN 
computes the antenna gains SPAR(25) and SPAR(26) and the beam solid angle. 
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FIGURE 9. peripheral REQUIREMENTS FOR PROGRAM ICHAR 




FIGURE 10. PERIPHERAL REQUIREMENTS FOR PROGRAM TDUPE 



Program WIDTH is seif sufficient when antenna pattern information is 
provided in the FUNCTION subroutine FUN . Presently functional representations 
of the main beam out to 6 degrees from boresight are provided. This routine 
may be replaced by a sampled version of the main beam. However, interpolotion 
between sampled values must be provided since a double integration routine requires 
values dictated by the Gaussian Legendre quadrature technique. 

Program GAIN requires both polarized and cross polarized power patterns. 
Functional representations of the polarized pattern Is now given in FUN! and the 
cross polarized pattern in FUN2. These may be replaced with sampled valued 
representation, if desired. Again, interpolation will be required. 

For more specific information see Appendix E. 
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APPENDICES 



APPENDIX A 


CONVERSION PROGRAM 


I, INTRODUCTION 


The conversion program (CONVERT) applies many reduclion algorllhms fo 
the raw RADSCAT dala lo provide ihe Invesligator wihh a se^ of comprehensive dafo 
producfs wi^h which fo inferpref ihe sca^^ering and emission characferisfics of fhe 
scene represented by the data. CONVERT consists of the Control Section and the 
following subroutines; 


1. IREAD 

2. READl 

3. WRITE2 

4. ABORT 

5. SEARCH 

6. CRUNCH 


7. CALIB 

8. RCONV 

9. ANGLE 

10. DOPCHK 

11. SCONV 

12. THERMO 


The theory and design of these routines are treated in this appendix. 


II. GENERAL OVERVIEW 


The Control Section reads and interprets external directives. On the basis 
of the directive, the appropriate files are retrieved and presented to the processing 
subroutines. Subroutine SEARCH executes entry into the data records of the file in 
search of the appropriate calibration records. When found, the calibration 
information is extracted and averaged in subroutine CALIB. When sufficient 
calibration data cannot be found, the file is bypassed. The actual processing is 
performed by subroutine CRUNCH which applies various algorithms and calibration 
information to reduce the data. Many of the algorithms are embodied In sub- 
routines, bearing descriptive names, which are called from CRUNCH. 

Subroutine ANGLE computes the incident and cross-track angles using 
aircraft attitude and look angle information. Subroutine THERMO checks Internal 
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receiver |■empera^tJres . The anfenna femperafure is computed in RCONV using fhe 
transfer funcHon of the RADSCAT radiometer, the temperatures developed in 
THERMO and RAD calibration parameters. When scatterometer data is present, 
the maximum on-scale measurement is selected from the four SCAT channels with 
an in~llne routine In CRUNCH . Using aircraft ground speed, DOPCHK computes 
the doppler shift In the scatterometer signal. The doppler filter gain Is chosen 
accordingly. The scatterometer transfer function is applied in subroutine SCONV 
to compute the scatterometer normalized input power. An in“line routine in 
CRUNCH applies antenna inversion parameters to the input power to compute the 
normalized scattering coefficient. 

When additional calibration records are encountered, the old calibration 
information is replaced by the new calibration information by subroutine CALIB. 

Subroutine IREAD and READ! perform the file read functions, whereas 
subroutine WRITE2 performs the write functions. 

Subroutine ABORT is employed by the CONTROL Section to abort the job 
when an abnormal condition is encountered. 


III. CONTROL SECTION 

A. Theory and Design 

The control section was developed to give the user a large amount of 
flexibility In the choice of the instruments characteristics and the files to be 
processed. As many or as few files of a particular type can be processed. Some files 
can be treated with one set of instrument characteristics and others with another. 
These features permit the user to (1) test a set instrument parameters, (2) group the 
files by type, (3) process different types of files with different instrument parameters, 
(4) etc. The flexibility in the choice of parameters was incorporated with the full 
realization that the Instrument would undergo modifications. 

Control of the program is executed by the choice of appropriate directives 
entered as data cards. The types of controls have already been discussed in Section II 
and III of the text and will not be treated here again. The method by which control 
cards are read' and files are retrieved is shown in the logic diagram of Figure A-1. 
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FIGURE A-1. 




(c) DESCRIPTIVE LOGIC DIAGRAM FOR 
THE CONTROL SECTION 
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The numbers on the connecting points coincide with the statement numbers of the 
text shown in Figure A-2. 

An instruction is read from the card reader. If it is POST;, the files on an old 
output tape are bypassed in statements between 39 and 1 of the diagram. Another 
instruction Is then read. The instruction is echo checked, file parameters are initialized 
and in the subsequent blocks of Figure A"l, the processing type is identified. 

Thereafter the specified instrument characteristic is read from unit 04 or 05. If the 
last set of characteristics on a file code 04 is desired and its file number is not known 
(ICU=999), the last file number is determined in statements between 6 and 5. If, on 
the other hand, they are to be read from cards (ICU=0), the appropriate action is taken 
and control is advanced to statement 15. If a specific file is designated, an appropriate 
number of files is bypassed until tape unit 04 is positioned on the correct file (see steps 
between statements 10 and 13). 

If the set is to be punched or printed, the- appropriate action is taken in steps 
preceeding statement 14. If NFP>1, NFP-1 subsequent sets are aiso punched or printed. 
(See steps between 112 and 13). When completed, the next control card is read. 

When the control card is a processing directive, the instrument characteristic 
tape will have been positioned to the correct file (set of parameters) or the set will 
have been entered from cards. Steps following 15 search for the appropriate data file. 
The search is governed by the type of directive. Bypassed files are cataloged in steps 
between 25 and 22. If the desired file is In the catalog and the current file Is strictly 
in the catalog, the appropriate number of files are bypassed in steps between 22 and 24. 
If the file is beyond the catalog, the catalog’ is updated as the program advances the 
tape In search of the file. 

Once found, the file label is lis'ted and the appropriate subset of instrument 
parameters is withdrawn from tape if they weren't entered from the reader. SEARCH 
is then called to find the calibration data. SEARCH in turn calls CRUNCH to 
process the measurement records. Once completed control returns to the Control 
Section where additional files are sought or another directive is read. 

B. Program Listing 

The Control Section Listing is shown in Figures A-2(a) through A-2(e). The 
variables are defined in the listing. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C ■ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c_ 

c 

c 

c 

c 

c 

c 

■ c‘ 
c 
c 
c 
c 
c 
c 
c 


t’KOoHAt) CLN\^tPT I IhPUr.UUTPUT.PUNCH. lA Pt 1 » T AP E<:i T APfc4, TAPEfl » 
TAPCi-=INPJl . IAPf6“0ljrpUTiIAPL4 J-PUNCm 
CCNVLKSICN CLNIBOL SECT 1 OH 

‘ this PRCORAH has PREPAkEO UY 

jChN P. CLAASSfcN 
CUES E. ELLIOTT 

,UMVEHSITY OF KANSAS CENTER FOR RESEARCH 


PROGRAM UfSCKlPTION. 

FILE CCOES USED. 

01 — INPLT RAh DATA. 

02 — LUTPUT PROCESSED DATA. 

04 -- INPUT INSTRUMENT CHARACTERISTICS FILE. 

05 — SltSTEH INPUT - DIRECTIVE CARDS. 

Oo — STSTEH OUTPUT - PRINT LINES. 

08 — SCRATCH file - USED BY CCHPLTATION SECTIGN. 

43 — SYSTEM CLTPUT - PUNCHED .CAROS. 

VARIABLE NAMES LSEO. 

LBLlltL YLAUEL TABLEa. TABLE CONTAINING 

LABEL INFORHATICN IHISSICN. FLIGHT, LINE, RUNl fCR ALL 
LATA FILES bypassed SO FAR. .. 


NWPST 

NCF 

NFF 

H ST Art 

_HJ L*. 

EQM 

lEOR 

IPTYPE 

JPTYPE 

__ 

ICU 


NCC 


aNER position.* array containing label INFO (H,F,L,R1 
OF next file to be PROCESSED. 


ANUHBEF of CURRENT FILS.* 
FILE CODE 01. 


POSITION POINTER FOR 


anumbcr of final 

*LBLT8L* TABLE. 


file.* index cf last entry in 


STARTING POINT FOR SEARCH IN *LBLTEL» TABLE. 

'relatTve pcsitiion on file CCCE 01. 


1 INPUT FILE TOJ.se 

^ ■ I 


CCCOCOOT 

OOOOOOOU 

oooooooy 

COOOOOlO 

OCOOOOll 

COOOODI2 

00000013 

00000014 

COOOOOIS 

CCOOOOI6 

COOOOOI7 

C000002B 

00000029 

00000030 

CCC0003I 

0CO00032 

C0000033 

00000034 

00000035 

C0000036 

00000037 

COOC0038 

C0000039 

CCOCC040 

C000004I 

CCOC0042 

C0000043 

00000044 

CC000045 

00000046 

C0000047 

C000004B 

C0000049 

00000050 

CD000051 

CC000052 

00000053 

CC000054 

00000055 

CC00C056 

CC000057 

CC00C059 


INCICATING CCOCOOoO 
CC00006I 


NLC 


RPAK 

SPAR 


PAR 


NFP 


NFAi* 


*ENO OF REEL INDICATOR.* LCGICAL VARIABLE, SET TRUE 
rFEN end of REEL ON FILE CODE 01 IS DETECTED, 

*LbLTiiL* ITABLEUS CCPPLETE. ’ 

CCCCC062 
00000063 
CG0000O4 
C0000065 
C0000066 
CCCOG067 
00000068 
CC00Q069 
CQ000070 
00000071 

oroooo72 
CC000073 
C000C074 
CC000075 
CCG0CCI6 
CC000077 

*N0H8ER of current CHARACTERISTICS*. REFLECTS PCSITICN CCC00C078 
FILE code 04. CCCC0079 

00000080 
CCCCC081 
00000082 
C0000063 
C0000084 
C0000085 
C0Q000B6 
CC00008T 

ccoccc&e 

00000089 
CC0CC090 
00000091 
CQ00C092 
C00CC093 
00000094 
00000095 
00000096 
CC000097 
CC00009B 
CC0C0099 
COOOOIOO 
CCOOOIOI 


INTEGER ARRAY CONTAINING kQROS *ENC OF REEL*. 

INTEGER *PR0CESSING TYPE* ARRAY, CCNTAINING RECOGNlCEC 

PROCESSING Types. 

INPUT DESIRED PROCESSING TYPE. 

'TND'e'x"6F‘>0S*ITIUN OF *JPTYPE* IN *IPTYPE*. 

*INSTR0M6NT CHARACTERISTICS TC USE*. NUMBER OF THE I.C. 
SET TO BE USED. REFLECTS THE PCSITION 0<= THE SET CN 
■ FILE CCCE 04. 

0 INOICAIES I.C. TD BE READ FRCH CARDS IFC*05), 

999 MEANS USE TFE LAST SET ON FILE 04. 


*NUMBEP OF LAST CFARACTERISTIC*. LSEC WHEN ICL=999. 

NLC = 0 INDICATES THIS VALUE IS YET LNKNORN. 

»raeicmeier paraketers*. 

*slatercmeter parameters*. 

*PARAHETERS* ARRAY CONTAINING DP TP 4 SETS OF RACSCAI 
INSTR PARAMETERS READ FRCF FILE CCDE 04 OR OS. RAC 
PARAMEIEHS compose the first 75 entries of par AND SCAT 
parameters OCCUPY THE NEXT 85 ENTRIES IN EACH SUBSET. 
THE APPROPRIATE SUBSET OF PARAMETERS IS SELECTED ON THE 
basis cf FPBDUENCY AND FEED, 

*NUHBIR of files TO PROCESS* LSING THE PROCESSING TYPE 
given by *JPTYPE* and DESCBIPTICN in *NrfPST». 

NpP IS ignored if JTYPL ■= SPECIF. IF ICO = 0, 

NTP IS SET TO 1. 


*NU.MBEK OF FILES ALREADY PROCESSED.* 
TO CHECK FOR TERMINATE CONDITICN. 


COMPARED TO *NFP* 


*^FIGURE A-2. (a) LISTING OF CONTROL SECTIGN 
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ccooo:} 

CC0O03 

COQ003 

CconOj 

(000C3 

C00003 

C00003 

C00003 

C00003 

C00003 

CC0003 

CC00C7 

COOOZ5 

COOOZ7 

000036 

CCOO'tO 

COOO<i2 

C00044 

COOO'tS 

000050 


C00053 
C00C71 
000071 
COOOl'V 
000 112 


C00112 

000113 

C00115 

000115 

000117 

000123 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

c 

c 

(. 

0 

c 

c 

c 

c 


LilPHUU XLI.Nfc HrllMlNO PCDDLUS#, OSEo U1 KEUUC6 UOTHUT 
PKCP tCPKiIPTUh iCtnuw. 0 IS AS 1. 


OEKR 


OF ERKOKSk. NOHOER Ot TIMES LOHRorATliJN SECIICN 
KEIURNtG ERROR CLNOITIUN. hCRR « A CAUSES ABORT. 


ERRFL6 etc FILE ERROR FLAO. 


PROCESS IMG 
TYPE 
All 


DIRECTIVE CARD CSAGE 
■pEANItlG “ 


START FRCP FIRST OF INPUT FILE ARC PRCCESS WlThOLT 
REGARD TO LABEL INFURMATICN. 


HISFLT CFECK LABEL ONLY FOR CORRECT PISSICN ANO FLIGHT. 
FLTLIN ChECiL FOR CORRECT MISSION. FLIGHT, AND FLIGHT LINE. 
specif CHECK FCR MISSICN, FLIGHT, FLT- LINE, ANO RUN. 

PRINT PRINT CUT VALUES 0£. INSTRUH ENT .CHARACTER 1ST ICS . 


CCOOU102 

CCCCCIOJ 

ooorolOA 

CC000105 

C00001C6 


OOOOOlC/ 

COOOOlOB 

cooooiov 

COOOOllO 
CCOOQlll 
00000112 
C0000113 
C0000114 
CC000115 
00000116 
COOOOllT 
C( 000118 
COOOOllV 
C0000120 
CC000121 
CO0001Z2 
00000123 


c 

PUNCH 

PPINF 

PUNCH 

l.C 

• 4S « 




COOCC124 

c 










CC000125 

c 

PU:>T 

HAf 

BE LSEL AS 

THE 

FIRST PRCCESSING TYPE. IF SO, 

C000013B 

c. 


Tllb 

PPCCRAM hILL SPACE 

DOhN 

TFE CUTPUr TAPE SEARCHING 

CC00C139 

c 


roR 

THE ENU OF 

REEL 

. POST : 

IS ILLEGAL EXCEPT ON THE 

00000140 

c 

f 


FIRST CIRECTIVE 

CARU. . 




CC000141 

c 


CIRECTIVE 

CARD 

FORMAT* 




C0000126 

c 










00000127 

c 

JPT.YPe 


MISS 

FLT LINE 

RUN 

LNPN 

ICJ 

NFP 

00000128 

c 

A6 

4X 


IS 

IS 

IS 

- 15 

.„I5 

- - 15 

C0000129 

c 










C0000130 

c 

ALL 


♦ 

» 

4 

4 

1 

3 

6 

ccocom 

c 

^il3FLT 


2C? 

12 

4 

4 

5 

3 

2 

C0000132 

c 

FLTLIW 


207 

12 

I 

4 

10 

3 

5 

CC000133 

c 

SPECiF 


2C7 

12 

i 

2 

1 

3 

1 

C0000134 

c 

PUNCH 



« 

4 

4 

4 

4 

6 


c 

PRINT 



4i 

4 

4 

4 

S 

7 


6 

POST 


4 

4 

4 

4 

* 

4 

4 



» INDICATES FIELD IS IGNUPEE EV PROGRAM. 


DI.lENSlcN NbPSTUJ. LBLTBLI 4,300 ) IECRI2) , RPARI75J, SPAR185), 
PARI160,4). IPTVPE16), IDAT131 
UlHFNSICN LdLIVI 

E«u I VALENCE (LBL18J ,FREU I , ( LBLI 9} , F EEC) 

EUJ 1 VALENCE I IDAT 11 ) ■ 101), ( ICATl 2 ) , 102 ) , ( 1 CA t (3) , 1 D3) 

LOGICAL EUR), ERRFLC 

LATA IEGR /IOfENC Of REEilHL/ 

U«TA 1PTYPE/6HF LTLI N,6HHI SFLT,6HSPEC I F, 3HALL , 5HPR1NT, ShPUNCH/ 
DATA 1PUST/4MP0ST/ 

LATA ECRl/. FALSE./, NL*C,NCC/0, 1/ , NCF ,NFF/ 1, 0/, NERR/0/ 

PtrfIND 1 
REMIND 4 

KEADIS.SCC) JPIYPE, NhPST, LNPHQD, ICU, NFP 

iF( jptypl.ne.ipcst) go to 18 

REA012) LBI.LB2 

ir iLul.NE. lEOiUl) I GO TO 40 

IFIHJ2.HE.1ECR12)) GO TC 40 

backspace 2 

GC TO 1 . _ 

REAL! 2) 

IFUUF.2) 39, 40 

READ lIRECTIVES 


hhAT MANNER TO PROCESS 

H«1 

HST,ART«1 
NFAP»0 
nu t. N>1,6 

IFI JPlYPt.E J. IPTYPE(N) ) GC TO 3 
«RI n (6.601) 




CC0n0l33 

CCC00136 

CCCOC137 

C00001A2 

CCCOCOld 

00000019 

CC0OCO2O 

C000002L 

CC000022 

C0000023 

CG000024 

CC000025 

00000026 

CCCC0027 

00000143 

C0000144 

OOOOOL45 

C0000146 

CG000147 

C0000148 

00000149 

C00C0150 

CfOCC13l 

000001^2 

CC000133 

C000013A 

C00001!>!> 

00000136 


1 

REAOIS.SnO) JPTYPE, 

NWPST, 

LNPMOD, ICU, 

NFP 

00000157 

:>oo 

format (A t. 4X, 7)5) 




C0000158 

m 

IFlFuF,5l 37,41 




C0000159 

41 

n'RITE(6.u00) JFTYPE, 

NhPST, 

LNPHOO. ICU, 

NFP 

C0000160 

600 

FOKHAKIHI ,//10X,17H 

CCNTRCL 

CARO IS , 


CC000161 


♦ A6,4X,VI5) 




CC000162 


^ ^ T 


6. 


COOOOltB 


C0000164 

CC000165 

COOOOlOO 

CG000167 

CCOCOloO 

G0000169 

CCOC017O 

C(000l7l 


FIGURE A-2. (b) LISTING OF CONTROL SECTION 
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OOOU6 

oOL 

1 im •lAt HOa ,30liL (66C06.4I 6fcl/ ('Ko4CiS16v (043/T( 

cfroci/2 

0001<!6 


Lj4LL 46CH 1 {261.H 1 

rooooi 73 

CQ0130 

3 

IMNlf'.or. 11 06 tc 59 

crooo 174 

000133 


ndl T 6 lo * 6 1 0 } NIP 


COOi<iO 

010 

^0■<MAT I // 1 {rXt »6f P LT l»i5x,*KPP » ■*.15//) 


COOl^O 


LALL AbLMT ( JriM-P 1 



C 


cccrci 15 


C 

Mho INStXUHENT CriAA. 

00000176 


c 


conoo(77 

000142 


IPl ICL.NE.OJ CC TC 4 

COOOO 1 7B 


C 

READ FRUM CARL'S 

CC000179 

C00143 


[LSET-O 

CCOCOIBO 

0001,44 


read (5, 4 30 3) ADKBERi ID l, IC2f 1D3 iRPAR{ 11 (RPAR12J ,RPAK(31 

OOOOOIOI 

OOOlbb 


■IHEGP»5) 12, (0 

CCCC01U2 

000171 

oO 

wRl'lEib’, 60 31 lC$ET,lUl,ll)2,lU3,HPARaj,RPARl2J, RPAR13) 

OOoOO UI3 

C00213 


UO 6l J»4,75,7j 

CCOOO lt4 

000215 


L=Jt6 

COOOOlbS 

0002 IT 


rF(LlCT.75» L=75 

COOOO 186 

000222 


REA015,4304) CPPAh (K1 ,K-J,L) 

cconoiu7 

000235 


IKECF.S) 12, 61 

COOOO lUB 

C00240 

61 

nRITH6,6a4l(PPAR(KI,K2j,L) 

CC000189 

000255 


REAiH 5,4 3031 MPBER , IDl , IC2, t03, SPAR (1) ,SP/R (2 1 , SPAR t 3) 

COOOO 190 

C00277 


IFIboF,5l 12, 62 

6CO0C151 

C00302 

o2 

nRl TLl6,t031 ICitT,l01,IC2,I03,SPAR(ll ,SPAR12J ,SPARt3 1 

COOOO IS2 

000324 


DO 63 J==4,B5, 7 

CC000193 

C00320 


L = J<-6 

CC0CC1S4 

C003i0 


IKL.Gl.t5J l = 65 

CCOO0145 

0003J3 


REAC15,4JC4) tSPARlKl ,R,= J,U 

CCP00196 

00034b 


1F(EI1F,5I 12. 63 

00000157 

CQ0351 

63 

nKl TL (6,6041 ( SPAR (RJ ,KsJ,Ll 

CCOOO 198 

C00366 


NFI> = i 

COOOO 199 

C00367 • 


GO TO 15 

OCOC02OO 

CG0370 


IM1LU-6CCJ 5.13,6 

00000201 

000373 

5 

remind 4 

cr. 0002 02 

C00375 


ilCL.* 1 

C0000203 

000376 

6 

1F( lLU.Nfc.9991 JC TU 10 

CC000204 


c 

FIND LAST SET 0,9 FILE. 

CCC002C5 

CQ0400 


IF1MC.Eu.0J go to 7 , 

00000206 

C00401 

t2 

ILu=NLC 

CCC002C7 

000403 


LU TU 4 

OOOC02CB 

C00403 

( 

NLC»NCC-2 

00000209 

C004C5 

a 

kEAD(41 

C(j000210 

000410 


NLC=HLC+i 

COO00211 

000412 


lF(cUF,41 45, 0 

r00('O212 

000415 


1CU=NLC 

00000213 

000417 


OU TJ 5 

COOOC2 14 


c 

POSITION TO THE SET 

conoo4i5 

C00417 

10 

.\CL=iTCCK 

CCC0021b 

000421 


uo 1 i 1 = NC 1 1 1 Cl 

CC000217 

000422 


HEAD14) 

.CCCCC21B 

COC425 


1(-(E0F,4I 12,11 

00000219 

OC0430 

11 

CCi'NTlNUfc 

CC000220 

000433 


uG TO 13 • 


000433 

12 

nRlT£t6,(,02 1 . , 

00000221 

000437 

6 02 

~ format (10X,J4rit,C SUCH lASTRUHENT CHARACTERISTICS! 

CC000222 

000437 


CALL AiiORT(2HICI • 

00000223 


C 

READ THEM 

CC000224 

000441 

13 

REA0(4J tCSET.lLl ,102,103, PAR 

C0000225 

000456 


lHcOF,4l 12,112 

CC00C226 

000461 

_ , 112 

NCC=1CU+1 

GC00Q227 

000463 


1KN.LT.5J GO TO 15 

CC000228 


C 


CC000229 


c 

WRITE [AND PUNChl INST. CHAR. 

00000230 


c 


C COO 02 31 

C00466 


DO 55 1=1,4 

00000232 

000467 


mRI T u (6.6031 ICU, 101, 102,103, PAR (1,11 .PARC 2,1 1, PAR (3, 11 

CC000233 

000515 

603 

FORHATlA/4 SET*. 15,* DATE* ,31 5, ElS.4 ,5X, AlO , E15.4 1 

C0000234 

C00515 


UU 54 J=4,75,7 

C0000235 

C00517 


L = J*6 

00000236 

000521 


IF(L.GT.75) L=75 

CC000237 

C00524 


MKI TE( 0,604 1 ( FAR (K, n ,N=J,L1 

CC00023B 

000543 

61)^ 

. FUR71AT17E15.41 

00000239 

000543 


mRI T£ ( 6,6031 ICU.ICl, 102, 1D3,PAR(76, 1 1, FAR I 7 7, 11 ,PAR( 7B.1 1 

CC00C24O 

C00572 


DO 5? J= 75,loC, 1 

00000241 

C00574 


L=J *6 

C0000242 

0005 76 


IHL.GT.1601 L=150 

C0000243 

000601 

5b 

riRlTC[6.6041(FARlK.lJ ,K»J,LJ 

COOD0244 

C00622. 


. 1I’(..ME,6! go .1 C_i4_ ... . 

C0000245 

000624 


UO 57 1=1,4 

00000246 

000626 


LL=1 

CC000247 

C00527 


MRI 11(43,43031 1CU,I01,I02,1U3,PAR (1, 11 ,PAR(2,1 1 ,PAR(3, 11 , 

C000024B 



* ICO.I.LL 

C0000249 

000664 

^303 

FORHAT(4l5,lOX,C10.3,A10,E10.3rl2X, 12,11,151 

CC0Q025O 

000664 


...DO 56 J=4 , ?5i 7 

00000251 

000666 


LL=LLF1 

CC000252 

C00670 


L=J*o 

00000253 


ORIGINAIj PAGE IS 
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(00672 

bb 

i,Kl TLl^j.9 309) (FARlR.l) ,X*J.LI .1 Cu . 

I, U 

C00715 


1 tlRHATl/LlO.i.ZXi 12,11.151 


000715 


LL’LLFl 


CC0717 


«RI 11193,93031 ICU,Ii;AT,PAKt76, 1 

1 , PAR (7.7 


* 

ICG. i ,Ll 


000797 


01) 57 J°79,16C,7 


000731 


LL»LL+1 


000753 


L" 3*6 


C00755 

57 

mRI IE ( 93,9309) (PAR IK, It ,K*J,L),lCUt 

I.LL 

C01002 


ICU>ICU»1 


C01009 


NFAP-NFAF+1 


C01005 


IHnfap.lt.nfpi go to 13 


COIOC7 


GO TO 1 



C 

c 

FINO Tt-E DATA RECORO. 



001010 

COiOlZ 


001016 

C01017 

COIOZZ 

C010Z5 

C010Z7 

COiOlZ 

001033 


COllOl 

001105 

C01105 

001107 

C0U15 


15 IH.N,EO.<t) GO 1C Z3 

IH«5IAKT.Cr.MF I 00 TO 17 

0 SIE IF »e/vE BEEN BT IT ALREAOV. 

1)0 16 N=>RSTAjIT.NEF 

IFCNhPST in.NE. LBL1BLU >M) ) GOTO 16 
IKhnPST 121 .mE.LBLTBUZ.HI ) GOTO 16 
JF( N.FU.Z) >.0 TC 23 
lr(OnP5TlAl.rg£.LBLTBl(3tN)) OO Tj 16 
IFIN.Fo.l) GO TC 23 

IFlNnPST(4).NE.UiLTBU5,K)l 00 TO- 16 


C01036 


, l>0 TO 23 



C01036 

16 




001041 

17 

1F(EURU CO TO 19 - 




c 

ggESS IT»S further OOdN 

001043 

c 

H=NFF+1 



0010^5 


IFINLF.GT.NFFJ CCTO 

24 


0Q1050 

20 

H1=K-1 



oolosa 


00 27 1=NCF,H1 



001054 

Zb 

READ! U 



C0LQ57 


IFlFOF.l) 27, 26 



001062 

n 

continue 



C0I065 


NCF-K 



001066 

24 

REAOIl) L8L 



001073 


IFI LBLID.I.E. lECR(l) 

1 GO 

TO 25 

C01075 


lFILdU2I.NE.IECR(21 

] GC 

TO 25 

colon 


eORI=.TKu£. 



COilOO 

19 

IFINFAP.NE.OI GC TO 

35 



o06 


35 
6 07 


«F ITE16,6C6) 

FORMAT ( 10X,15FNC‘SUCh' CATA' RECORD) 

CAU. ABORT tZHPFl 
hRlTElotCC?) NFAP 

FORMAT (lOXf 13 »32Fi FILES FOUND. CAN FINC NC HCKE.-I 


CD12B2 


oolzo; 

C01203 

C01Z05 

001215 

001221 

001223 


60B FURHAT(lHl,51X,I7hCCNVEBSlCN REPORT. 

» //lOX.dUHISBlCN-, 13.0H FLIGHT-. 13, 

* 7H 0ATE-.2I3.I5.2X.11H FLT LINE-. IX, 15. 6H 
_ _ 13H FREi.UENCY-,F5.1,7H F£EC-,2X,A6» 

* ’ 6F 1PAR=, I2.3IJ////I 
IF| ICU.CL.O) GO TO 33 

C GET CORRECT INST, CHAR. 

DO zy K=l,5 

IF IAoSIFKEC-P 7R(‘ .KI/l.O EOl.GT.l.Ol CO TO 2V 
IFll LLD.C w.PAR(2,k 1 ) GO TO 30 
29 tCNTINuE 

hRIIE(6.cC5) FREu, FEEU 


C00002I5 

00000255 

CC0CC256 

C0OGO257 

conoozsa 

C0II00259 

00000260 

00000261 

COOOOZeZ 

00000263 

(0000265 

00000265 

00000266 

C0000267 

C000026B 

CC000269 

00000270 

CC000271 

C0000272 

CCOOC273 

C0000275 

COOOQ275 

C0000276 

CC000277 

CCC0027B 

00000279 

CC000280 

CCCOCZBl 

00000262 

C00002B3 

CC0002B4 

CC0002fc5 

000002B6 

C0000287 

C000028S 

00000289 

00000290 

00000291 

00000292 

CC000293 

CC000299 

CCCCC2S5 

C0000296 

CC0QQ297 

C0000298 

CC000299 

COL00300 

CO0UQ301 

C0000302 

C0000303 

C0000309 

OC000305 


001115 


GO TO 1 





00000306 

C0U16 


LBLTdL(l,NCF)=LELUl 




(CCCC3C7 

C0I121 


LdLTBL(2,NCF)=LaL(21 




000003C8 

C01124 


HiLTbU3.NCF)=Laci6 1 




00000309 

001126 


LBLTdLI9,NCFl=LaL(7) , 




CC000310 

COlUO 


HFF=MAxrHNFF,NtF 1 




CC0003U 

001133 


1FIN.EG.9I 

oil TC 36 




CC000312 

001135 


IFtNdPST (1 

) .NE.LbLTdL(L,HI 1 

GO 

TO 


CC00C313 

001140 


‘iFdJhPSTiZ 

).NE.LdLTBL(2,K)) 

00 

TO 

22 

00000319 

C0I143 


IF(N.Eu.2) 

GC TC 36 




C0000315 

001144 


IH NWPST13 

).NE.LBLTBL(3,H)) 

GG 

TO 

22 

00000316 

001147 


iFlN.Cu.l) 

gC TC 36 




CC000317 

C01150 


1F(NhPST(9 

).NE.LBLTBLI9,Hn 

GO 

TO 

CM 

00000316 

CQ1153 


GO TO 36 





C0000319 

00U53 

22 

H=M+1 





C0000320 

001155 


GC TO 20 





CCCCC321 

001155 

23 

IFIH-NCF) 

28. 29, 20 




00000322 

001160 

20 

remind 1 





C0000323 

001162 


NCF*^1 





(CC0C329 

001163 „ 


eu TO 23 





00000325 


C 

FOUND 

itL 

ITj^S This ONE. 

0CO00326 

001164 

p6 

V>RITE(o,60d) LBL, ICSET, IDl, 

102 

» 103 

CC000327 


RUN-,13, 


00000326 

COOC0330 
00000331 
00000332 
00000333 
cr.ooc3j5 
00000335 
C0000336 
CC00033 T 
Cl00033d 
00000339 


FIGURE A-2. (d) LISTING OF CONTROL SECTION 
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C012JJ 

605 

lUKJlAT (10;c,171-KL SUCH IKKo-f EEU,F10.5, lX.AlO) 

001233 


CALL A()0KT{2Hfi J 

001235 

30 

DO 31 1=1,75 

001237 

31 

«PAR( I l=HAK( I > 

001246 


DU 32 l=l#fc5 

001247 

32 

SPAlt(U=f'AK(l+75fK» • 


C 

CALL COMPUTATION SECTION 

C01256 

33 

tKKEL6 =. false . 


C 

LABEL OUTPUT TAPE 

00125/ 


riKlTE(2) LCL 

001264 


«pnE(2J ICStT, lOl, ID#:, ID3, RPAR, SPAR 

C013U3 


CALL iEARailSPAR.KPAR.LAPHUO.ERKFLCJ 

C013C6 


CALL ENU(-2 

C01307 


NCE=NCF+1 

0013 11 


lE(FhRPLG) NERP=NERR+l 

0013 13 


IKHbkR.CT.3) CALL AUCRTIZHER) 

001317 


M=M + L 

001321 


MbTArsT=H 


c 

MURE FILES TU PROCESS 

001322 


lElN.Eu.S) oO TC 1 

001324 ' 


NFAP=-NI-AP+1 

001325 


IF1NFAP.CE.NFF> CG TO 1 

001327 

r- 

GO TU l6 


L# 

c 

terminations. 

C01330 

3 1 

NKITH2I lECR. (1, 1=3, 9) 

001342 


..RI IF(o, /00> 

001346 

7 00 

FORMAT ani,///lUf 10 I1H4 JjlUHNORMAL TERMINAT ION, 

001346 


STOP 

C01350 


cNC 


roooo 

cooooa^i 

OOQWOi'ki 
OOOOOJ^j 
CC0003^A 
0( COOS'#!? 
(1000346 
00000347 
COOOOJ4d 
tt0O0349 
COOOCi!>0 
0000036 1 
00000 352 
CCC0C353 
00000354 
00000355 
C0000356 
C0000357 
00000358 
COOOOJ59 
00000360 
C0000361 
00000342 
CC0G0363 
00000364 
CC0O0365 


C0000366 

000003o7 


FIGURE A-2. (e) LISTING OF CONTROL SECTION 
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C. Peripheral Requirements 


UNIT 01 

RADSCAT Raw Data Tape 

UNIT 02 

Output Tape 

UNIT 04 

Characteristics Tape 

UNIT 05 

Card Reader 

UNIT 06 

Printer 

UNIT 08 

Scratch Unit (Perfeiably Disc) 

UNIT 43 

Punch 


IV. SUBROUTINES 

A. Subrou^^ne Search 

1 . Theory and Design 

This roufine serves as an enfry and exit point for the selected file. One of 
its primary purposes is to find calibration data of both kinds. It anticipates that 
both NORMAL and BASELINE calibrations are in the opening records of the file. 

If this is the case the processing of the remainder of the file is given to CRUNCH. 
Otherwise, a search for the missing calibrations is conducted deeper into the file. 

The bypassed files are stored on a scratch unit (08) for subsequent processing. When 
both types of calibrations are found, the bypassed records on file code 08 are processed 
by CRUNCH and the remainder of the records on file code 01 are also. 

If the return from CRUNCH is without an error flag, control is simply returned 
to the Control Section. If not, an error message Is given and the remaind er of the file 
bypassed before returning. 

See Figure A“3 for a descriptive logic diagram of SEARCH. 

2. Program Listing and Variables 

The listing for SEARCH is shown in Figure A-4(a) and A'"4{b). The 
definitions of the variables are given In Table A"1. The entries in DATA 
were defined in Section III B. 

Table A~1 Definition of Variable Used In SEARCH 

Vector Containing Raw RADSCAT DATA 
Integer Equivalent of DATA 
SCAT Calibration Vector (Chonnels 1-4) 
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DATA 

IDATA 

SCAL 


























csearlh 
c 
c 
c 
c 
c 
c 
c 
c 


SLiSROUTlNe SEARCIi(SPAR,t.p/!H,LNPMUO. Ef-IAGJ 

CALIU bEAKCM KUUI INE 

SLOKOLIINE SEAKCHC SPAR, r{ PAH, LAFMCU.E FLAG* 
“this* FHCCHAH has PREPARED BY 


JCHN P. CLAASSEN 
glen E. ELUIQTT 

UNIVEHSITY OF KANSAS CENTER FOR RESEARCH 


C 

C 

c 

c 

c 

c 


This SUBROUTINE EXECUTES ENTRY INTC THE CONTENTS OF THE 
FILE ANC SEEKS CALIBRATICN INFORMATION, BOTH NORMAL AM) 

baseline, the program anticipates that butsi types cf 


00OOOA3B 

COOOOA31 

oonoo'.AO 

CCCOOAAl 
00000',', 2 
COOOOAA3 
00000'.A<j 
CCOOOAA5 
OOOOOA'tA 
CCOOOA47 
00000448 
00000449 
C00004SO 
00000451 
00000452 


CALIBRATICNS ARE PRESENT IN THE INITIAL RECORDS. H0REVERCC0004S3 



C 


WFE'I This NOT THE CASE, THE FILE IS SEARCHED UNTIL BOTH 

C0000454 


c 


TYPES ARE FOUND. BY-PASSED RECURDS ARE STORED ON UNIT 8 

00000455 


c 


SLBSEClEHT PROCESSING. 

C0000456 


c 



C0000457 

C000Q7 



DIMEilSIUN 0ATA(5C). 1DA1AI50), 

00000458 



* 

SPARI85), RPARIJ5), SCALI4), RCALC21 

00000459 

C00007' 



CCMMDN /INPLT/ HCFl, DATA 

00000460 

CC000 7 



CUUIVALfcNCE (CAT A (1 J , ICATAl 1 1 ) 

COOOC461 

CC00C7 



LOGICAL EFUAG, ECFl, ECF8 

00000462 

tCOOC7 



LGuluAL NORMAL, BASE 

C0CQ0463 


c 



C0000464 


c 


READ CATA RECORD FRCH TAPE UNIT 1 

00000465 


c 



C0000466 

(C00C7 



IPASS=P 

00000467 

000010 



nRITE (&,6G) 

CG000468 

C00013 

60 


FORMAT! 12F RECCRD NO.*, 13X, 7HHES SAC6, 19X , Ih* ,20X, 

00000469 



* 

10KPARAMETERS///1 

OO0OQ47O 

000013 



CALL READl 

C0000471 

C00014 



IF(EOFl) GC t'c 130 

C000C472 


c 



00000473 


c 


CHECK IF RECORD IS CALIB DATA 

CC000474 


c 



C0000475 

CC0020 

10 


IF (1DATAI41 .EC. 0 .CR. ICATA(4).EG-2J GO TO 50 

00000476 


c 



C0000477 


c 


PREPARE CISC FOR BY-PASSEC RECORCS 

00000478 


c 




C0000479 

C00027 



RErfir^U 6 

CO00O4B0 

C00031 



IPASS=1 

CCCC0401 


c 



000U04B2 


c 


ERROR MESSAGE 

CC0CO4S3 


. -C._ . 



00000464 

C00033 



hRltE(6,2CI 1CATA(41 

CCCC04B5 

C00040 

20 

F0RMATI5X,1H*,5X,27HFILS STRUCTURE INCCMPAT 1 BLE , 7X, IH*, 5X , 

00000486 



« 

33HFIRST RECORD NOT CALIERATICN CATA,5X, 6HM0DE= ,I3/J 

C000C467 

0000^0 



HRITE(6,120) data . 

CC0CQ4B8 


C 



0COQ0489 

_ , 

c _ 


SEARCH FOR CALI BRAT ICN RECORD 

CGC00490 


c 


- 

00000491 

C00046 

aO 


hRlTEIB) CATA , 

C0000442 

CQ0053 

35 


CALL READl 

000004S3 

C00054 



IFIEOFU GO TC 100 

0000045 A 


c 



00000495 

_ 

c 



CHECK IF CALIB DATA 

00000456 


c 



0000049 7 

000060 

<tO 


IF nCAlAI4J .NE. 0 .AND. 1 DATAC4 J.NE.2) GC TO 30 

CC0Q049B . 

C00067 

50 


•IC0UNT=0 

0C0O0499 

CC0070 



u.=o 

00000500 


c 



0000050 1 


c ^ 

- - - 

HLTMDRAn. CALIBRATION INFORHATICN 

CCC00502/ 


c 



0000050^ 

CQOOfX 



NORMAL = .FALSE, 

CC0005(^ 

C0Q072 



BASE = .FALSE. 

ooooosffy 

C00073 



CALL CAL U (SCAL, RCAL,RPAR, NORMAL, BASE 1 

coooo-5^t& 

C0007I 



IP (.NUT. NORMAL .AND. .NOT. BASE! GO TO 35 

0000<5S1C7$ 

C001C4 

. 



, IF (.NOT, NCRMAL) GO TO 160 ... 

0000050^ 

000106 



IF (.NOT, BASE) GC TO IBO 

00000309 

C00107 



IF( IPASS.EL'.O) GC TO 70 

CCCC0510 

ooouo 



ENDFILE c 


C00112 



REWIND 0 

0000051 1 

000U4 



CALL CRUNCH(B,LAI A, lOATA.IL, ICCUNT , SC AL ,RCAL , 

00000512 




_ SFAR,FPAR,LNFM0C,EFLAG,60FB) 

C0000513 

000131 

70 


IFIEOFU GO TL £0 

00Q00314 

C001J5 



CALL CHUNLKU ,CATA. (CATA, IL, ICOUNT , SCAL iRCAL , 

C0000S15 



4 

SPAR,RPAR,LNPMU0,EFLAG,EUF1) 

C0000516 

000152 



IF(EFLAG) go to 190 



c 



00000517 

. _ 

_ c ... 


. .. NCRMAL TERMINATION MESSAGE 

CC00051 8 


c 



CC000t> 19 

C00156 

uO 


WRITE (6,90) IL, ICOUNT 

CO0C052O 



FIGURE A-4. (a) SOURCE LISTING FOR SUBROUTINE SEARCH 
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000166 

VO 

rUKMAr (U,6Xf Ih»,2X,36HCCHPlErEL) CCMPL lA 1 1U.96 FUR THIS flLEt 

C0O00521 


« 

IX, lM*,5XtlOHTllERi: wERt,l5, 

C0000522 



24h HARMNG FIAOS OeNERArECJ 

00000523 

C00166 


RETURN 

C0000525 


c 


00000526 


c 

ERROR KESSACe 

C000052/ 


c 


C0000528 

0001.67 

100 

V.RITE 16,1101 

00000529 

000173 

1 10 

RORHAT t HX,lh«,5X,27HFILE STRUCTURE I ACCHPATI bLE, 7X, Ih* , 5X, 

00000530 



25FUNACLE TO FIND CALIU DATAl 

00000531 

0001/3 


Vikl Tr«6,l201 DATA 

CCOO0532 

000201 

120 

FCkHAr(/lQ(5022/)/l 

00000533 


C 


C0000534 


c 

AbURT FILE AND SEER NEXT INSTRUCTION 

00000535 


c 


C0000536 

000201 


EFL A6= • TRL E » 

00000537 

COOZO^ 


RETURN 

C0000538 

C00205 

1 30 

WRITE16, 1401 

00000539 

000211 

140 

FORMAT lllH empty FI LEI 

C0000540 

000211 


in A6=.TRLE. 

00000541 

00021^ 


RETURN 

CC000542 


C 


C0000543 


c 

SLA8CH FOR HISSING CALI6 DATA 

C0000544 


c 


00000545 

000215 

150 

WRITE (8) data 

CC0DC546 

000222 

loO 

CALL RcADl 

C0000547 

C00223 


INEUFII GC TC 100 

C0000548 

000227 

163 

IF nCATA(4) .NE. 01 GO TO 150 

C0000549 

C00230 


GU TO 65 

C0OC055O 

C00231 

170 

nSI Tt (81 data 

00000551 

000236 

laO 

CALL REACT 

CCOOC552 

000237 


IF(F.uFl) lO TC 100 

COOOC553 

C002-'i3 

185 

IF nCATA(41 -AE- 2} GO TC 170 

00000554 

0002-iiS 


CO TO C3 

C0000555 

C002A6 

190 

i.RITE(6,l9Sl 


0002S2 

195 

F0RMAT(///llA,lM'>,5X,2dhexCESSIVE NO. OF BAD RECORDS, 6X, 1H« 1 


C0O252 


IflEOFll RbTOHA 


000257 

200 

•vEADl 11 


C00262 


IFtEUFjll 250,200 


CC0270 

250 

RETURN 


C0C271 


END 

CCOO0556 






FIGURE A-4. (b) SOURCE LISTING FOR SUBROUTINE SEARCH 
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(Confinued) 

Table A"1 DefIni^ion of Variable Used in SEARCH enfries in DATA 


RCAL 

EOF! 

EOF8 

EFLAG = 
IPASS 

NORMAL = 
BASE 

ICOUNT = 
IL 


RAD Cal ibraf ions (NORMAL, BASE) 

End of Fife for Unif 1 

End of File for Unif 8 

File Error Flag 

Flag for Dafa on Unif 8 

Flag for Normal RAD Cals 

Flag for Baseline RAD Cals 

RADSCAT Performance Error Counter 

Record Counfer 


B. Subroutine CRUNCH 

N - '■ ■■■■- ' 

1 . Theory and Design 

The reduction algorithms are applied in or from subroutine CRUNCH. The 
routine was designed to process RADSCAT data whether it consisted of radiometer 
measurements (RAD only mode) or both scatterometer and radiometer measurements 
(alternating angles, fixed angle, or short SCAT modes). In proce ssin g the records 
on an independent basis, there Is a tacit assumption that coupling between 
polarization is negligible. It was also necessary to equip the routine with the 
ability to read records from file code 01 (row data tape) or file code 08 (scratch 
unit). A description of the Internal operation of CRUNCH Is shown in the logic 
diagram of Figure A~5. 

In the opening steps of CRUNCH the data record is extracted 
from the file and verified by an in-line routine. The first sixteen entires In the 
data record are checked for type and magnitude. If an invalid entry is present, an 
error message identifying the bad entry is given and a bad record counter is 
incremented (NGR). If more than 10 bad records are encountered, an error flog (EFLAG) 
is set and control is returned to SEARCH. If not, the next record in read (statement 
100) and the process repeated. If the record is valid, it is further examined to 
determine if it is a calibration record. If It Is, CALIB is called (statement 220) to 
extract and average the calibrations . The next record is then read and the process 
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INITIALIZE 
CH\N. IND. AND 
OUTPUT ARRAY 



STATUS 

message 


WRITE OUTPUT 
PARAiV.HERS i 
2 


TRANSFER RAD 
AND OTHER 
PARAMETERS 




it 


transfer scat 

PARAMERERS 

I 


INVERT 
SCAT DATA 


FIGURE A~5. DESCRIPTIVE LOGIC DIAGRAM FOR 
SUBROUTINE CRUNCH 


Driginal Page is 

op POOR QUALHT 
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repeated. If if is nof calibraHon dafa, fhe mode enfry Is ^ransIated and fhe SCAT 
and RAD receiver polarization are checked for agreement when the instrument is not 
in the "RAD only" mode. If they disagree the record is declared invalid and treated 
as a bad record.* If the measurement occurred in the "RAD only" mode or the 
polarization agree, the record is then processed in the steps following and including 
statement number 119. 

Initially the record counter is incremented, the SCAT channel indicator 
initialized and the output array cleared. CRUNCH then calls THERMO to compute 
and check receiver temperatures. An Interpretation code is updated after which 
subroutine ANGLE is called. ANGLE computes the incident and cross-track angles 
and again the interprefatlonal code is updated. Subroutine RCONV (radiometer 
convert) is called to compute the antenna temperature. If no scatterometer data is 
present, control advances to statement 200 where the computed results are 
transferred to the output array. Otherwise when SCAT data is present, DOPCHK 
is called to compute the doppier shift and set the doppler filter band index. The 
interpretational code is then updated. 

In the following statements (up to number 170), the maximum on“scale SCAT 
.measurement is selected from the four output channels (DATA(8) through DATA(ll)). 
To select the appropriate channel output, the channels are tested for an output in 
excess of a value' that would cause the output to saturate on the next higher 
sensitivity channel . Note that channel 1 (DATA(8)) is least sensitive whereas 
channel 4 (DATA (11)) is most sensitive. The selected measurement is stored in 
DATA(8) after scaling it with some channel dependent parameters. If the SCAT 
return fell below or above the channels, an interpretational flag is set In statements 
between 170 and 175. 

Following the channel selection, the range gate setting is compared with the 
aircraft altitude. If the altitude is not within 300 feet of the range gate setting, an 
interpretation flag is set. 

Subroutine SCONV (SCAT convert) Is then called to compute the normalized 

input power. Once computed, the scatterometer measurement Is Inverted for antenna 

pattern effects with an in-line statement (number 190). The antenna inversion is 

based on the notion of an equivalent pencil beam replacing the actual main beam 
T- — 

The implication is that the instrument didn't operate correctly. 
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an^enna paftern. Computer simula^ions have shown ^ha^ fhe method is extremely 
accurate (errors ^0.1 dB) for narrow beam^antennas . Briefly the equivalent beam- 
width notiqn is based on the fact that the normalized input power 

2. —o 


ice.) = J/ P C9-) dA 


where 

= normalized scattering coefficient 
= radar range to elemental area dA 
= normalized antenna pattern 
= elemental area 

can be approximated by an algebraic expression through the following considerations, 
For narrow beam” antenna 


lce.)= JJ 


sins R d© 

(cos ©') 


where cos 9' Is a projection factor, may be approximated by 


ice.)- 


0-"(6o) 


R (cos ©c 


) J j p^(6, 4') sin 0 d© 


The Integral may be computed numerically from a description of the antenna pattern 
and equated to the effect produced by a pencil beam of width 9eq , i.e.. 




sin 0 d0 



%/s 

StU 0 d0 d<^ 


z 

4 
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Therefore t'he normalized scaH'ering coefficieni is given by 


(T COo) = 


A I C0o) R cos e. 


ir 0 




Bu!" no^e that R = k / cos 0„ so thai 


^ (e«) = 


4 I (Qq) k 

TT* 0g^ COS 0, 


where h is ^he aircraft' al^itude. The enfry 0^^ is compufed by t'he engineering 
rouHne WIDTH described in Appendix E as well as SecHon IV F of i-he fexi". 

The scafl'eromei'er and radiomej-er product's are ihen sfored in Ihe out-put' 
array and Ihe output' array Is transferred to the output tape (unit 02). 

The next data record is then read (statement 100) and the sequence of 
events repeated. If an end of file is encountered on the input tape, control is 
returned to SEARCH, 

2. Program Listing and Variables 

The source listing for CRUNCH is shown in Figures A"6(a) through A~6(d). 
The variablesemployed in CRUNCH are defined in Table A~2. Entires in SPAR 


TabI 

e A-2 

Definition of Variables Used in CRUNCH 

DATA 

= 

Input Raw Data Vector 

IDATA 


Integer Equivalent of DATA 

SPAR 

= 

SCAT Parameters 

RPAR 

; 

RAD Parameter 

SCAL 


SCAT Calibration Vector (Chonnel 1 to 4) 

RCAL 


RAD Calibration Vector (Normal and Base) 

VERI 

= 

Data Validation Parameters 

IVERI 

= 

Integer Equivalent of VERI 

EOFIFC 

= 

End of File Indicator for FC 01 

EFLAG 

= 

File Verification Indicator 

IL 

= 

Record Count 
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C00017 

C00017 

C00017 

OOOOL7 

C00017 

C00017 

C00017 

C00017 

000017 

C00017 

CQ0017 

C00017 

000017 


■ .{>tUT i i\r CKurtCMlI CtUAlA, lUATA, IL ,[CCJNT ,bCAL,KLAL, ’ 

* il AR.i'P/K.L.NFHUL', EPLAU, EOH FC » 

CCKONCM OKUHCH 

(. 

C ins FKCGBAH HAS PREFAKEU IlY 

C 

C JCHA P. CLAASSEH 

C CLtN E. ELLIOTT 

C 

C _ __ IIMVEFSIIY OF KANSAS CENTER FOR RESEARCH 

c 

C 

L THF PRIHABY FUNCTION OF THIS ROUTINE IS TO CONVERT SCATTER- 


C01JOOSS7 

(cnccssa 

ccooosav 

C0QC0560 

C0000‘j6l 

CC00CS6<! 

OGOOOS03 

C000056A 

COOOOboS 

CC0CC5U6 

COOOOSbT 

CCOOCS6B 

C0OS1OS6V 


C 

C 

C 

C 

w 

c 

c 

c 

c 

c 

V 

C 

c 

c 

c 


ORETtf ANO PACICFETER OUTPUT VCLTACES INTO THEIR RESPECTIVE 
NORMAL UGU INPUT POHERS. THE CCNVERTEO SCAT UATA IS 
SLliSLLUENILlt INVERTED FOR THE SCATTERlNu LOEFFICIEHT. 

PRIMARY COMPUTATIONS ARE PERFORKEC IN SUORQOTINE CRUNCH 
WHICH CALLS A NUMBER OF IMPORTANT ROUTINES. FROM CRUNCH 
INTERNAL CALIBRATION INFORMATION IS klTHURAWN FROM THE JATA 
RELORCS AND APPLIED TO SUBSEOUENT MEASUREMENTS. Th£ VIEW 
angle is CtPRECTED FOR AIRCRAFT ORIENTATION. A WARNING FLAG 
IS GIVEN WHEN THE RCLL ANGLE CAUSES EXCESSIVE DEPULAPI ZAT ION 
A FLAG IS ALSC generated WHEN THE DOPPLER SHIFT EXCEEDS 
TFL UCPPLEP FILTER BANOHlDTH. ThBRMISTER DATA IS 
checked FCR REASCNABLE TEMPERATURES AND TEMPERATURES OF 
CRITICAL CCMPCNENTS ARE HITriURAhN FOR SUBSfcJUENT LSE IN 
CONVERTING RACIOHETER MEASUREMENTS. 

UIHENSICN DATAISO). ICATAISO),' SPAR105J, RPAR175Jf SCALI'i). 

® RCALI2I, AkRAYl2Q)i lARRAYtZOl. ALT ( A) 

OiMENSICN ICMGI3) 

UIMENStUN VERmStZIi IVcRUIS,2I 

UIHENSICN VHAXllSJt VMINIIS), IVMAX(15)» IVMIM15) 

Ego 1 VALENCE ( VMAXUlf'IVKAXi 11 ) i ( VH I N ( 1 1 , I VM I N I U I 
EUOI VALENCE I VEPl (1,1) .VHlN(l) 1 , tVEPI (1,2) ,VMAXll) J 
EQUIVALENCE (VERI(1,1),1VERI(1,1)I 
EQUIVALENCE lAPRAYll) , lARRAYIlH 
LOGICAL ECFlFC 

DATA ItVM.IM I )■ I=l,6)/1, 0,0, 1,0,0/ 

...DATA IlVMAXt I1,I = 1,6)/3,A,2,6,1,1/ 

data (vm1MI).I=7,IO),(VKAX(I),I = 7,10)/A«0.0, A*10.5/ 
p^TA 1VERI(11,.1),1VERU11,2)/0,1/ 


CCCCC370 
00000571 
00000572 
00000573 
CC00057A 
C0000575 
00000576 
CCC0C577 
.00000576 
CC000579 
CC00C5B0 
00000581 
C0000582 
COOOC5E3 
0000058A 
00000585 
CC00C586 
C000C567 
COOOOS88 
C0OCO589 
CC000590 
CC00C551 
. 00000592 
00000593 
CC000594 
CC0CO5S5 
00000596 
C0CCC597 
00000598 


000017 



oATa (VMlNm ,1 = 12, 14), IVKAXU), 1 = 12, 141/3+O.Q, 10.5, 10.5, 

10.5/CC000599 

C00017 



DATA IVERKlS.l), 1 VERl (1 5,2 ) /O, 3/ 

CCOOOeCO 

C00017 



DATA ICriG/1,0,2/ 

CC000601 

000nl7 



CATA ALT/2000. Oi 5000. Or lOOOO.Ot 200CQ.0/ 

CG0CC602 


c 



CC0C0603 


C 


keau in next record 

00000604 


c 



00000605 

000017 



NGR=0 

CC000606 

C00020 

100 


1F(IREAC( IFO.NE.C) GO TC 230 

0000C6C7 

000026 

101 


1 = 12 

OGOflOoCS 


c 



00000609 


c 


VERIFY CONTENTS OF RECORD 

CC000610 


c 



CC0C0611 

000027 



IF( ICATAI 1 ).LT. IVERl ll-l, D.OR. 

00000012 



* 

lCATA(n.GT.lVERl(I-l,2)) GO TO 115 

0000PC13 

000041 



II 

C0000614 

0000 -«2 



IFUOATAIII.LT. IVEFI (I“1,11.0R. 

CC00C615 




ICATAt D.GT, lVERIU-1,21 ) GO TO 115 

CC000616 

CC0054 



00 111 1=2,7 

C0000617 

000055 

111 


IFt IDAtAI I) .LT. iVER III-l, ll.OR. 

COOOO&ld 




ICA7AI I) .GT. IVERHI-1,2) ) GO TO 115 

COOOC619 

C00071 



Oil 112 1 = £,11 

C0000620 

C00072 

112 


IFUiAIAl I) .LT.VER I( I-1,1J .OR. 

CC000621 




CATA(1).GT.VERUI-‘1,2)) Gu'TO 115 

C0000622 

C00106 



OU 113 

CC000623 

COO 107 

il3 


IFICATAI 1 ) .LI .VER I( 1-1, 1) .OR. 

C0000624 




CATMII.GT.VERK 1-1,21) GO TC 115 

C0000625 

000123 



GO TO iia 

GOOOC626 

000123 

lib 


HRirE(6,Il6) I 

00000627 

COOUl 

116 


FCkHAT(//llX,lF*,llH 8A0 REC ORO, 28X, IH®, SX,9HCOMPONENT , I 3 , 

cooao62ti 



* 

12H I S oARCAGE. // ) 


C00I31 



hRlTEI6.lL7) DATA 

COO 00029 

000143 

117 


FURHAT(/10(5022/)/) 

C0C0O630 

000143 

105 


NGR « NOR + 1 

C000063I 

0001 4 1» 



1FINgR.LE.10) GO TO 100 

00000632 

OOOLi>2 



EFLAG = .TRUE. 

OOOOOC33 

000153 



GU TU Z:>0 

00000634 


C 


. 

COOOOE 35 


c 


BRANCH If CAL IB DATA 

C0000636 

' 

C 



CC0006J/ 

000154 

lid 


IF (lDAfA14) .EC. 0 .OR. I0ATAI4).EJ.2) GO TO 220 

OOOOOC 38 


L 



COOUO640 


C 


TRANSLATE CER FA IN. ENTR Ifc S IN FILE 

C000064L 
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C00163 


000 lb3 
C00175 
000212 

000212 


000226 


000234 

000236 

000242 

000246 


000246 


000252 

C00262 


1000 


C0000042 
CC0C0643 
00000644 
00000645 
CC000646 
00000647 
0000064U 
C0000649 
C0000650 
CO0C0651 
00000652 
00000653 
G00006S4 
C0000655 

IMTI3LIZE CHANNEL INUICATOR FLAG ANO CLEAR OUTPUT ARRAY - 00000656 

00000657 


iOATA{ 3 1-1 UAl M3 1-1 

iJN“IUAIAl4) 

lOAT A14 I -I CHLl UK3 1) 

JLlATA(6l«l-ICATA(61 

CFECK FOR AGREEMENT IN RECEIVE PCL 

IF(10AIA(7).E0. I0ATA(12).0R.10ATA(3I.EC.0) GU TO 119 
WR1TE(6,1000) 1CATAI7). IEATA(121 

F0RMAT(//llRilH*t31H RECEIVE POLAR I CAT IONS DISAGHEEr BX, lh« , 5X i 
9HSCAT KfcC-i 13*5X»aHRAO KEC=iI3//l 
GU TO 105 


C00216 

1 19 

I CHAN =0 

C0000658 

000217 


00 120 1 = 1 »20 

C0000659 

000221 

120 

I ARRAtYI I J = 0 

CC000660 

000224 


IL = IL * 1 



CHECK TEHFtFATURES CF CRITICAL ELEMENTS 

CALL ThERNC tOAT A, I DATA ,RPAR, I L, LNP KC C, 1 FLAG, JWARH ) 

FORM ERROR CODE FESSAGE FOR SCAT AND RAD 

1ARRAYU2) = IhARM 
ICOUNT = ICOUNT + INARM 
IARRAV(12> = lARRAYl 121*10 ♦ IFLAG 
1C0UNT=ICCUNT ♦ IFLAG 

COMPUTE LCCK ANGLES. 

CA_L_LjANGLe tUATAiSPARi ILfLNPHODt IFLAG) 

AMEND ERRCR CODE 

1C0UNT=1CCUNTHFLAG/10 
1ARRAYI12) = IARRAYI121 *100+IFLAG 

“CCNVERT 'rACICMETER LATA 


00000661 

00000662 

00000663 

C0000664 

00000665 

COOOt)666 

C0000667 

C0000666 

00000669 

CC000670 

00000671 

C0000672 

00000673 

CC000674 

00000675 

CC000676 

C0000677 

0000067B 

00000679 

C00O06B0 

00000681 

CC0006B2 

C0000683 


000264 


0002 70 


000275 


000 301 
000304 


C0000684 
OC000665 
C0000686 
C0000687 
C00006BO 
C000Q689 
00000690 
C0000691 
00000692 
CCCCC593 
00000654 
00000655 
C0000696 
00000697 
00000698 
C0000699 
00000 700 
CC000701 

THIS SECTION SELECTS THE MAXIMUM CN-SCALE SCAT CHANNEL UUTPUT00000702 
ANU TRANShCRHS IT TO A NURHALI2ED RECEIVER INPUT PCWER, 00000703 
hHERC READING FALLS CUTSlOE OF DYNAMIC RANGE OF THE FCUR 00000704 
CHANNELS. IhE DATA IS YET SELECTED. HOWEVER. A FLAG CCOOC7C5 
ACLUKPANIES THE VALLE. TFE TRANSFCRKEO VALUE IS FOUND C0000706 


CALL RCONV (OATA.IUATA.RPARiRCAL. ID 
BRANCH IF RAO ONLY DATA 
IF (I0ATAI31 .EC. 0) GO TO 200 
CHECK DOPPLER 

CALL DCPCHK (C AT A , SP AR 1 1 L ■ 18 AND. LNFMCD. I FL AG ) 
,AME_WO INDICATOR 

IAkRAYI12J = IARRAYI12I *10*IFLAG 
1LUUNT= ICCUNTi IFLAG 

SELECT SCAT CHANNEL OUTPUT 



C 


IN CAIACBK 




00000707 


C 






00000708 


C 


IMll ALIZE FLAG 




CC000709 


C 






00000710 

CC0310 


IFLAO =0 




COODO 711 


C 






C00007 12 


c 


DETERMINE HHAT CHANNEL 

MEASUREMENT FELL IN 

CCQ0C713 


c 






C0000714 

C0O311 


I F 

[DATA <81 ,CC. SPAR(25) 1 

GO 

TO 

140 

COOD0715 

000314 


IF 

(DATA 191 -GE. SPAR(291) 

GO 

TO 

150 

00000716 

000317 


IF 

(DATA (10). GE. SPAH(29)) 

CO 

TO 

160 

CC000717 

000322 


.IT 

(DATA (iU.GE, SPAR(291) 

UO 

ru 

130 

C000071& 


c 






COOQQ719 


c 


DATA FELL BENEATH RANGE 




00000/20 


c 






CCOOC/21 

C00326 


DATA111I=0ATA(11M1.0 E-10 




00000722 

CC0330 


I FLAG =1 




CC000723 


c 






CU000724 


c 


APPLY SCALE FACTORS 




CC000725 


c 






COOUC726 


3 PA.GB IS 

Q-aALTEYl 
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CQOiH 

1 JO 

OATC. (uj= 6AT/. ( llM'6(>Am21J /SCAU4 ) 


CC0C0127 

C0033C 


ll,hAIJ = 4 


C0000728 

COOJJ t 


00 ru I/O 


CC000729 


c 



C0000730 


c 

UAlA 1 ELL AOOVE RANGE 


CC00C731 


c 



CC000732 

000337 

1^0 

IF IDATAIO) .LE- SPAR(27)1 IfLAO*!. 


OOQ0O733 


c 



C 0000734 


c 

APPLY scale factors 


00000735 


c 



00000/36 

000344 


OATA (6l» LATA ( b )»SPARau l/SCAL U1 


OC000737 

0003&0 

“ 

ICIIAN-1 


C00C0738 

000351 


GO TO 170 


CC000739 

000351 

1^0 

LAIA (Ul= LATAl SJ »SPAR1 191/SCAH 21 


00000740 

000356 


I CHAN =2 


C0000741 

000357 


U) TO 170 


00000742 

C00357 

160 

OATA (U)= CATA( 101*SPAR(20J/SCAL131 


CCPOC743 

C00304 


I CHAN =3 


00000744 


c 



CCCC0745 


c 

ENG SFLECTION of SCAT CHANNEL. 


00000)46 


c 



CC000747 


c 

AMINO ERROR INCICATCR 


C000074rf 


c 



CC0C0749 

000365 

170 

IARRAY( 121 =1 APBA^ ( 12I*10+IFLAG" 


C0000 750 

CCC370 


1 COON 1 = I crONT < ] F L AC 


C0000751 

0003/1 


1F( IFLAb.NE.l 1 GO TO 175 


CC000752 

C00J73 


If (HDU( 1 Li LNFPCDl.Ev.O) riRI T E<6, 1700 1 IL 


OCOCC753 

000410 

1700 

FURNATl U.,5X, 2F* 22HEXCEE0E0 DYNAMIC RANGE 

,16X,1H*1 

00000754 


c 



CC0CC755 


c 

CHECK RANGE GATE SETTING 


C0Q00756 


c 



C0OOO757 

000410 

I7b, 

IFLAG = 0 


CC000758 

CC0411 


IALT = tCATAdbJU 


C0000755 

000414 


CHFCKsALTUALI 1 


00000760 

C00416 


IKoAT^( 4-1) .LI. CHECK-t-300. .ANO. 0ATA(411 

.GT. CHECK-300.1 

CCOOC761 



GC TO 177 


00000762 

C00431 


I flag -= 1 



000432 


»BITEl6i 17601 check, CATA(41I 



000445 

i'/oO 

FGRHAriI6,5X,2F« ,32HNOT nlTHIN 300 FT. OF 

RANGE GATE, ex, IF*, 



« 

5a, llriRANGE GaTEajF 7,0i5X,gHALTlTuDB=,F7.0) 


C00445 

177 

lARRAYI12r= 1AFRAY1121410 + IFLAG. 


CCrOOC/63 

000450 


ICou'.T = I CLLNT+ IFLAC 


CC000764 


C 



CCC0C765 


C 

CCIPLETE CCNVEBSION of scat DATA 




C 



C0Q00767 

000454 

leo 

call SLUNV [OAIA, 10ATA,SPAR,13AN01 


0000C766 


c 



00000769 


c 

SCAT INVERSION. 


C0000770 


c 



COOOO/71 


c 



CCCCC772 


c 

CCHPARE SCAT PCLARIZATIQNS- 


CC000773 


c 



C0000774 

C00457 


IXH T=IUATAi6) 


CC000775 

C00464 


IBEC»I04TA17) 


CC00C776 

000466 


1F( IaMT.EO.IRECI GC TQ 190 


00000777 


c 



CC00C778 


c 

error message 


C000O779 


c 



000C0760 

000467 


IFIHUOl [ L.LNFRCCI .ED.OJ nRITE16, 1800} IL, 

IXrtrt IREC 

CC000781 

000511 

IdOO 

F0RMAT(Ie,5X,2F» ,37HDNA8LE TO HANDLE CROSS POLARIZED CATA, 

C0000782 



IX, 2H* ,4X,26HwILL ASSUME PULARIZEC DATA, 5X, 

C0000783 



lOHXMIT PUL.=, I2,5X, lOHREC. PCL.=, 

121 

C00Q0784 


c 



CC000785 


c 

INVERT SCAT CATA 


C0000786 


c 



CCOOC787 

OOOSll 

IVO 

LiATAId) * 4. 0»DATA( 8)ADATA141)4DATA(41)/(CCS(0. 017453253* 

ccono/bB 



CATAI 14) i *3. 141 59265*5 PAR U XHT+32) *SP AR nREC + 32) * 

CCCCO/89 



5.290304) 




c 



CC000790 


c 

transter scat INI period 


CC000791 


c 



COOCC792 

000533 


1ND*1DATA( 5) 


00000793 

C00535 


AKK4y(71»SPAK( INO + 31 


CC00C754 

000537 


IF( lUATAOl.FC.l) AKRAY(7)»SPAR(10) 


COOOO/95 


c 



OCOCC 756 


c 

TRANSFER SCAT LATA TO OUTPUT ARRAY 


00000797 


c 



OC00079S 

000544 


(ARRAVfSI “lOATA I6 1 


0G00C799 

C00546 


1aRRAY(6) =IUATA (7) 


COQOQdOO 

C0054 1 


ARRAY 18) • CATA (81 


CC00080 L 

000551 


ARKAV(9) « l0.O*ALL'G10(DATA(0)J 


cocuoao 2 

C00562 


GO TO 210 


ccrooaoa 


c 



00000804 


L 

TRAMSi CK RAC CATA TC CUTHUi ARRAY 


CCC0Q805 


c 



C0000806 

C00562 

t'OO 

1 ARRAY! IZl^lANRAY (12 1 *1000 


CC00Q8C7 


IS 
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C00565 

210 


6KRAY (101= DAT/ (13) 

CCCC0SC8 


C 



CC000609 


C 


r TRANSIER CTI-EB SELECTED DATA 

CO0CO310 


C 



OOOOOU 11 

_C00567 



ARRAY(1)=LATA(1J . .. 

CC000012 

C00570 



1ARRAY(2J=I0A TA(3> 

OOOOOS13 

C0C572 



ARRAY (3)=LATA1 M ) 

COOOOO 14 

000573 



AKftAY(A)= UATA115) 

COOOOttlS 

0005/5 



AKKAY(ll)=CATAt40) 

CCCC0616 

C0057t> 



ARRAY {la)= CATA141) 

oooooaiT 

CC0600 



ARRAY (141= DATA (471 

COOOOtUU 

C00601 



ARRAY (15)= DATA (461 

cooooaio 


C 



ccocoazo 


c 


WRITE TC OUTPUT TAPE 

00000821 


c 



CC000d22 

000603 



CALL HRITc2( ARRAY) 

00000823 


c 



C0000824 


c 


IATERHECIATE PRIAT-CUT 

OOOOOS25 


c 



C0000fi26 

C00604 



IF (HODIIL, LKPRCDJ .EU.O) WRI TE(6f 2000) IL, ARRAV(IO), ARRAytd), 

00000827 



* 

IChAN, lEANC 

CCCCOaZa 

000653 

aooo 


FORMATUcf 5X,2H«<‘f4X,29HCOHPLET60 BECCRO CCHPUTAT I ONS i4X 1 2 H**t 

C0000629 




4Xtl4F AATENMA T£HP= , F6. 1 . 2X t IbhSCATT ER I NG C0EF=,Ei2.6, 

C 000 0830 



* 

ZX, 8HDHANNEL=,12t2X,9HD0P BAND®, 13) 

COOO0631 

C00653 



GO TO 100 

CtC0Cb32 

00Q654 

2 20 


CALL CAL)d(SCAL t RCA L t RPA P iNORpiAtf BASE ) 

00000633 

CC0650 



IF(.NUT.ECFIFC) GC TO 100 

CC000634 

000656 

230 


WRI TE(6i2C00 ) 1L« ARRAYilOJt ABRAV(0)t ICHANt IBAND 

00000635 

. C00702 



.RETURN 

CC0CC636 

C00703 



EWO 

00000637 
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Table A“2 DefinUion of Variables Used in CRUNCH (ConHnued) 


ICOUNT 

= 

Number of Interpretationa! Flags 

LNPMOD 

— 

Line Print Modulus 

IFC 

= 

Unit Code Number ,(1 or 8) 

NCR 

= 

Number of Bad Records 

ALT 

= 

Range Gate Settings 

CHECK 


One of Gate Settings 

IXMT 

= 

Transmit Polarization Indicator 

TREC 

r: 

Receive Polarization Indicator 

IND 

= 

SCAT Integration Selector 

ARRAY 


Output Array 

lARRAY 

~ 

Integer Equivalent 


and RRAR are defined In Section III D; whereas, the entries in ARRAY are defined 
in Section III C. 

C. Subroutine CALIB 

1 . Theoiy and Design 

Subroutine CALIB extracts, averages and displays calibration parameters from 
calibration records. CALIB was designed to anticipate calibration records in a format 
generated by the RADSCAT instrument. Although additional flexibility is provided 
so that baseline records, which are entered by hand,may be inserted in almost any 

fashion. The baseline records may occur anywhere in the file except within a 

* 

group of four norma! calibration records. In regards to sequences of normal 
calibrations, the subroutine will orient itself on a set of four by interrogating 
whether the first normal calibration record causes SCAT outputs to saturate in 
channels 2, 3 and 4. This feature was embedded in the routine when it became 
apparent that the RADSCAT instrument would occasionally drop calibrations on 
the first few SCAT channels. A descriptive logic diagram for subroutine CALIB is 
shown in Figure A-7. 

The baseline (BCAL) and normal calibration counters are initialized upon 
entry into the routine. Also the SCAT channel (CALI (I), I = 1,4} and RAD 
accumulators (CAL 1(5), RCL) are cleared. If the type of cal Ibration record (normal 
or baseline, see Section III B) is normal (IDATA(4) = 0), processing Is directed to 
statements Including and following 4. Otherwise, if It is a baseline record, it is 
accumulated in RCL, the baseline counter is Incremented, and the next record is 
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read. In sfatement 4 and the one following, the type of record is again determined. 

If additional baselines records are present (IDATA(4) = 2), the above processing is 
repeated; if it Is normal calibration data, the record is examined to determine whether 
the calibration indeed occurred on channel 1 (see DO loop terminating In statement 
5). If not, the next record is read and the type is again determined. Once the 
program is aligned on the calibration for channel one, it extracts it In statement 7 
and pulls the SCAT calibrations from the subsequent records for channels 2, 3 and 4 
and the normal RAD calibration from the last SCAT calibration record (Channel 4). 

The normal cal counter Is incremented and the calibration accumulator update. 

The next record is read and the type determined again. If it is a calibration 
record the above steps are repeated. If not (measurement record), the unblocking 
routine READI is called to backspace a record (statement 10). The steps 11 through 
60 determine what kinds of calibrations were present and how many. The accumulated 
calibrations are averaged and scaled accordingly. The types of cal Ibrations present 
are reflected by setting the logical variables NORMAL and BASE to TRUE. The 
calibration parameters are then printed and control is returned to the calling routine, 
either SEARCH or CRUNCH. 

2. Program Listing and Variables 

The source listing for CALIB is shown in Figure A~8. The definitions of the 
variables are listed in Table A~3, 


Table A”3 Variables Used in CALIB 


DATA 

— 

Input Raw Data Vector 

I DATA 

— 

Integer Equivalent of DATA 

SCAL 

- 

SCAT Calibration Vector 

RCAL 

= 

RAD Calibration Vector 

RPAR 

- 

RAD Parameters 

NORMAL 

-- 

Normal Cal Indicator 

BASE 


Baseline Cal Indicator 

BCAL 

= 

Baseline Cal Record Counter 

CAL 


Normal Cal Record Counter 

RCAL 

= 

Baseline Cal Accumulator 

CALI 


Normal Cal Accumulator 
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COOOlO 

COOOlO 

CCOOlO 

COOOlO 

COOOlO 


OOOOlO 

COOOll 

C0Q012 


iUilKOOl IHE CAI m(SCAL,l<CALiRPAI',.\li)KMALiriASl » 

CCAUli CALtUHAnoN SUROUTIME 

C 

C This THCORAH HAS PREPAHEU tJY 

C 

C JChA p. CLAASSEN 

C GLbN E. ELLIOTT 

C 

C ■ UNIVERSITY OF KANSAS CENTER FCR RESEARCH 

C 

C TIUS SUUROUTINE EXTRACTS CALIURATICN IKFORMATION FROP 

C THE RAOSCAT LATA FILE AND DISPLAYS IT. 

C 

CCHMOA /INPUT/ EOFl, LATA 
LOCICAL NCRHALi UASE, EOFl 

OIMENSION DATAlPOl, ICATMII, SCAUll, RCAUII, RPAKdl 

OIMfHSlCN CALUSl, CAL2IS1 

ELUIVALEACE (CATAUlflCATAIlIl 
C 

C IMT1ALI2E CAL COUNTERS 

C 

bcal=o.o 

KCL = 0.0 
CAL=0.0 


COOOUdJU 

CCOOOBJV 

COOOOOAO 

COOOOUAl 

OOOOOOA.! 

CCOOOUAS 

OOOOOB4A 

C00008AS 

00000046 

C0000U47 

COO0OB4O 

CC00C84O 

OOOOOtfbO 

C0CCC651 

00000652 

LC000853 

C0000S54 

C0CCCt55 

COOOCdSC 

CC000657 

cooooosa 

CCC0C85S 

00000860 

C0000861 


C CC000862 

C CLEAR CAL ACCUMULATORS 00000663 

C CC000864 

C0D013 00 1 1=1,5 ■ " C0000865 

C00014 1 CALl(l'l = C,0 " C000C866 

C C0000867 

C BRANCH HFEN NCT BASELINE CC0C0868 

C CC000869 

C00017 2 IFf ICATAtAl .NE,2I GO TO 4 O00C0870 

. _ . C _ ' 00000871 

C EXTRACT BAS'ElTnE ' ’ CC00C872 

C ■ COOOCB73 

C00021 3 RCL=RCL-*LATA( 13 1 ‘ COOOC674 

C ' 0C00087S 

C UPDATE ECal count C0CGC876 

C 00000877 

C00023 8CAL=aCAL<1.0 COOOOB78 


C 


00000879 


C00025 

C00026 

COU033 


CC0035 


000036 

C00037 

C00044 


C00044 

000045 


C00053 


C00055 

CG0056 

000057 

C00064 

COOObS 

C00071 

C00073 


000075 
COOO/6 
COO 102 


C READ NEXT RECORD 

C 

CALL REAOl 

IFlFOFll G(j TC 11 ' / 

4 IF(ICATA(4).EL.2I GO TO 3 

c ■■■■ ■ ' 

C EXIT hHEN NO MORE CAL DATA 

C 

IF(IDATA(41.NE.C1 GO TC 10 
C 

C.„ .. CHECK IF RECORD CONTAINS CAL CN 

C FIRST SCAT CHANNEL 

C < 

DO 5 1=2,4 

5 IFIOATAI I* 7) .LT.S.OI GO TO 6 

GO TO 7 

C „ . ... __ 

C FIND FIRST SCAT CAL 

C 

6 CALL READ! 

IFIEUFll 11,5 

C 

C EXTRACT SCAT 1 

C 

7 CAL2U1 = CATA( 61 
C 

C EXTRACT remainder OF CALS 

C 

.UO 8 1=2,4 

CALL RfcACl 
IFIEUFll GOTO 11 
15 1F( IDATA141 .NE.Ol GO TO 4 

8 CAL2I I1=LATA( 1471 
LAL2I5l=CATAa3) 

.. CAL=C.AL4U0. 

C 

C ACCUMULATE REPETITIVE CALS 

C 

DU 9 1 = 1,5 

9 CALHn=CALl( I)«CAL2(U 
CO TCI 6 

C 

C REI'OSIllCN DATA RECORD 


OQOCO880 
00000881 
CCQ00882 
00000883 
00000884 
CCOOOoES 
C0000886 
0CC00B87 
G0000883 
00000 839 
C00CC890 
C0000891 
C0000892 
fQ000893 
0000CH54 
00000895 
CCC00896 
C000C897 
COOOObSe 
CC00CO99 
0O0CO9OO 
CCQC09Q1 
00000902 
C0CGO9O3 
C0000904 
CCC0C9G5 
00000906 
CC000907 
00000903 
OOQOQ909 
00000910 
CD0C0911 
00O0O912 
C0000913 
00000914 
CC0CQ915 
00000916 
C0000917 
CC00091B 
00000919 
C 00009 20 
C0000921 
COOOQ ,22 


OBIGINAli PAGJB IS 
OJ’ POOR QtJALnY 
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000102 

C0O103 

coom 

00011 1 


C0011& 

C00120 

CCO120 


000122 


000126 


000132 

C00167 


000167 
0001 fO 


c 

10 


c 

c 

c 


11 



c 

c 

c 

13 
C 

c 

c 

c 

c 

c 

14 
100 


60 


CALL OCKSPl 


CHECK FCR PRESENCE OF 6ASBLINE CAL 

IFtOCAL.LT.l.O) CC TO 12 
BASE = .TRUE. 


AVERAGE AND ’TRANSLATE BASELINE 

KCALI1)= (RCL/tCAL 1«RPARU9)/RPAR(201 

CHECK FCR PRESENCE OF NORMAL CAL 

IMCAL.LT. 1.0 ) GO TO 14 
NORMAL = .TRUE. 

DO 13 1=1.4 

AVLRACE SCAT CALS 

scALi n=CALU n/CAL 

AVERAGE AND TRANSLATE RAO CAL 

RCAL C2) = tCALl(51/CAL )«RPAR( 19)/RPAR(20J 

DISPLAY CAL DATA 


♦ 

*■ 


WRITE (6,100) ( I.SCALII ),l = l,4l ,RCAL(21 , RCAL(l), CAL, BCAL 
FORMAT eilX,lH*10X,16FCALIBRAT10N C AT A , 13X , IH*» 3X, 

4 (2X,4HChAN , 1 2 , 1H= , F 7.3 > , 

/51X, lH»,5X,5HRCAL=F7.3,2X,6HReAS£=i F7.3, 

7F NCAL=,F6-C, 9H NdCAL=,F6.0) 

RcTLRN 

END 


C00C042i 
CCOOOV24 
CC000926 
00000926 
(C0CC927 
00000928 
CC 00 09 29 
C0000930 
C0000931 
(10000932 
CC0C0933 
CC000934 
CC000935 

CCC00936 

C0OC0937 

CC00Q938 

OOOOC939 

CC000940 

CC000941 

C0CCC942 

C0000943 

CC000944 

CC0C10945 

CC000946 

C0000947 

COOOO940 

C0000949 

■CC0CC950 

00000931 

C0000952 

CC000933 

COOC0954 

00000955 

CC000956 

00000959 




FIGURE A-8. (b) (conHnued) 
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D, Subroutine ANGLE 


1 , In troduction 

This routine computes the incident angle and cross track angle both of 
which are required to define the beam position on the surface with respect to the 
aircraft flight vector. When roll or pitch becomes excessive, as the case may be, 
so as to make the polarization difficult to interpret a flag is generated. Irrespective 
of the flagged condition the data is reduced in the standard way. The theory by 
which these angles are computed and the method by which the angles are applied to 
determine the polarization break-up at the surface are described below. 

2. Theory and Design 

g. Incident and Cross-Track Angles - Suppose the aircraft is vectored along 
the positive x axis of an unprimed coordinate system where the z axis corresponds to 
the local vertical. The orientation of the aircraft is represented in a primed 
coordinate system where axis is located along the fore-aft axis of the aircraft. 

They are related to the unprimed system by a drift angle about the z axis, a 
pitch angle 9^ with respect to the x, y plane and a roll angle 4>. about the x'axis 
(Figure A-9).* 

The antenna points at an angle 0^ with respect to the -z' axis. 0^ Is 
assumed to be corrected for the relative angle between aircraft and antenna platforms. 
The true incident angle 6 on the ground is desired. The incident angle may be 
derived from relation 

♦ 

COS 0 = - X,. • <j, (1) 

where is a unit vector in the boresight direction. Also of Interest is the cross 

track angle 4> measured in the x,y plane with respect to the “X axis as shown In 
Figure A- 10. Clearly we have 


t*c 


oun. 




( 2 ) 


It should be noted that these definitions of roil, pitch and drift coincide with those 
provided by the Litton Navigator LTN-51 . 
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1 



Figure A“9. Alrcraf}- Fllgl 1 ^ Geomeiry 


z 



y 

Figure A-10. Geometry Defining the Cross 
Track Angle. 


Now in the primed coordinate system the boresight dxis is described by 


- sm 0 ^ - cos (3) 
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The primed and unprimed coordina^e systems are related by the following successive 
transformation 



where 


A 


B 


ABC 


A 


y 



Substitution of (3) and (4) into (I) yields 

COS 0 5= sln0p si.u0(x. ®p 


(4) 


(5) 


which reduces to the expected result when 4>_ =0, namely, cos (0^ + 9^). 

When the drift angle is considered zero, a similar approach will yield 


icin C = 0) 


sIti cos & 0 L 


SL-nS^ cos 0p- sln0p cog cos 


( 6 ) 
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When <(> , =0, is clear from Figure A"!! fhaf fhe correcfion for <l\ is simply 
d d 

given by 


f J - C^j=o)] 


(7) 


or 


= - <#.J + 

t»n‘' [ ^ 


( 8 ) 


SlTl ©c 


si'll 9 ^ C6S0p ~ Sinfip COsSp, COS 


1 



\ 

\ 

y 

Figure A-11 . Cross-Track Angle Correcfion for Drift, 

• 

Equations 5 and 8 thus define the beam position on the sea with respect to 
the aircraft trajectory. 

b. Polarization Decomposition - When the aircraft pitches and rolls, the 
Incident polarization can decompose Into vertically and horizontally polarized 
components. When the unwanted component becomes excessive, the measurement 
becomes difficult to interpret; and, consequently, the reduced data should be flagged.* 
To determine the size of the undesired component, the degree of "de-polarization" 
is computed from considerations similar to those above. 

Correction for this situation is possible when polarized -measurements are considered 
jointly. 
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Wifhouf loss of generalify we may assume ^l^a^ ^he aircraff is vect-ored along 

fhe -X axis wifh zero drift-. If is pifched af an angle of 6 wifh respect to the -x 

P 

axis and a positive roll of<}>,. is induced about the -x' axis. The x,y plane forms 
the local horizontal and the z axis is pointed at nadir. See Figure A“12. 

The antenna is boresighted in the x‘ , z' plane with an angle 9^ from the z' 
axis. The vertical polarization emitted by the antenna is described by 

- cos - sin (9) 

along the boresight axis. The horizontal surface polarization at the boresight point 
Is described as 


= ' - sin "t cos 


c 


( 10 ) 


where 4>c corresponds to the cross“track angle for the zero drift condition. Now 


the percent power de-polarization is approximately given by 100 1“^ 
points Illuminated by the main beam. 


for 



Figure A“12. ''Depolarization" Geometry. 



Now I'he f■ransfo^na^ion be^ween fhe coordinaf'es is given by 



or 



/ cos e p 

0 

sl-n 9p \ 


j sin 0p 

coS<|>^ 

I 

^sin<kt. CC3 0 

r j 

K 


sItv 

cos4*j. cos 0p/ 





(12) 


Wi|-h fhe help of (12) we may fhus wrife 

I I = I sin ( CCS &p + stTi 8g, cos sin 0p ) 

+ cos slu 0^ sin I 

Now a suitable criterion for excessive polarization may be formed by requiring 

that 


a 

I ‘ I < ^ 


(14) 


where 0< ^ < < 1 , However, when the aircraft roll becomes excessive so that 
> 1 2 

‘ ^ 0.5 the interpretation of the incident polarization should 

be reversed from that defined by the instrument. Under this circumstance the 
polarization criterion should be written as 


2 
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where 0< ^ ^ 1 . This ^a^^er cril-erion should be useful when RADSC^T 

observafions are conducfed from an aircraff in a roll (blank ) at'Hfude» 

From SecHon IV.B.T we nofed fhaf for zero driff 

^ _ sin COS 

^ sin ©o, cos 9p i- sin 0p cos 9a 

The drifr angle does nof Influence fhe polarizafion decomposifion as a liffle 
fhoughf will substanf iafe . This is forfunare since driff angle is nof well defined 
in a roll maneuver. As a consequence fhe criterion is now formed on the basis of 
measured and computed parameters. 

Relationships (13)^ (14), (15) and (16) establish the depolarization parameters 
and the polarization reversal criterion. 

3. Program Listing and Variables 

The FORTRAN listing for ANGLE is shown in Figure A”13. The comment 
statements are sufficient to establish the logic of the program. The variables defined 
in Table A- 4 should be he Ipfu I . 


Table A- 

•4 

Definition of Variables in Subroutine ANGLE 

DATA 


Input Raw Data Vector 

SPAR 

= 

SCAT Parameters 

IL 

= 

Record Count 

LNPMOD 

= 

Line Print Modulus 

IFLAG 

- 

Excessive Depolarization and Reversal Flag 

DEGRAD 

- 

Degress to Radians Conversion Factor 

ANTA 

= 

Antenna Angle 

PITCH 

= 

Aircraft Pitch 

DRIFT 

- 

Aircraft Drift 

ROLL 

= 

Aircraft Roll 

THETA 

= 

Incident Angle 

ARG1,ARG2 

= 

Intermediate Storage 

PHI 


Cross-track Angle 

DEPOL 


Depolarization Factor 
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COOOlO 

CCOOlO 

000010 

COOOlO 

C00013 

C00015 


CANGLE 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c , 


SUOktlUllME angle (UAIA,SI>AK,lL.LKrt'CL,lELAU OOOOIHA 

CUMPUTE INLlbEAJCE ANGLE CCOOll^i 

bLOROLil INE ANGLE lU AT A>$ PAK « 1 L . LAPPOUt 1 FLAG) 000011V6 

OP001197 

TKIS FACGRAH WAS PREPAftEU BY OOOOUSB 

00001199 

JOHN P. CLAASSEN COOOIJOO 

GLEN E. ELLIOTT OOOOlJOl 


0000 1^02 

UNlVERSnjf OF KANSAS. CENTER, FEB RESEARCH OOOOU03 

0000 IJOA 

THIS ROUTINE COMPUTES THE INCIDENT ANCLE THETA FROM AIRCRAFT OOOOUOlJ 
AND ANTENNA PARAMETERS. CORRECTION IS MADE FUR THE RELATIVE OQO01JO6 
angle OETREEN the aircraft and ANTENNA PLATFURH. 000012C7 

THE bores IGhT CROSS-TRACK ANGLE PHI IS ALSO COMPUTEU. OOOOlJOO 

WHEN the RCLL CR pitch ANCLE BECOMES EXCESSIVE, AS THE CASE OC001J09 
MAY BE, SO AS TO MAKE THE RETURN CR EMISSION PQLARUATION OOOOIJIO 


DIFFICULT TC INTERPRET, A FLAG IS GENERATED. 
DtMENSICN OATA(l), SPARU) 

LtCREE-RAOIAN CCNVERSION FACTOR. 

DATA OEbRAD/0 .017 9532925/ ' 

INITIALIZE flag 

1FLAG=0 


00001211 

00001212 

0G0C121J 

00001219 

00001215 

00001216 

00001217 

00001218 

00001219 

00001220 

00001221 

00001222 


AVERAGE SCAT AND RAC ANGLES 00001223 

00001229 

DATA(19) = ICATA( 19 1+DATAI15H/2.0 0000122 5 

00001226 

, ..^convert to degrees , ... 00001227 

0000122B 

0ATA(19 J=CATAU91»SPAR(30) 00001229 

00001230 

COMPENSATE FOR ANTENNA PITCH RELATIVE TO AIRCRAFT 00001231 

. . 00001232 

CATAIi9J»LATAa9)+SPAPI311 . 00001233 

00001239 

C.ONVFRT. T® RACIANS 00001235 


OjIGINAL PAGB JH 
0|POOR QUALIT5J 


C00017 

C00021 

CC0022 

C00029 

C 

C 

c 

C00025 

C00093 


C 

C 

C 

C00096 

c 

c 

c 

C000 5 2 
(00057 
C00075 
C00109 
CQ01C7 
0001 12 

000115 10‘ 

000116 

C00121 

000122 30 

C00125 90 

C. 

C 

c 

C00130 


C 

c ■ 

C 

000162 
000163 
000172 
COOl/2 50 
C 
C 


anta= gataii9)*cegrad 

P irCH=OATA( 93 )*EECRAD 

drift=oatai991»cecrao 
roll =DATA{95)*CEGRA0 


COMPUTE INCIDENT ANGLE '.'/.LIT " ’ 

theta = STNIAMA1*SIN(PITCH1*-C0S lPITCHl + COSlANTAl*CaS (ROLLl 
theta = ACCSITFETAl 

■ -QtNVEFT"TO DEGREES. I ^ 

DATAtl9)=ThETA/CEGRAD ' ' ' 


COMPUTE CROSS TRACK ANCLE ", “ " ' 

ARGl=Sl.NIR0'LL')iC0Sl ANTA) “ 

AkG2 = SInUNTA)*CCSlPITCH)-SlNtPITCH)*caSlANTA)*C0S(RQLL) 

IF (ABSIARC2) .LT. 1.E-07J GO TO 30 
IF<ABSIARC1).LT.1,0E-Q7) GO TO 10 
PHI*ATAN2( ARG1,ARG2) 

GO 10 90 

■pHI=0. ■ - “ ’ • — 

I H ARC2 .LT .0.0) PHI=3. 19159265 
GO TO 90 

PllI = S I GN 1 1.57079633 ,ARG11 
DATA! 15I=PHl/OEGRAO-DATA(94) 

COMPUTE 'fFTETlEPOLARlZATICN ' 

OEPOL« 'sIn'iphI M 

. . „ ICOSI ANTAJ-SCOSIP ITCH I*-SI Nl ANTA » *CQS ( ROLL) *S I NI PITCH) )t 
LCSfPHI )*SIN(ANTA|*SINIR0LU 

' CHECiTFCR" POLARIZATION REVERSALS'. 

OEPOL = 0EP0l*0EPQL 
IFIOEPUL.GT. 0,500) GO TC 50 
GO TU 60 

OLPIIL ■= 1.0 - OEPOL 

I LAO REVERSAL 


}-IGURE A“1j. (a) FORTRAN LISTING. FOR SUBROUTINE ANGLE 
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00001236 

00001237 

00001235 

00001239 

00001290 

00001291 

00001292 

00001243 

00001299 

00001295 

00001296 

CC001247 

00001298 

00001299 

OOOC1250 

00001251 

00001252 

C0001253 

00001259 

00001255 

00CC1256 

00001257 

00001258 

00001259 

CCCC1260 

00001261 

00001262 

00001263 

00001269 

00001255 

00001266 

00001267 

00001268 

00001269 

00001270 

00001271 

00001272 

00001273 

C000I279 

0QQ01275 

00001276 

00001277 

C0001278 



00001^^79 

OOOOldoO 

OOOOliJbl 

00CC12fi2 

00001283 

0000128 ^ 

00001285 

00001286 

OOOOlZcY 

00001286 

C00012S9 

0000129Q 

00001291 

C0001292 

C0001293 

000012?^ 




FIGURE A- 13. (b) FORTRAN LISTING FOR SUBROUTINE ANGLE 




c 

C00179 lFtHUU( I LtLSPMCOI .EU.OJ hRITEl6,500l IL 

000210 500 f-URMATC16,5X,2H# 21HPULAR UA TION RE VE RJAL t ITX , IH* » 

000210 !FLAO=l 

C , . , ... 

c ' flag excessive oePOL. 

c 

C0021L oO 1HOEPOL.GT.0.O2J GO TO 70 

C00215 kETUkN 

000215 /O IFL AG = lFLflGt-10 

000216 iFlHLiuUL.LAPRCCJ.tO.OJ WRIT£C6> 700) IL, OEPOL 

000235 700 F0RHAT(It,bX,2F'*= ZAhEXCESSIVE DEPOL AR 1 2AT ICN , lAX , IH* , 

* 5X ,6H0EPCL=,F7.3) 

000235 UATAl AOJ^ULPOL 

000237 kETUkiN 

COO2A0 fcNO 
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E, Subrou^ine THERMO 


Temperal'ures within the RADSCAT microwave assembly are monitored at 
various points. Certain of these temperatures are essential to the reduction cf the 
radiometer data. All of the temperatures should be checked to determine whether 
the instrument is operating within its temperature limits. 

The temperatures are monitored by means of thermistors. Above 283*^ the 
thermistor output voltages are linearly proportional to temperature. This routine 
converts the recorded thermistor voltages, DATA{17) through DATA(33), to degrees 
Kelvin by using linear conversion constants, RPAR(23) through RPAR(56). Once 
converted the computed temperatures are examined to see that they appear in the 
operating temperature range specified by RPAR(21) and RPAR (22). Those- that fail 
outside this range are flagged with a message in the Conversion Report. A flag is 
also generated in the S/R validation code appearing in the output files. Routinely 
the temperatures critical to the conversion of radiometer data are reported by this 
routine. 

The above steps are reflected in the FORTRAN listing for THERMO shown 
in Figure A~14, The variables employed in the routine are defined in Table A”5. 

Table A~5 Definition of Variables Used In THERMO 


DATA 


Raw Data Vector 

IDATA 

= 

Integer Equivalent of DATA 

RPAR 


’ 'Radiometer Parameters 

IL 

i . 

- \ 

Record Number 

LNPMOD 

1 

Line Print Modulus 

IFLAG 

= 

Abnormal Receiver Temperature Flag 

IWARM 

= 

Warm Load Temperature Estimate Flag 


F. Subroutine RCONV 
1 . Introduction 

This routine converts the radiometer output voltage into an antenna (input) 
temperature. Radiometer baseline calibration RCAL(I) and normal calibrations 
RCAL(2) are employed to properly scale the output voltage DATA(I3) to a receiver 
Input temperature. Thermal noise contributions from front-end elements are subtracted 
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C00012 


OOOOl^ 

COOOLE 


C00013 

ccoru^ 

tCOOl.6 

C00020 

CC002b 

C00027 

000031 

C00033 


bUUKUUTl.NE IhERf'O t DAT A, I OAI A, R»> AK i 1 L . LtJPf'U D , I FLADf 1 WARf' I 
C 

CT HEKKU CHECKS reUPEPATbRE 

C 

C * " Thfs PROGRAM hAS PREPARED tlY 

C 

C JCHA P. CLAASSEN 

C GLEN E. ELLICTT 

C 

C UAUERSITY OF KANSAS CENTER fCR RESEARCH 

C 

C THIS ROUTINE EXAMINES ATM VOLTAGES AND CONVERTS THEM TO 

C TEMPERATURES. THE TEMPERATURE VALUES APE CHECKED TO SEE 

C IF THEY OCCUR IN THE PERMISSIBLE OPERATING RANGE. «HEN 

C THEY FALL CLTSIOE THE RANGE A FLAG IS CENERATEB ANC AN ERROR 

C MESSAGE DISCLOSING THE BAD VALUE ANJ KCNITORING LCCATICN 

c 'is printed- the voltage tc temperature parameters are 

C ENTERED FROM RPAR. THE STARRED STATEMENT IS TO uE 

C USED uPEN TEE HOT LCAD TEMP IS AVAILABLE. 


COnOlDGB 
0000106V 
CC001070 
00001071 
00001072 
00001073 
0000107A 
OQ00107S 
00001076 
0000107 7' 
0000107(5 
CCD01C79 
0000 10«0 
ooooioei 
OCC010E2 
OOOOlOtB 
OCOOlOflA 
OO001O86 
GOGO 1CE6 


C 

DIMENSICN CATAdlt RP AR ( 1 > , lOATAd) 

C IMTIALUE FLAGS ' ■ 

C „ _ 

IFLAG=0 

IhARH = 0 ‘ _ . . 

C __ _ 

C CONVERT TC TEMPERATURES ~ 

C 

J = 22 

DO 20 1= 1,17 
K = 1 +J 

UATA( 1+16J=RPAR(K1*DATA<I+16J+RPAR(K+1)+273.18 

T = OATA( I«i6 J . 

J = J5-1 

' (FI I.EV.17) GC TO 20 

IF d.LT.RPSRiai J .AND. T.GT. RP AR 1 22 I ) G'O TO 20 
C 

C FLAG AND PRINT MESSAGE WPEN NCT IN LIMITS. 


00001087 
CCOClCaB 
0000 1089 
D0001090 
00001091 
00Q01OS2 
00001093 
00001099 
C0001095 
00001096 
OOC01C97 
00001098 
00001099 
00001100 
OOQOllOl 
G0001102 

00001103 

00001109 

CCC011C5 


C 


COOC93 

10 

IFL 

C0Q099 


IF( 

C00066 

100 

FCR 






* 

COOOCfc 

20 

CUN 


AG=1 

MODdL.LNPMCC) .Eg.O) WRITE ( 6 , iOO J IL 

PAT(l6,SXf91F9 rECVR temp NOT MITHIN OPERATING RANGE A, 
5X, 12HLCCATICN NO. 13, 5X, __ 

13HTEMPERATURE =,F5.1,7H KELVIN! 

T I N UC 


C 

c 

c 


CHANGE AIR 


TEMP^O KELVIN 


CCC070 


DATA19U) = DATA196) + 273.18 


C 

C CHECKS FCR ESTIMATED WARM LOAD TEMPERATURE 

C 


C00072 
C00079 
CC0079 
OOOll 3 


IF (IDATAI36! .EC. 0) GO TO 30 ' =” 

IwARM = 1 

IF IHOD( I L, LNfPCC 1 .EG. 0) HR1TE16,1501 IL 

1‘jO FORMfiT CU,5 a,2H# .31HWARM LOAD TEMPERATURE EST I MATEO ,7X , , 1H9 ) 
C 

C LIST critical TEMPERATURES 

c ■ “ 


000113 

cooip; 

C00167 

C00170 


iO IFIMODdLiLNPHCDI.Eu.O) WRITE (6,200! IL, CATAI271, DATA1331, 

* CATA(22!, DATA(28!, DATA(96!, 0ATAI25! 

20U FORMAT tI6,SX,th» 21HCR1TICAL TEMPERATURES , 1 7X, 

* 1H*,3X,5HWARM=,F6.L,2X,9HH0T=,F6. 1,2X,9H0HT=,F6.1, 

* 2X ,3HSW=, F6. 1,2X,9HAIR=, F6. 1,2X,9HGDE=, F6. II 
RETURN 

END 


00001106 

000011C7 

00001108 

O0CO11G9 
OOOOlllO 
0000 dll’ 
00001112 
0COO1113 
00001119 
OOOOillS 
00001116 
00001117 
00001118 
00001119 
00001120 

C0001121 
00001122 
00001123 
00001129 
OOOOU25 
OOQQU26 
OOOOU27 
0000 1 128 
00001129 
00001130 
00001131 
00001132 


FIGURE A-I4. FCRjRAN LISTING FOR SUBROUTINE THERMO 

85 




ou^ to refer the receiver input temperature to the antenna terminals with the help of 
conversion cons tan ts 'RPAR and temperatures of input elements,. DATA (17) to DATA (33). 

2. Theory and Design 

* 

■ The radiometer conversion model was developed in TM 186“3. In view of the 
involvement In the derivation, only an outline of the derivation will be presented 
here . 

During baseline calibration an output proportional to a zero receiver intpu 
temperature is generated 


^ [ V, + “Tjr- c -T„ -T„ )] 0) 

where t = calibration integration period 

c 

RC= output Integrator time constant 
, Q= receiver gain established by AGC circuit 

V = output reference bias 

T.j^= hot load temperature referred to receiver input 
^2r~ temperature referred to receiver input 

During a normal calibration sequence the output integration produces a voltage 
given by 



( 2 ) 


This result yields the unknown gain factor 

Q = C V,- V, ) 

c 


( 3 ) 


where (1) has been employed. A measurement cycle yields an output voltage given 

- G(eT^,.-T,^-T^ r) 

■4(T,,-Ti:r) - ( 4 ) 


V = 

''m R c 




] 


Ciaassen, J. P., "The RADSCAT Radiometer Transfer Function and Its Application 
to the Reduction of RADSCAT Dafa," University of Kansas Center for Research, Inc . 
Lawrence, Kansas, December 1971. 
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where T 


ar 


is ^he I'eniperature a[-the antenna port of the Dicke switch. This 


result may be written as 

^ T,». ^ T 


gr 


a,r 


m 


C T,r - T,, ) 


2ir 


(5) 


where 




m 


Vm - Vb 
<■ - V. ) V-r, 


(6) 


and where (1) and (3) have been employed. From pages 8“11 of TM 186-3 it is 
noted that certain leakages at the Dicke switch require that certain terms in (5) be 
modified with reflection factors and so that the modeling equation actually 
becomes 




'ir +• >7.r- 


Our 


C"^tr s.<*} 


(7) 


Now the temperature at the antenna port of the Dicke switch can be described by 

(l-R.,)Tar +f 2 T 3 + f“,T, 

^t’sT> + t's Tt * fn fs"'''- 


( 8 ) 


where 


= desired antenna temperature 
Tj^ = leakage from opposite port 
Tg = Internal guide temperature 
Tj = temperature of polarization switch 
Tp = temperature of feed between switch and OMT 
. = temperature of OMT 
T^, = temperature of cutler feed 

= referral constants i= 1,2 ^..,,8 

It should be noted that the factor simply accounts for small but unknown contri- 
butions originating from front-end elements not well defined physically. Now let 




2 , 


(9) 


87 



and 


f m ' - T„ ) - e„ ( V, T, - Ts ) 


( 10 ) 


Then (7) may be wrlH'en as 


'f, Tu = ■^1', + c Ui T, - v; Te') 

-J>^T3 -y^Tp (H) 

- faTa -S't 


as 


02) 


or after division by J, 

Tn.+ 5, = ^, T, + <r^ V ( <r,T, - (^T^> 

-’.■J 

Now since Tp is no^ monitored If will be assumed fhaf Tp = -s~ Q.^ I ^ 

where a + a = 1 . In fhis case fhe conversion model simplifies fo the form 

<r f 


= o;T, t cr,r, - cr, X - - cr-Tj 

-''kTs -<Tjr„- <i|Ta - 


(13) 


The leakage term Tj^ is actually the antenna temperature ot the opposite polarization 
This term can be estimated from measurements at that polarization » There is no 
guarantee that it is always available; however, over the sea it can be estimated 
from T^ with reasonable accuracy when the sky Is clear. Hoi linger* has shown that 
the polarized temperatures are related by 


T^v = 


I ~ c 
I -t* c 


Tak 


where 


C = 0.00012 6? 

9 = observation angle in degrees 


Holiinger, J. P., “Passive Microwave Measurements of the Sea Surface," Journal 
of Geophysical Research , vol . 75, no. 27, pp. 5209-5213, September 1970. 
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The lahter es^ima^ing fechnique is employed by fhe subroufine since femperai'ures a^ 
f-he opposit-e polarization are not always available. Fortunately ? is small 
(estimated at .004) so that under some circumstances the leakage term can be ignored 
(by setting =0). 

The above model is incorporated in RCONV. The listing for the subroutine 
is shown In Figure A~15. The variables are defined in Table A~6 . 

Table A“6 Definition of Variables In Subroutine RCONV 

I = Receive Polarization Index 

J = Frequency Indicator 


RPAR(5) = 

, cr, 

RPAR(6) = ' 

Vi, 

RPAR(7) = 


RPAR(8) = 


RPAR(I+9) = 


RPAR(I+n) = 

^7 

RPAR(I+I3) = 


RPAR(I+15) = 


RPAR(17) = 


EST 


DATA (28) = 

's 

DATA (32) = 

’’o 

DATA (46) = 

Tc 

DATA (1 8) = 


G. Subroutine 

DOPCHK 


1 . Introduction 

This subroutine computes the doppier frequency shift induced by the relative 
motion between aircraft and sea. The doppier shift is compared with the doppier filter 
characteristic to determine whether the shift is within the bandpass. When it is, a 
frequency index IBAND, which is employed to locate the filter gain, is computed. When 
the doppier shift lies outside the defined doppier filter characteristic, the band Index 
is set at the appropriate extreme and a flog is set. 
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COOOlO 


CCOOiO 

COOOll 


C.KQJNV 

c 

c 

c 

c 

u 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0 

c 

& 


SUJKCUT Hb RCLNV IDAf/l, JCATA^KPAKt ftCAL , IL I 

CCAh/hRT KAUlGPfcTEK 

SOURlUIM: ftCCNV (UAT A. 1 UAT A, RPAK , PC AL . 1 U 


C 

C 

c 

c 


This PPCOKAM hAS PKEPAKfiU BY 

JCHA p. CLAASSEN 
GIEM E. ELUCIT 

UMVcPSHY OF KAASAS CENTER FCR KESEARCH 


THIS PKJCFAP CCNVEBTS RADIOMETER CUTPUT VOLTAGES TO AMENNA 
TEMPERATURE . 

UlHENblLN pATAi 1) , 1 1 AT A 1 1 J RPAR U > » RCALllj 

(-CRM POLARIZATION ANu f-REOuENCV INDEX. 

i = ii;ata( IZ J 
J = 1 LATAt 2) 

EXTRACT LCAU TEMPERATURES 


C000I2 


U = CATa(S3I 

CnOQI4 


T 2= DATA 1 2 2 1 

CC0015 

. 

1FIJ.E«.3) T2=C'ATA(27) 


c 

r" 

HCRH PEASUREMENT PARAMETER. 

000021 

V 

DATAll3t=-CLATA( U J-PCALU) I/IRCALI 


V» 

/' 

c- 

r 

APPLY lEPPERATURE FACTCHS. 

C00C2O 

V 

iJATA(t3) =LATa( 13 )*( RPAR I31+TI-RPaR( 

* KPAR( I+S J«DATA( ZM)-KPAR( 1*1 

* CATAI46J - BPAR 1 15)*CATA (10 


APPPCXIMATE KISSING CATA 


CCCCCUtC 
CCOOOVoI 
CCOCOSL2 
CC000VU3 
CCCCC'jfcA 
C0000D65 
CC000S66 
ooooovgt 
CCCCCSGd 
CCOOO'jt.'J 
OCOCCS70 
CC0C097 1 
COOO09 72 
COQCQ<> 73 
COOOOS74 
CC0CCS.T5 
CC000976 
CO000S77 
CU00097B 
0CCCC979 
OOOUO 90 O 
CCQOQVBl 
ccoooyei 
CQOQO<;t3 
00000'JB4 
ccQaQSdL- 
COOOOV&T 
CC0CCS68 
coooo 9 a 9 

Cf.0C099a 

C 00 fi 0 y 9 l 

C0000992 

OL000993 

CCCCC99A 

0CQ00995 

OCOOC99t 

CCOCCI997 

CC<l0399a 

00000999 

CCCCICCO 

OOOOIOfU 

C0C01002 



C 

APPkUX IMaT If^o parameter ' 

oncoLooA 


c 


OCCOIOO', 

C00062 


C=O.OOOi^*UATA 114 )#[}ATA 114) 

ccrciocb 


L 


OOOOlOOtj 


C 

FOR tiCfUACMAL PCLAPUATIJN 

CC0010C7 


c 


CCOOiOOU 

CC0065 


£ST=DATA(l3)*a.0tC) /(1.0-C) 

01001009 


c 


OOOOIOLO 


c 

FOR vertical PCLARUATION 

0000 10 11 


c 


00001012 

000071 


IF (1 ,Eu.l> EST = DATA ll3 ) # 1 1- 0-C ) / (1. 0*C J 

C0001013 


c 

hOKM ANTENNA TEH PERiJ T UK ES 

0000101*1 


c 


CCOOIO lb 

COOlOO 


UATAC13) = 0ATA113J ’ E S T *R PAR 1 1 +3 ) 

COOOlOlU 

COO 104 


RETURN 

OOOOlOlT 

C00105 


END 

0000 JOIB 


DEIGINAi; PAGE) ® 
lOF POOR QUALTO 


FIGURE A“I6'. FC.RTRAl'J LISTING FOR SUBFOUTINE RCONV 



2. Theory and Design 

For non-relaHvisfic speeds hhe fwo way doppier frequency sh iff Is known fo 
be given by 

r -- - f 

+ _, = 2 < IT. •«., ) — (1) 


where = boresight unit vector 

V = aircraft velocity 
^ = operating frequency 
C = speed of I Ight 

Denote as the cross-track angle and 0 as the Incident angle (see subroutine 
ANGLE). Then the boresight vector is given by 



= “ sin & cos 4 stn & 

3C?1 - COS & 

(2) 

Also 

■ 



So 



( 3 ) 


Zr ~ OOS 4*^ 


( 4 ) 

and 

ft "f V . p, 

t J - 2 ^ StTl 0 cos 


( 5 ) 


Expression (5) is the desired result. 

Now the computer routine determines whether the doppier shift falls in the 
doppier bandwidth as defined by the upper and lower limits, SPAR(23) and SPAR(22), 
respectively. If f^ is between SPAR(22) and SPAR(23), f^ is divided by the doppier 
filter gain sample interval SPAR(24) and the result rounded off to determine the band 
index (IBAND) from the centeV frequency. When the shift is outside the bandpass, 
the band index is set at the appropriate extreme and a flag is set. An error message 
is also appended to the Conversion Report. A listing of DOPCHK is shown in 
Figure A~16. The variables employed by the routine are defined in Table A"7. 
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cooou 

COOOll 

cooou 

cooou 

00001.3 

CC0015 

cooou 


iODRUUTINe DOKHK ( DATA iSPAK t 1 Lt t BAND t LNPMCU , IFLAGl 
COUPCHK CHECK DOPPLER SHIFT 

C SUBROUTINE DOPCHK (0 AT Af SPAR, IL , I H AND, UNPHOO, If LAG 1 

C ^ ^ ^ 

C This PRCGR.AH WAS PREPAREU 8Y 

C 

C JOFN P. CLAASSEN 

C GLEN E. ELLIOTT 

C 

C UMVEPSITY OF KANSAS CENTER FOR RESEARCH 

C 

C 

C THIS SUDRCITINE COMPUTES THE DOPPLER SHIFT TNOUCEC BY THE 

C REUAIIVE MOTION BETWEEN AIRCRAFT AND SEA. THE DOPPLER SHIFT 

C IS CCMf-ARED WITH THE DOPPLER FILTER CHARACTERISTIC TO 

C OETERMI-NE WHETHER SHIFT IS WITHIN THE BANDPASS AND WHEN 

C IT IS TO SELECT THE DOPPLER FILTER GAIN RELATIVE TO THE 

C CALlbRATlCN SIGNAL. 

C 

UIKENSIUN DATA! U , SPARUJ 
C 

DATA DEGRAD/0. 017^532525/ 

C 

C INITIALIZE FLAG 

C 

IFLAG=0 

C 

C CONVERT aircraft SPEED TC HETERS/SEC 

c 

VEL=DATA IA2J/3.0AB 
C 

C CONVERT TC RADIANS 

C 

THETA = DATA UA)«CEGRAC 
PHI=DATA(15UCEGRAD 
C ' . \ 

C 'CCHPUTE THE UGPFLER, SHIFT 

C . , ’ 

D0Pf-2.0«5PAR( I l<‘VEL/2SS793000.0*SIN(ThETAU'COStPHI ) 

C ■ ' 

C IF JCFPCER SHIFT EXCESSIVE FOkH FLAG 


OCOCU33 
OOOOU3A 
OOOOIU5 
OOOOU36 
OOOOU37 
00001138 
C000U39 
OOOOUAO 
OCOOUAl 
C000UA2 
00 00 U A3 
CCOOllAA 
OOOOUA5 
OOOOUA6 
OOOOIIA7 
COOOUA8 
OOOOUAV 
C0Q0L150 
COQOU51 
OOOOU52 
OOCpOU53 
C000115A 
0000U5S 
00001156 
00001157 
C00CU58 
OOOOU59 
00001160 
00001161 
0000U62 
00001163 
OOOOU6A 
00001165 
CG001166 
OOOQU67 
00001168 
00001169 
00001170 
00001171 
00001172 
C0Q0U73 


C 


OOOQU7A 


C00036 

CCOO'il 


OOOOAA 

COOOaT 

C00053 

C00(15A 
C00057 
C00060 
C00063 
C00O6 A 
OOOU5 


COO 115 
000 U6 


lF(DOP.LT.SPAR(22n GC TC 10 
I F ( DDP.GT.SPAR(23II GO TC 20 
C 

C FORM DCPPLER BAND INDEX 

C 

IBAND=-DGP/SPAR 12A) +-0.5 

IF tDOP.GT.0.0) iaAND= -COP/SPARt 2 A J-0 . 5 
RETUkN 

C fCft.'l FLAG AND ASSIGN BAND 

10 IbAND = -SPAR(22I/SPAR(2A) 

GU TO 30 

^0 IdANu = -SPAR (23) /SPARI2A) 

IFLAG=1 

IHIMOLlll.LNFHCCl.eD.O) WRITE (6,A0J IL',‘C ATA ( A2) , DATA (lA ), OOP 
AO FORMAT ( 16,5X,lh<‘rlX,i7h£XC£SSlV£ DCP F L ES , 2 1 X , 1H«‘ , 5X , 

« 13HAIPCFAFT SPD= , F 5 . 1 , 5X , 11 HV I 6 W ANG LE=, F 5 . 1 , 5X. 

* lAHDUPPLER SHIFT-, F6.0) _ , 

RETURN 

END 


00001175 

00001176 

CC00U7/ 

00001178 

00001179 

OOOOUdO 

OOOOUBl 

00001162 

00001183 

CI000116A 

OOOOUE5 

O0G01186 

00001187 

OOOOLlatj 

00001189 

00001190 

C000U91 

00001192 

00001193 


FIGURE A- 16. FORTRAN LISTING OF SUBROUTINE DOPC^IK 
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OBIGINAII PAGE XS 
OP POOB QUALITV 



Table A~7 Defini 

IFLAG 
VEL 
THETA 
PHI 
DOP 

IBAND = 

SPAR(22) = ' \ 

SPAR(23) = ' 

SPAR(24) = 


fion of Variables 
Flag for Excessive Doppler 
Aircrofl- Ground Velocify 
Incldenf Angle 
Cross-frack Angle 
Doppler Shiff 
Doppler. Band Index 
Lower Doppler FlHer Lim^^ 
Upper Doppler Filter Limit 
Doppler Filter Sample Interval 


H . Subroutine SCONV 
1 . Introduction 

SCONV converts scatterometer output measurements into a normalized input 
power. The scatterometer transfer function is employed to perform the conversion. 
The transfer function is developed below. Certain significant simplifications arise 
when SCAT normal calibrations are employed. 


2 . Theory and Design 

The block diagram for the RADSCAT scatterometer is shown in Figure A“17. 

The front~end elements have been decomposed into a calibration element G , 
transmission elements • HT'^ ^VHT reception elements 

^VR' ^HVR'^ ^HR ^VHR analysis and reduction of data. The transfer 

functions account for mis-match and loss. Calibration attenuators A^ through A^ 
are employed to calibrate the output channels 1 through 4, respectively. 

During calibration transmitter power is routed sequentially through A.j, A 2 , 

A^ and A^ and the respective channel output is noted. The resulting output voltage 
for the i channel is thus given by 



p 



’( - INTEGRATOR 
CONSTANT 


FIGURE A- 17. BLOCK DIAGRAM FOR SCATTEROMETER 




















where d = ^he effecfive dufy factor 

P^= transmitter peak power 

calibration integration period 

The factor d occurs in the above expression since only the center frequency of the 
power spectrum is sampled by the receiver. The calibrations thus serve as a measure 
of the transmitted power 


vf 

* ' A^, % 


( 2 ) 


'During measurements at VV polarization and incident angle 8^ the power 
arriving at the receiver Input is given by 


r, C9) r, (5 


VT 


ce)r„(3„„ ] (3) 


VT 


where 


I.. = Ji 




■4'Tr' 


R 


dA (4) 


a = normalized scattering coefficient for transmission at 
polarization n and reception at polarization m 

R = range to scattering element dA 

Pm = normalized antenna pattern 

T = antenna directivity at polarization m 

X = operating wavelength 

Similarly for observations at HH polarization the power at the receiver input is 
given by ^ 

~ (-47ri^ ^ 

**■ ^VHR ^ 

^ 3 
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vHr 


(5) 



The laH'er fhree ferms in expressions (3) and (5) may be neglected since products such 

as Ghvi^^YH' ^VHR^VHT'' ^HV^VHT' etc,, are extremely small . Therefore (3) 
and (5) simplify to: 




Ad O % 
<4ir) ^ 


Q P I iQ)P Q 

VR y vv ^ V ^ vr 


( 6 ) - 


and 



2 ^ 
A d 




( 7 ) 


The corresponding on-scale output voltage for these measurements occur on say, 
channels j and k, respectively. The outputs are given by 




‘^R“j 






'Tty 


and 


( 8 ) 


where ^m ~ fi^sasurement integration" time 

m — 1,2, ...,y 

A subscript is employed on to identify the different integration periods associated 
with six angle indicators (1 through 6 ) for the ALTERNATING and FIXED ANGLE 
MODES and the integration period ^7 unique to the SHORT SCAT mode. 

Nowequations( 8 ) and (9) may be employed to express the receiver input 
power. Thus we may write 




a" d' P, 




^VR^HT 


( 10 ) 


V 


U m 


a' d' P, 

(4 'rr i ^ 
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(n) 


iO bi- Q-i 1- 'v 



Now P may be replaced by expression (2) for ^!^e corresponding out’puf channel fo 
yield ^he desired resuI^s 







(4^ 


( 12 ) 


and 




_v, 

V 


k m 

tc) 


'm 


^ fe 

^HR ^HT 


i4'7T ) 

2.^ 


(13) 


If is no^ed fhaf fhe receiver gain ferms have dropped out. 

Subroutine SCONV implements the algebraic expressions for lyy and 
Actually the SCAT channel selection routine in CRUNCH applies the factors A^ and 
(m = j or k) since they are channel dependent. Missing from these algebraic 
expressions is the doppler filter gain factor which is a function of the doppler 
frequency shift. The gain is normalized with respect to the gain at zero doppler 
since the calibrations are made there. An equal ly spaced sample version of this 
filter characteristic must appear in the lower part of SPAR. The filter gain at zero 
doppler must appear at SPAR(50), Down doppler filter gains must be stored in entries 
between SPAR(51) and SPAR(85) whereaS' entries for the up-doppler gain must appear 
between SPAR(34) and SPAR(50). Presently the gain function is sampled at 250 Hz 
intervals. To compensate the measurernent the relative doppler gain must divide 
the measurement, DATA (8) . 

The FORTRAN coding for SCONV is shown in Figure A~18. The variables are 
defined in Table A“8. 


Table A“8 

DATA{8) 

INT 

I DATA (3) 

SPAR(3) 
SPAR(INT+3) 
SPAR{IBAND'!-50) = 
SPAR(17) 


Definition of Variables 

Selected Channel Output, V. 

jm 

Integration Period Index 

RADSCAT Mode 

Calibration Period 

Measurement Period 

Relallve Doppler Filter Gain 

Calibration Attenuation, G 

' c 
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iOtiRUUTlNC SCCSV (OATA, lUATA.SI'AR, IBftND) 

oooo ioi‘y 


cscunv 

SCAT COMVCRSION 

00001020 


c 

SUJRCCTlNfc SCONV <IJATAi lUATA, SPAft.I bANt)» 

00001021 


c 


C000);022 


c 

“" ‘TMS PHC6RAM HAS PR6PAKE0 BY 

COC.01023 


c 


0000102A 


c 

JCHA P. LLAASSEN 

OCOO 102B 


c 

GCtN e. ELLIOTT 

00001020 


c 


0Q00L027 


c 

LNIVERSnv OP KANSAS CENTER ECR RESEARCH 

ODC0102U 


c 


O0O0LO29 


c 


00001030 


c 

THIS SUBRCUTINE CONVERTS A SCATTEPO-HETER NORMALIAEO COTPUT 

00001031 


c 

VCLTACe TC THE SCATTERING INTEGRAL. THE RESULT IS LEFT 

CC0C1032 


c 

IN 0A1AI8). 

00001033 


c 


C000103A 

C00007 


UIHENSKjN UATA(l), ICATAU), SPAR(l) 

O00C1035 

CC0007 


CATA PI/3.LA1S927A 

CC0C1C36 


c 


00001037 


c 


0000103B 


c 


PCOQ1039 


c 

APPLY CCNVERSIUN PARAMETERS 

OOOOIOAO 


c 


OOOOIOAI 

COO 007 


LiATAlul= CAIA(U)t(A.*Pl ) < 2S9 793 COO . /SP Aft < 1) )*«=2 

000010A2 


c 


CCC01043 


c 

SET INTEGRATION TIME INDEX 

C'0001044 


c 


C000104S 

CC0014 


INT = luAIA (6) 

CI0001046 

C00015 


IF tlOATA (J| .60. i) INT =7 

00001047 


c 


0C00104d 


c 

COHPENSATE FOR INTEGRATION TIMES 

00001049 


c 


Q000105D 



LATAIu)= LATAlU 1=^5PAR(3)/SPAR( INr + 3 ) 

CTiOOiOSl 


c 


CC0010S2 


c 

SET PCLARliATION INDICES 

00001053 


c 


0000 1054 

C0002^ 


IXMIT =1UATA 16) 

C00ni055 

C00025 


IREC = ILATA 17) 

00001056 


c 


OCOO 1057 


c 

LCMPENSATE FCR DOPPLER FILTER C HARACTERI S TI C 

CC.C01054J 


c 


rcooio59 

C00027 


UAlA(tS)= CATAie )FSPAR ( I6AND + S0) 

CC001060 


c 

* 

CC001061 


c 

APPLY GAIN- EACTCkS 

00001062 


c 


00001063 

toooi 1 


DAT A(d)=CATA(d) ♦SPARl ir )/ ( SP AR ( I XH I T-i- 1 l')*S FAR ( I XMI T + 25 ) -*> 

C0001064 


* 

SPAR t IREC-fl tij «SPAR( I REC+25 1 ) 

OCOO 1065 

coon^ i 


kF,TURN 

Cf CC1Q06 

C000^2 


END 

00001067 




FIGURE A- 10. SOURCE LISTING FOR SUBROUTINE SCONV 
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Table A-8 DefiniHon of Variables (ConHnued) 


IXMIT 

IREC 

SPAR(IXMIT+11) = 
SPAR(IREC+15) = 
SPAR(IXMIT+25) = 
SPAR(IREC+25) = 


Transmit Polarization Index 
Receive Polarization Index 
Transmit Transfer Function 
Receive Transfer Function 
Transmit Antenna Gain 
Receive Antenna Gain 


V. SPECIAL ROUTINES 

A . Tape Read Routines 

FUNCTION subroutine IREAD and subroutine READl are employed to read 
the Input Raw Data records. IREAD permits CRUNCH to read records from either 
file code 01 or 08. IREAD colls READl to unblock records from file code 01 . The 
records on the scratch unit (08) are not blocked and so can be read from IREAD. A 
second entry point in READl permits subroutine CALIB to backspace one physical 
record (DATA). Read functions on unit 1 may be performed by simply calling READl . 
The source listings for IREAD and READl are shown in Figures A" 19 and A“20, 
respectively. 

B. Tape Write Routine ^ ' 

Subroutine WRITE2 performs WRITE functions on file code 02. Output records 
(lARRAY) are accumulated in IBLOCK until It is full at which time IBLOCK is written 
to file code 02. A second entry point in WRITE2 performs the END OF FILE function. 
However, the remainder of IBLOCK is first filled with zeros before putting an EOF 
mark on file code 02. See the listing for WRITE2 in Figure A~21 for further insight. 

C. Subroutine ABORT 

To purposely stop the execution of the program and abort the job, an abort 
function was provided. The circumstances underwhich an abort is executed are 
described in Section III B. If the system has an abort routine available that dumps 
more registers and memory than this fatal execution will, It is advisable to replace the 
one provided. A listing of subroutine ABORT Is shown in Figure A~22. 
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FUMOTION IREADttK) 


00000J79 


C 

C 

c 

c 

c 

u 

L 

C 

C 

c , 
c 
c 
c 

C0CPC3 

CQ00U3 

COOOOJ 

tCOQC3 

CC00C6 

C00007 

000013 

C00015 1 

CC0Q22 

C00026 2 

C00027 3 

GOG031 


THIS r-PCOPAH WAS HRlPAKEC bV 
OLPN E. ELUOTT 

ONUEBSllY OF KA^SAS CHNTEP PCR RESEARCHi INC. 

This RCtTlNE PERMITS SUbRCUTIRE CRLNCH TO REAO FROK FILE 
CCOe 1 OK II. RECCRUS ON FILE UNIT 1 ARE bLCCKED IN 
ORLOPS CF 10 PHYSICAL RECOPUS- EACH RECORD CONSIS'TINC. 

OF SO hOBOS- SLilROUTINE REAUl UNBLOCKS THE RECORCS- 
CN TFE CIFERHANC THOSE RECORDS CN FILE JNIT 0 ARE NOT 
OLOCKEC AND AS A LONSEUUENLE A SIMPLE BINARY READ 
IS EHFLCYEO. 

CC0G03U0 
00000301 
CC0C03b2 
00000363 
C00003fc<i 
CGOOQ3E5 
CCCCCiBfa 
00000367 
CCOOOifct 
00000309 
COO0C390 
GOQ00391 


LOGICAL ECFl 

COMMON /INPUT/ fcO FI , C AT A ( 5 0 1 
IREA0=0 

IFl IFC-Etf-t» GC TO I 
CALL REAUl t , / ‘ 

IFIEOFl) IREAC=1 , 

return'. ' ■ • 

REACIdl' LATA ■ ' 

IFtEUF,U> ' I 

IREAU-l'j 

RETURN 

cNU 




? 



r 

bOOKUUI l^r Rt/lL 1 


CCCfO J'/2 


c 

This I’RCGRAM WAS hREPARGU BY 




c 

(jLEN E. ELLIOTT 




c 

UMVPRSUY UF K.AN5AS CENTER FCR 

RESEAKCHt INC. 



c 

c 

This SLEAOUTINE REAUS BLOCKS CF 

10 PtlYSlCAL RECORDS ANO 



c 

SEPARATES THEM, RECORDS FILLED w 

ITH 99999. ARE 6Y-PASSEC 



c 

UNTIL AA ECF IS EACOONTEKtC. A 

SECONDARY ENTRY FCINT 



c 

IS PRLVIUED TO PERMIT BACKSPACIAG INUIVIUUAL RECORDS 



c 

/» 

WITHIN A BLOCK. 



CC0002 


LOHMUN/ INPUT/ lECFfCATAI 501 


COOOCJS3 

C00002 


LOGICAL lECF 


O0OCO394 

C00002 


OlMENSlcA BLOCK (50t 10) 


00000395 

cooooa 


DATA NDAtAUT/C. 1/ 


ccoQcasc 

C00002 


I EOh= .F ALSr . 


COOOC397 

C00003 


IFINIJX.NE .G) GL TC 10 


OOOQC390 

C00004 


NOX = l 


C0000399 

COOOCb 


feadinctj block 


CC000400 

COOOL2 


!FIEOr-,NCTI 100.10 


00000401 

coco L 5 

10 

IFIrfLCCKI 1 ,WJX) .EO. 99999.1 GO TO 90 


CC00C4C2 

C00022 


DO 20 I = 1,30 


C0000403 

C00023 

20 

CATAIU = ULCCtC (IfNOXJ 


cccnc404 

C00033 


NOXstiOXtl 


00000405 

C00034 


I FINOX.GT.IO) NCX=0 


coooo4oe 

000037 


RETURN 


C00004C7 

CC0040 

90 

RFAG(KUT) 


CCCCC40B 

COOQ43 


IFIEUP.NUT) ICC, 90 


00000409 

CCC046 

ion 

IEOF=.TRUE, 


CC00041C 

C00047 


NDX=0 


00000411 

conoso 


RETURN 


COOC0412 

C00061 


Entry bCkspi 


00000413 

COOOSo 


IFINDX.LT.ll AOX=ll 


CC00C414 

C000c2 


iMjX=NDa-1 


00000415 

OCOOC3 


KETUkN 


CC00041t 

C00064 


END 


00000417 





FIGURE A-20. LISTING FOR READl 
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SUbU((UTINE nKnE2 (IAUkAy> COOOOAIO 

c 

C This I'RCGKAH HAS PREPARED BY 

C GLEN E. ELLIOTT 

C UNWERSnV OF KANSAS CENTER' FCR RESEARCH, INC. 

C 

C This ROUTINE BLOCKS 20 HORO RECORDS INTO 20 RECORD 

C GROUPS teLCLKS) AND WRITES TEE ULOCK DOT TO FILE UMT 

C 02. when INSUFFICIENT RECORDS EXIST tU FILL A BLOCK 

L TEE REPAIMNL PORTION IS PiLLIC WITH 2ER0ES AS PRGYIDEO 

C AT THE SECLNCARV ENTRY POINT. AN EOF IS APPENCEC 



L 

TC ThE FILE 

TFEREAFTER. 



C00003 


UlllENSiLN IBLCCKI20, 

20), lARRAYCU 


00000410 

C00003 


DATA INULX /!/ 



C0000420 

(C00C3 


DC 20 I = 1,20 



CC0C0421 

C00005 

20 

I BLUCMI . INDEX) = lARRAYllI 


00000422 

OCCOI5 


INDEX = INDEX + 1 



' 00000423 

CC0016 


II- ilNUL-X -Lt. 20J 

RETURN 


00000424 

C00021 


rtKlTE (2) lisLGCK 



00000425 

C00026 


INDEX = 1 



GC000426 

000027 


RETURN 



CC0Q0427 

C00030 


ENTRY ENDF2 



C000C428 

C00CI36 


IF (INDEX .LE. IJ 

GO TO 30 


C0000429 

CCOOA2 


DO 25 I s INDEX, 20 



C0000430 

CCOOAJ 


DO 25 J- 1,20 



C0000431 

COOOAA 

25 

liSLDCKU.lJ = 0 



00000432 

COOOS^f 


«RITE I2J IBLCCK 



0000C433 

COOOLl 


INDEX=1 



00000434 

C0CD62 

30 

ENDFILE 2 



00000435 

OOOOE4 


RETURN 



C0000436 

C0006B 


END 



CC000437 



FIGURE A-21 . LISTING FOR WR1TE2 
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SLBKObTIfie *iUU[<T(A) 

C000036U 


C HALT PkiJbRAH IlM CASE OF ERROR 

00000369 

(OOOOJ 

UlMENSJOR 1ECRC2] 

00000370 

COOOOi 

DATA lEOR/lOHEAC CF fte6,lHL/ 

00000371 

C00003 

bkt TE (OtIODJ A 

CCCCC372 

cooou 

UlO FUEHAT ( ICXi lOfABCkT COOEf 1X,A10) 

00000373 

Cf.OOU 

CALL EHLIF2 

CCCCQ37<, 

C0001<' 

HRHr(2) IE0R,< I ,l = 3,9J 

CC00037S 

0000^4 

ti=A/U.O«^1.0 

CC0CC376 

CC0030 

STOP 

C00003/7 

CC0032 

E.-IO 

C00OG37e 



FIGURE A-22. LIST FOR SUBROUTINE ABORT 
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APPENDIX B 


OUTPUT PROGRAM 


L INTRODUCTION 


The oulpuf- program (OUTPUT) lisis groups of oulpul files produced by ihe 
conversion program. Exf-ernai direclives permil ihe user lo seiecf file lypes lo be 
lisfed. This appendix describes fhe operalion of OUTPUT, In ihe descripHon below 
if will be helpful fo refer fo fhe logic diagram of Figure B-1 and fhe lisfing of Figure 

B“2. 


II. DESIGN AND OPERATION 

Upon enfry Info OUTPUT a confrol card is read. The formal and confenf of 
fhe confrol card Is treated in Section III C. The card content is printed and its 
validity is checked. An invalid command terminates the program. The tape on unit 
2 is also rewound. The file label is read and compared with fhe validated directive. 
If fhe file label does not agree with the label, fhe file is bypassed and a bypass 
message is printed. The search continues until fhe appropriate file is found; if in 
fhe search the end of reel is encountered/^ message to that effect is printed and 
another control card Is read (statement 1). 

Specified file types are listed 80 records at a time. Each record is identified 
with a record number. The listing continues until fhe file is exhausted (EOF) at 
which time another file of the same type is sought. 

The variables employed in OUTPUT are defined in Table B-1. 


MISSD 

Table B-1 

Definition of Variables 

Directed Mission Number 

FLTD 


Directed Flight Number 

LIND 

= 

Directed Line Number 

RUND 

= 

Directed Run Number 
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Pa 





FIGURE B-1. DESCRIPTIVE LOGIC DIAGRAM 
FOR OUTPUT PROGRAM 
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COOOOJ 

C00CC3 

tOOOOi 

<00003 

<00003 

(0C0C3 

000003 

<00003 

C 00003 

C 00003 
COOOc 1 
COOO^l 
C0002<t 
C00026 






COUTOUT 
C 
C 

c 
c 
c 
c 
c 
c 


l■H 110 kAH CLTPUUINPUT.OOIPUT tfAHEa, TAP C 3 > INFO r.TAPeo-UUrPUr J 

OUTPUr PKCORAK 

This program has PREPAKEO by 


I 

100 


11 


JChh P. CLAASSEN 
GLEN E. ELLICTT 

LNIVERJIIY OF KANSAS CENTER FOR RESEARCH 

'iNTEOtR PITO, BLNDr FLT, RUM 
UlMENSION CATA<20) F DATE (3 J 
UlMENSION LSTYPl'tli IHOUE(A)f 1PL12) 

CIHEhblCA IE0P(2» 

LOGICAL tCFZ 

EgUlVALeNCE ( PCCE,DATA(2I 1, 

( tPl.CATAI 51 ) , ( IP2,DATA(6J ) 

CATA IbUR/lOhtnC OF KEEfIHL/ 

UATA iH3L)E/5MR,t.f'iHS.S, t‘tHF.A.,APA.A./ , IFL/1HH,1HV/ 
DATA LSTYP/6HFLTL IN , 6HM ISFLT ,6HSPEC IF t3HALL/ 

get directive and file id. 

REAL-(3tl00J LTYFE, HI SSOi FLT Di LI NO,RD ND 
FORllAl I Ab iAX,'i 151 
lMEUf,5J 10, U 
REWIND 2 

WkITEl6,200I LTYPE, MISSD, FLTD, LIND, RUND 


Cf 0012 VS 
00001296 
00001257 
0 G 00129 B 
00001299 
00001300 
O 0 C 013 O 1 
00001302 
00001303 
C 0001305 
0 OC 01305 
0 OOG 13 CC 
000013 C 7 

£CC 013 CU 

00001309 

C 0001310 

00001311 

00001312 

CC 00131 J 

0000131 A 

00001315 

00001316 

00001317 

00001318 

00001319 


COOC44 

200 


FORMAT (Ihlt///1CA13HCCN7P0L CARD- ,3X i A6,4X ,4 152/ 1 


0000 44 



IjD 20 i s 1,4 


CC0C46 

20 


IFI LTYPE. EC. LSTYFII n CO TO 2 


C00052 



HRI TE(6,7r0) 


C000 55 

700 


FDRMATIlhl ,/// 10X22HUNRECCGMZED CIRECTIVE) 


000055 



STOP 



C 


FIND THE FILE. 

00001320 

C00057 

2 


LONTINLiE 

00001321 

C0OO57 



wRITEIo.iOO) 


000063 

300 


hURHATUHU 


C00063 



kEAC(2) PISS, FLT, DATE, LIN, RUN, FREC, FEED 

00001222 

000104 



IFIMSS .NE . IECRI 1) ) CO TC 13 

00001323 

C00106 



IHFlT.NI. lEOR (2) ) GO TC 13 

00001324 

OOOILO 



Mr. 1 Tit ( 6 1 4G G ) 


COOllw 

4 00 


FORMAT {///10X10(1F*1,11FEND OF R~ EL , 10 ( IH* 1/// 1 


C00114 



GO TO 1 . 

„ 00001325 

C00113 

13 


CCNTINJS 

00001326 

C00115 



IF(LTYPE.E.i.LSTYK4n GC TO 4 

00001327 

000117 



IHMl SS.NE .31550 GO TO 3 .. 

O00C132B 

000l2>l 



IF< FLT.NE.FLTOI GC TO 3 

00001329 

000123 



IMUrYPE.Ew.LSl »Pt2) J GC TO 4 

00001330 

<00125 



IFiLIN.NE, UNCI GC TO 3 

00001321 

C00127 



IKLTYPE.NE.lSTYPO)) GC TO 4 

00001332 

000131 



IFCKUN .Ej.RJNDJ GC TO 4 

00001333 

C00133 



hRI TE( 6,500) KISS, FLT, LIN, RUN 


0001^6 

5 00 


FORMAT (//IOX13FEYPASSEO FILE, 4151 


coo 146 

3 


READ(2J 

00001334 

000151 



1F(E0F,2) 2,3 

00001335 


C 


WiUND IT, OUTPUT MEADINGS. 

0000 X 336 

000154 



CONTINUE 

00001337 

oooi:;^ 

V 


HRIIE( 6 , 604 ) 

00001638 

flOOlbO 

604 


FOkfiAT ( in , 5 JX, 14 MP ICRChAVE CATAI 

00001339 

OOOL^O 



LIN£N =0 . . 

00001340 


c 


BYPASS 1C RECORD. 

00001341 

000161 



READI 2 ) ICSET. ICl, 102 , 103 

00001342 

CQ0174 

5 


nk I TE( O.eco) MISS, FLT, DATE, LIM, RUN, FREQ, FEED, 

00001343 




ICSET, 101, 102, 103 

00001344 

cn 02^6 

bOO 


FGR( 1 AT(/ 9 P MISS ICN-, 13 , Sh FLIGHT-, 13 , 

0000134 b 




7 H OATE-, 213 , I 5 , 2 X, 11 H FLT LINE-, 1 X, 16 , 6 H RUN- 

1 13, 




13 H FRECUENCY-, F 5 . 1 , 7 H PEEL'-, 2 X,A 6 , 0 H IPAR-, 

I2.313/)CC001347 

C00226 



WR 1 TE( 6 , 605 ) 

OD001348 

C00232 

6 03 


FORMAT (/ 

00001349 




62H ktCJRJ TIME MODE INCID CROSS XHIT REC S 

INT 00001350 



4 

53h$CAT SCAT RAD DEPCL S/R ALTITUDE, 

OC 001351 


UH FLT 

62 H no. 

52ti 

IIM DIRECT 

6 2h 

63H 

1211 (DEG) . 


DB 

/, 


//} 


ANGLE ANGLE 
FACTCR 


POL 

COCE 


I DEG) 
IDEG K) 


(DEG) 


POL TIPE 


(SEC) 
IFEET) , 


000232 

000234 

C 00236 

000240 

000242 

000243 

C 00245 


DU e I-l,tO 
CALL RFAD 2 (UATA,E 0 F 2 ) 
IHi;aF 2 l GC IC 7 
12 L 1 (IFN»L 1 NCNT 1 

rtUUC^IMOCeiHODE*!! 
ii‘i»iPLnpui) 
1 P 2 »IPL( IF 2 UJ 


OUTPUT ONE DOUBLE COMPUTER PAGE. 


00001352 

00001353 

00001354 

COOOIJ 55 

00001356 

00001357 

0000135 U 

00001359 

00001360 

C 0001361 

00001362 

00001363 

00001364 

00001365 

00001366 


FIGURE B-2. SOURCE LISTING FOR PRCGRAAA OUTPUT 
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i t 14 ) 


C00^46 

C00<ia4 


6 

rtPl Tr-(0.60 

U LI 

6 01 

FURI-IAT 


C 



in 

* LINE 

TIME 


* NO. 


«■ 

* ' 



* (IX. 14, 2X 

. F6 . 1 

c 



in 

*SCAT 

SCAT 


* 

DB 

* 

* 



< Ell. 4. 2X, 

FO » 1 . 


^E^. <L/iT /U) t I = 


MATCH HbAUlNG 


MODE 

INCH) 

CPCSS 


AMULE 

ANGLE 


(OEGl 

(DEG) 

A6, 

F5.1.2X, 

F5.1 

RAC DfcPOL 

S/R 


f ACTUR 

CODE 

lOEG 

K) 


F 7.1 

, FB.3, 

IS, 


hITH FURHATS. 

XKIT KEC S IM 

FCL PUL TIME 

(SEC! 

6Xf A5f AI|4 Xt F5»3|2Xi 

ALTITUDE*. 

*. 

( FEETJ*. 

FIO.O . 



c 



f LT */, 


« 



DIREC */■> 


in 



(UEu( *) 





1) 


C 

c 



FOOT PAGE WITF CODE KEY. 

000264 



hR17E(6.6C2) 

CC0270 

602 


FORMAT (//14K KEY TC CODES-/, 





25X.32H 0 NO FLAG *, 5X, 



* 


,:1F iOOO POLAR REVERSAL/, 





'■ ' 23X,* 1 POSSIELV OUTSIDE RASGE GATE*,5X, 





22H lOOOO EXCESSIVE DEPOL/, 



in 


25X,32F 10 OLTStOE DYNAMIC RANGE ,5X, 



*Jc 


24H1CC0G0 REC TEMP ABNORMAL/, 



in 


25X,* ICO EXCESSIVE DOPPLER '»»4X, 





#1000000 HARM LOAD TEMP ESTIMATED#/, 





25X.37H CCHBINATICNS OF FLAGS CAN OCCUR./, IHIJ 

C00270 



GU 

TO 5 


c 



FOOT PARTIAL PAGE 

C00271 

7 


.jKI 

TL (6,l02) 

CO0275 



1F( 

LTYPS.EC.LSTYPIBJ ) CO TO 1 

000277 



GO 

TU 2 

C00300 

10 


STOP 

C00302 



END 



r CC0C1367 
n00Q136U 
00001369 
*,000013/0 
COOCIJ 11 
*,Q00Q13/<I 
00001373 
00001374 
Of 0C1373 
00001376 
000013/7 
000013 70 
CC0C1379 
OOOOX3UO 
00001361 
00001302 
CCCC1363 
00001384 
■ C00C13fc!) 

00001306 
OCOC13t7 
000(11380 
C0001304 
000013>10 
OOOOljSl 
0C0C1392 
00001393 
O00Q13S4 
00001395 
fC0C13S6 
00001397 
00001396 
00001399 
CC0C14r,C 
00001401 
CC001402 
aC0C140J 



igure B~2. (Continued) 



Table B-1 Definifion of Variables (ConMnued) 


MISS 


File Mission Number 

FIT 

= 

File Flighf Number 

LIN 

- 

File Line Number 

RUN 

= 

File Run Number 

LTYPE 

= 

Control Directive 

LSTYPE 

- 

Program Recognized Directives 

lEOR 

- 

lOH End of Reel 

DATA 

= 

Output Data 

EOF2 


EOF Indicator on File Code 02 


III. SUBROUTINE READ2 


READ2 performs fhe record refrieval funcfion ‘for OUTPUT. Oufpuf files • 
are blocked in 20 record groups. READ2 rehrleves fhese blocks and wifhdraws 
records as required by OUTPUT, The end of fhe file indicafor EOF2 Is generafed 
by READ2. A lisfing for rhis rouHne appears in Figure B“3 . 
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iuiikUul IM' ia*L<il ILATA, EOI-(; ) 


COOOl'iOA 


c 



c 

Thii PKCGRAH rfAS PRL'I’AKtU BY 


c 

GLEN t. ELLIOT I 


c 

UNlVtMSnv OF KAN3AS CENTER FCR 


c 

THIS ROLTtNC UNllLCCKS 20 KECORlJ 


c 

r 

AN ENL CF FILE iNUiCArOR, E0F2. 

COOQO? 

u 

OIHEflSICN IUATA<20)» 1BLOCK(20,20) 

toooo 


LOGICAL LCF2 

COOOC5 


DATA NuEX/0/ 

coooos 


E0F'2= .FALSE . 

(OOQO& 


IK INDEX. NE.OJ CC TC 10 

C00006 


N0£X=1 

COQOC 7 


REALtZ) IBLOCK 

oooo 14 


IF(EUF,21 ICO. 10 

CC0020 

10 

IFC 1 JLl-CK( 1 ,NCEX) .Fa.O) GO TO 90 

CG0Q23 


DO 20 1^1.20 

CC002& 

20 

IOATA ( 1 1 =-1 BLOCK ( I ,KDEX 1 

fOOOjS 


NUEX=-NDEX+ 1 

C00036 


IFt.NOEX.GT .20 1 NO£x = 0 

00004 L 


RETURN 

COOOtZ 

VO 

REA012I 

C0004S 


1F1SUF,2) 100,90 

COOOsl 

100 

tUF2=.TRLb. 

COOO!52 


NUtX=0 

C00053 


RET JKN 

00004 


END 


BEStARLHf INC. 
yiCCKS ANU PROVIOei 


COOOlACb 
OOOOlAOt) 
CCOO iAC? 
OOOOiACtJ 
0CCO14CV 
COOOl-tlO 

OCOOJAl 1 
OOOQ14U 
0 GO 014 13 
C00C1A14 
OCOD1415 
0000141G 
C00014 17 
CCCC141* 
0000141^ 
00001420 
00001421 
CCCC1424 
00001423 
00001424 




•FIGURE B-3. SOURCE LISTING FOR SUBROUTINE READ2 
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APPENDIX C 


INSTRUMENT CHARACTERISTICS MAINTENANCE PROGRAM (ICHAR) 


I, INTRODUCTION 

Program ICHAR mainlains a historical file of RADSCAT instrument characteristics. 
New sets of instrument parameters may be submitted to ICHAR for storage; old sets 
may be printed or punched. The design and operation of ICHAR is treated in this 
appendix. The reader will find It helpful to refer to the logic diagram (Figure C-1) 
and listing (Figure C”2) for ICHAR. 


II. DESIGN AND OPERATION 


The files on the historical parameter tape are identified numerically. So upon 
entry into the program the file pointer is initialized to zero. An instruction is read 
and echo checked. The instruction Is then identified as one of three types. If the 
instruction is INITIALIZE, the file appendage counter (NNC) is set equal to one and 
four subsets of instrument parameters are read and transferred to file code 02. The 
format and entry for the parameters were described in Section III D of the text. The 
program is then halted. 

If the historical tape had been initialized already and the cornmand was 
APPEND, the old historical tape would have been rewound (Statement 5) and the contents 
of that tape transferred to another one (Statements between 5 and 18). Once 
transferred, the new set would be appended to the new file (Statements 26 through 23). 
When Insufficient parameters are presented to fill the nev/ file, the program is aborted. 
Otherwise the next control card is read; if it is APPEND, the process is repeated until 
all APPEND commands have been executed. 

If the (next) command had been PRINT or PUNCH, the program would have 
advanced to statement 6 and executed the instruction. The .program -continues until 
all control cards have been exhausted. 
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INITIALIZE 
FILE COL/.-; ][.{ 
(NCS * 01 


INITIALIZE 
APPENDAGE 
COUNTFR 
INCC - NSC) 


REV.'IND 1 

INITIALIZE 
FILE COUNIER 
(NCS - 0) 


/ READ 
/ INSTRUCTIC-J 
HJOBT, NUMBER) i 



i 

ENDFILE 

} 



INCREMENT 
APPE.S'DAGE 
COUNTER 
(NCC-NLC + 1) 


/"read subset i 

^OFINSTR.CHAR./ 


/READ FILE OF! 
IlNSTR. CyAR.j 


/ WRITE 
f instruction) 


identify 

INSTRUCTION 


initialize 


luT )' 


APPEND > -{ 5 

7 /^YES ^ 


WRITER 
OR PUNCH 

\ 7 




/ PRINT SUBSET 
/of INSTR, CHAR, i 


ENPFILE 2 j f^Q 
rr— ^FOUR SETS> 


. append 


[error /Message/ 

(UNRECOGNIZED / 


/ APPEND 
flNSTR. CHAR. 

2 


hVITIALlZr 

^0 


i liir.Tr'**''/'/' I 

j Hkfll I 

FILE COUNTER j 
NCS* NCS +1 I 


PRINT FILE 
6 


[APPEND FILE 

2i 


YES / FILE 


'i.nitiauzeo/ — -(a)op; 

0 / 


^bort up) 


OEIGINAL PAGE JS 
OP POOR QUA 3 M 3 


EOF. 5 
I YES 


'RD\0 INSTR/ 


FIGURE C“I. DESCRIPTIVE FLOW CHART 
FOR PROGRA/v\ ICHAR 


/ ’ — ^ 

(/A BORT IC^ 






































PnOGKAH 1 C HflP ( I NPUT . ODinU F . PlJUCH • I APE 1 • T A(’E2 . 

OOOOI4V5 


o 

tapes* INPUT, 1A('E6=0UT(>UIi tape/. I^PUMDII 

IJ000laV6 


Cl CHAP 

HAINIAIN ichar file. 

f^00014V7 


C 


00001498 


c 

TKIS PROGRAM WAS PHEPARFO HY 

00001499 


c 


00001500 


c 

JOHN P. CLAASSEN , 



c 

OLEN E. ELLUm 

00001501 


c 


00001502 


c 

UNIVERSITY OF KANSAS CENTER FOR RESEARCH, I NC. 

00001S03 


c 


00001505 


c 

THIS PROGRAM WILL STORE RADSCAT INSTRUMENT 



c 

characteristics on a historical tape, the HISTORICAL 



c 

TAPE IS INITIALLY prepared BY USING THE DIRECTIVE 



c 

“INITIALISE” . TO append new sets of characteristics The 


c 

COMMAND “APPEND* WILL TRANSFER OLO SETS FROM TAPF UNIT 



c 

1 TO TAPE UNIT E AND APPEND THE NEW SET , NEW SETS ARE 



c 

READ from CARDS. PREVIOUSLY STORED SETS MAY 



c 

PE PBINFEO OR PUNCHED BY USING THE CONHAf DS PRINT OR 



c 

PUNCH RESPECTIVELY, THE SET NO. HJ3T 8E DESIGNATED . 



c 

THE SETS WILL HAVE BEEN STORED NUMERICALLY IN ORDER ON 



c 

THE TAPE. TO maintain the instrument characteristics 



c 

IT IS advISAOLE to retain The old tape when a 



c 

NEW ONE ONE IS prepared. OLDER versions HAY BE DISCARDED. 


c 

it is advisable also to use a tape comparison utility 



c 

program to verify the transfer. 



c 

NCS 3 file COUNTER FOR THE OLD TAPE ALSO SERVES AS 



c 

A file pointer on unit i 



c 

Ncc » file counter fob the appended sets 



c 

JOBT = EXTERNAL DIPCCTIVE. SEE DATA STATEMENT 



c 

NUMBER = number of THE REQUESTED FILE 



c 

PAR s INSTRUMENT PARAMETERS 



c 

IDAT = DATE 



c 

NUH? s SET NO, AS READ FROM TAPE 


000003 


DIMENSION PARUSO.A). IOATI31. IOPTTA) 

00001S06 

000003 


data 1OPT/10HINITIAL1ZE.6HAPPEN0.SHPRIHT.5HPUNCH/. NCS/0/ 



c 

READ DIRECTIVE AND SET NUMBER. 

OOOOISOB 

000003 

1 

READI5.S00I JOPT. number 

00001509 

000013 


IF (FOF.SI AS. 2 


000016 

500 

FORMATIAIQ.ISI 

OOOOlSll 

000016 

2 

WRITEI5.600) JOPT.NUHHEP 

00001512 

000026 

600 

FORMAT UHl.///“ DIRECTIVE IS “.AlO.lSI 



c 

CHECK DIRECTIVE 

00001S14 

000026 


no 3 N=l,A 

00001515 

000030 


ifijopt.eo.ioptinm go toia.s.s.gi >n 

00001516 

00004] 

3 

CONTINUE 

00001517 

000043 


WRITE(6.601) 

00001516 

000047 

601 

FORMAT!//* UNRECOGNIZED DIRECTIVE*//) 


000047 


CALL ABORT IZHUD) 

OOQ01S2O 


C 

find FILE 


000051 

6 

NCS = NCS.I 


000053 


IF(NUMBER-NCS) 7.8.9 


00005S 

7 

iriNUHBER.GE.ll GO TO 10 

00001525 

000060 


WRITE (5.602) NUMBER 

00001326 

000065 

602 

FORMAT!//" IC SET NUMBER*. 15.* DOES NOT EXIST.*, 



« 

« REQUEST IGNORED.*//) ' 


000065 


GO TO 1 

00001529 

000066 

10 

REWIF'O 1 

00001530 

000070 


«.CS = 0 


000071 


GO TO 6 

00001532 

000072 

12 

NCS = I-l 


000074 


WRITE (G.GOa) NCS, NUMBER 


000104 

608 

FORMAT l//“ THERE ARE 0NLY“I5.“ FILES.* 




• CANT find FIlE*I5//) 


000104 


REWIND 1 

00001534 

000106 


NCS = 0 


000107 


no TO 1 

00001536 


C 

FORWARD SPACE TO CORRECT SET. 

00001537 

000110 

9 

NUM=KUMHER-1 

00001538 

000112 


00 11 I=MCS,NUH 

00001539 

OOOlU 


BEAD(I) 

00001540 

000117 


IFIfOF.l) 12.11 

00001541 

000122 

U 

CONTINUE 

00001542 


c 

READ INSTRUMENT CHARACTERISTICS 

00001543 

00012S 

B 

! = NUMBER 


000127 


REAO(l) NUMZ. lOAT.PAR 


000140 


IFIpdF.U 12.13 

00001545 

000143 

13 

NCS=NUHUER 



C 

WRITE INSTRUMENT CHAHACTFPI SIICS. 

U0001547 

000145 


write 16.610) 


oooiso 

610 

FORMAT (1MI.///I 


OOOISD 


no 15 1 = 1 . A 

(10001548 

000152 


WMITr!6.603) NUH? , IDAT.PAPil . 1 1 .PAH(2. 1 1 .PAR(3, 1 ) 



FIGURE 0-2, 




SOURCE LISTING FOR PROGRAM ICI-IAR 
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0001 7^ 

60J 

rOHMA/<//« 5Et«#Ibi“ UATt**3J*>*tl5.4t5XtAl0.E15,4l 


0001 tit 


no |4* J=4f7b*7 

1)000 1 SSI 

000176 



Il0001S^>2 

000^00 


ir(t.OT.7b> L=75 

0000 1 ss 1 

000^0 i 

14 

Tl' (6*604} (PAR(t(* t 1 tKsJ*U 

UOOOISS'I 


604 

FOttMAT (7E1S,4> 

00001S55 

000232 


rfPITE <6f603) NU 12 i lOAT »PAR (76 f I) * PAR ( 77« IJ «PAR < ]) 


0002AS 


no 15 J»79*160*7, 

0000JS57 

0002A7 


L»J*6 

OOOO JSSrt 

0002S1 


ir(L.GT.l60) L=160 \ 

000 01 ■'SO 

0002ba 

IS 

WRITM6«604MPAH(KiU fKsJiLl ► 

OOOOIS60 

00027b 


jr(N*EO,3) CO TO 1 

0000 ISO I 


C 

PUNCH IC><5. 

00001Sb^ 

000277 


no 17 1=1.4 

00001S63 

000301 


LL*l 

1)0001560 

000302 


write (43.43031 NUM2* IOA7 .PARO *1) .PARCatl) «PAP(3« I) . 




HUHRERt If LL 

00001566 

000333 

4303 

FOR«AT(4I5.10X*E10*3»A10.EIO*3. l2Xf 12.11*15) 


000333 


no 16 J=4.75.7 

00001568 

00033S 


LL=LL*l 

00001S6O 

000337 


L=J.6 

00001570 

OOOJ41 

16 

write (43*4304) (PaR(K. I } fKsJfL) . NUMBER. I *LL 

00001571 

000364 

4304 

rOR*.AT (7E10.3.2X.1 ?*n .15} 


000364 


LL»LL*1 

00001573 

000366 


wRiTf (43*4303) NU42. lOAT .PAR<1 . I > *P AR (2 • I } «P AR ( 3. I } • 

000 


» 

HUMHER* I. LL 

1)0001575 

000416 


no 17 J=79. 160*7 

00001576 

000420 


LL=LL*1 

00001577 

000422 


L» J*6 

00001578 

000424 

17 

WR1tF(43.4304) (PAR(K.l) .KsJ.D .NUHBER.J.LL 

O0001S79 

0004bi 


GO TO 1 

00001580 


C 

, 

00001583 


C 

COPY OLD IC^S 

OOOOlSao 


c 


00001581 

000452 

5 

REwi*TO 1 

0000158S 

000454 


»1CS s: 0 


000455 


WPItE(6.605) 

00001587 

000461 

605 

FORMATUMlf* THE FOLLOWING IC SETS WERE COPIED FROM *. 

00001588 



<»THE OLO file®*/ 

00001589 


a 

® SET NUMBER 10 NUMBER OATS*} 

00001590 

000461 

15 

R£AD(U NUW2*IOATrPAR 

00001591 

000472 


IF(EnP*l) 4.19 

00001592 

000475 

19 

NCS = NC5-*1 



000477 


WR1TEI6.606) NCS.NUM2.IOAT 


000510 

606 

FORMAT (3X.IS*9X*I5*3X.3I5) 

00001599 

000510 


WHITEO) NUH2. JOAT.PAK 

00001595 

OOOS21 


GO TO 18 

00001597 


C 

READ NEW SETS, 

00001598 

000522 

4 

NCC = NCS 


000524 

26 

‘nCC « NCC*1 


000526 


wP1tE(6.607> 


000531 

607 

F0RHAT(1H1/« the following IC SETS WERE APPENDED TO ® 



« 

®THE OUTPUT file.®) 

00001601 

000531 

25 

00 23 1=1 .4 

00001802 

000533 


READ(5»4303> number, ID AT .PAR (1 . I) .PAR (2. I ). PAR ( 3 . 1 1 

00001603 

00055S 


ir(FfiF.5) 40.20 


000560 

20 

WRITE (6.603) NCC.lDATtPARd.I) ,PAR(2.I) *PAR(3*n 


000603 


no 21 J=4.7S*7 * 

00001606 

0006QS 


L=J*6 

U0001607 

000607 


I^(L.GT*75) L*75 

00001603 

000612 


READ (5*4 304) (PARIKt n*K=J.L) 

00001809 

000527 


inroF.S) 40*21 


000632 

21 

WRITF(6«604} <PAR<K.I) .FsJ.L) 

00001611 

0006S1 


READ<S.4303} NU«13£R.I0AT.PAR(76*I) ♦PAR(77.n .PAR(78.1) 

00001612 

000674 


IF(E0F*5} 40.22 


000677 

22 

wRITE( 6.603) NCC*lOAT.PAR(76flJ .PAR (77* I) . PAR(78 . I) 


000722 


no 23 J=79.16 o*7 

00001615 

000724 


L=J*6 

00001616 

000726 


IF(L.OT.160) L=160 

00001617 

000731 


read ( 5.4304) (PAR(K*I) .Ksj.L) 

00001618 

000746 


IF(FOF.S) 40.23 


000751 

23 

write (6.604) (PAR(K.I) «K=J.L> 

0000Lo20 


C 

APPEND IC*S 

00001982 

000772 


WRITE(2) NCC.IDAT.PAR 


001003 


!F(jOKT.NE.IOPT(n ) GO TO 29 


001005 


WRITE(6,800) 


00)011 

»00 

FORMAT (//I 0X®FILE INITIALIZED®//) 


oolon 


STOP 


001013 

29 

REAO<5.SOO) J08T. NUMBER 


001023 


IF (FOF.S) 45.30 


001026 

30 

IF( JOBT .EO. lOPK?) ) 26.35 


001033 

35 

EmijfIlE 2 


001035 


GO TO 2 


00)036 

40 

rfPlU (6,609) 


00)042 

609 

format 1//0 unsatisfied read,®. 



a 

« incomplete slt of CHARaCIFPISTICS*//) 



Figure C~2. (continued) 
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001 0A2 

call abort. (2HIC) 

001 O'. A 

05 ENOriLE ? 

001 OAfj 

STOP 

001 OSD 

END 


00001622 

00001623 


000003 

000011 

oooon 

000013 

000017 

000021 


SUBROUTINE ABORT(A) 

halt program in case of error 

W’lTE (6.100) A 

100 format CIOX. lOriABORT CODE.lX.AlO) 
ENOFILE 2 
R = A/'0.0*1 .0 
STOP 
E'Jt) 




Figure C~2. (Continued) 



III. DEFINITION OF VARIABLES 


The variables employed by ICFIAR are defined in Toble C-1 . 


NSC 

NCC 

JOBT 

NUMBER = 

lOPT 

IDAT 

NUM2 

PAR 


Table C~1 Definilion of Variables 

File Poinfer for File Code 01 
File Poinfer for File Code 02 
Job Insfrucfion 

File Number fo be Prinled or Punched 
Array of Recognized Insfrucflons 
Dale in Paramefers 

Sef Number as Read from File Code 01 
Sef of Insfrumenf Paramefers 
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APPENDIX D 


PROGRAM TDUPE 


1. INTRODUCTION 


Program TDUPE Is a translalion and duplicalion program which prepares an 
inleger version of ihe oulpuf" lape, II is intended ihat copies of fhe inlegerlzed fape 
be senf lo the Invesligators upon requesf for use on olher computing machines. A 
utllify program may be used to prepare additional copies. This appendix describes 
the design and operation of TDUPE. Helpful to the reader are the flow charts of 
Figure D-1 and the program listing of Figure D^2, 


II. DESIGN AND OPERATION 


Action In TDUPE is initiated by an external control card specifying the 
files, by mission and flight, which are to be copied. A file counter is initialized, 
file unit 1 is rewound^ and the instruction, is listed. The labels on the files are 
examined and compared with the instruction. Files not having the designated 
mission and flight are bypassed until one is found. At this time a record counter is 
initialized. Then the instrument parameters are read, scaled, and integerlzed . The 

t 

file heading and Instrument parameters are transferred to file unit 2. The file heading 
is also listed on the printer and the fife counter is Incremented. 

Consecutive records are then read, scaled and integerized. In each instance 
the record counter is Incremented. The records are accumulated in blocks by 
subroutine WRITE2 and transferred to tape unit 2 when the block is filled. The 
terminating block is filled with Integer zeroes when there are insufficient records to 
fill it. When the file has been exhausted, an EOF is appended to the transferred 
file and a message is given to indicate ihe number of records which have been 
transferred. 

The search for additional files continues until the end of'the reel is 
encountered. At this point another instruction is sought. The program terminates 
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FIGURE D-K DESCRIPTIVE LOGIC DIAGRAM 
FOR PROGRAM TDUPE 


OiifGINAri PAGE IS 
OF PCX)R QUALITSy 
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READ DATA 
RECORD 
1 






000003 


000003 

000003 

000003 

000003 

000003 

000003 


000003 

000011 

000011 


OOOOIA 

0000^6 


000030 

000031 

000033 

OOQOAI 

00004! 


0000*»6 

OOOObO 

0000b4 

000061 

000061 


000062 

00006<i 

000070 

000U70 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PM0(J<<AH TDUPf ( ♦IHMI'UT * lAPUl « 1 AI»L?* T AK '»= INPUT* I A|*F TirOUTpUT > 
TAPT. TWAMjLATION AND DUPL I CA T MN POUTiNK 
THIS PWOGHAM WAS PPEPA«CD HY 
JOHrj P. CLAASSEM 

uNivewsiTY or Kansas 

CE'^'TtP row ttESLARCH 


THIS PROGRAH prepares A MASTER OUPI.ICATE TAPE OF THE 
RADSCAT OUTPUT TAPE FOR USE OH OTHER MACHINES, 0THF;P 
COPIES OF THE rtASTfcH DUPLICATE TAPE HAY RE PREPARED 
BY USING AN IN-HOUSL UTILITY PROGRAM* THE OUTPUT DATA 
RECORDS ARE TRANSLATED TO INTEGER FORMAT AFTER 
HULTtPLICAnON BY AN APPROPRIATE SCALE FACTOfi. 

MESSAGES ARE LISTED TO INDICATE THE STATE OF THE 

transfer. 

control Parameters 

HISFUT = MISSION AND FLT OaTA TO BE DUPLICATED 

INPUT/OUTPUT parameters 

IHEAD s FILE HEADING 

lchar = Instrument parameter label 

ichar = instrument characteristics 

DATA * INPUT RECORDS 

IDATA = Integer eouivaeent of data 

lARRAY OUTPUT RECORDS* 

WORKING PARAMETERS 

NOF = NO* OF FILES 

IREC = NO. OF RECORDS IN A FILE 

leoR = •end of reel* 


DIMENSION lHEAD<<9)t ICHAR(l603t DATA(20), 1DATA(20)» 


UOOrn iJfi 
COOOl JJO 
II00 01J40 
UOOOl J4l 

ijOOOl 342 
(lOOOl J43 
[lOOOl J44 
60001 34S 
00001346 
OOOOl 347 
0000 134H 
00001349 
OOOOUbO 
oOOOUbl 
OOOOl JS2 
00001353 
00001354 
00001355 
U0001356 
UOOOlib? 
OOOOUSH 
00001359 
00001360 
00001361 
00001362 
OO0O1363 
00001364 
00001365 
OOOOU66 
00001367 
00001368 
00001369 
00001370 
00001371 
00001372 
00001373 
00001374 
00001375 
00001376 
00001377 
00001373 
00001379 


IARPAY4201 * 1EOH12). MISFLT(2J* LCHAR<4) 

DIMENSION RPAR(75)* SPAR(35N SFACTUS) 

LOGICAL EOFl 

EQUIVALENCE < ICHAR I 1 ) *RPAR ( 1)1 ♦ C ICHAR< 761 i SPAR C 1 )J 
equivalence iOATAUl* lOATA(D) 

DATA SFACT/I.OEIOi 1.0» 1-OESi 1.0E7 t 1.0£*7t 3«l.0t 4*10. 0» 
l.OE3*2*100.0/ 

DATA leOR/lOHEND OF REEilHL/ 

READ the directive 


10 

1000 


?0 


30 

2000 

< 

35 


40 

3000 


SO 

60 


read (SilOOOl 
format 12110) 

IF (EOFiSl 2000 


MISFLT 


WRITE *ENO OF, REEL 

WRITE (2) IEOR* a*I=3*R» 

STOP 

initialize file count 

NOF = 0 

remind 1 

WRIT? (6.2000) H-ISFlT 

FORMAT (IHI.lU* »W1LL PREPARE COPY OF DATA FROM HI SSION* IS* « 

FLT*13///) 

o£A0 (11 IHEAO 

CHECK FOR END OF REEL 

no 40 1 s 1 ,2 

irilHEAOa) lEoR(I)) GO TO 50 

MRlTF (6*3000) NOF 

FORMAT (//IIX.IS* files have BEEN TRANSFERRED *■///) 

GO TO 10 

J 

CHECK for file name 
DO 60 1=1 «? 

IF (IHEAOm *NE. HISFLKl)) GO’ TO 70 ' 

GO tO 90 ' ^ ' 

WPI tF ( 6*4000) (iHEAni n « 1 = ! ,?) « MISFLr 


000013B0 

00001361 

0D0013B2 
00001363 
00001384 
00001365 
00001386 
00001387 
00001388 
00001389 
00001390 
U0001391 
00001392 
00001393 
00001394 
90001395 
00001 196 
00001397 
0000139R 
00001399 
OOOOUDO 
OOOOWOl 

00001^02 

AN000001403 

OOOOU04 

OO0Dl*«O5 

00001<*06 

00001407 

ooooi^oa 

00001469 
00001'*10 
dO0Ol««li 
oDooHia 
OOOOl^ 1 3 
O000l‘*14 
UOOOI-IS 
U0001<* 16 
00001'* 1 7 
0000)418 
U000l4 19 
00001420 


FIGURE D“2. SOURCE LISTING FOR TDUPE 




120 



000104 

4000 

FORMAT (11)C. “FILE HEADING •> 215 " DOES NOT AGREE WITH DIRECTIVE 00001421 


0 

“215) 

00001422 


C 


00001423 


c 

BY-RAS5 UNWANTED FILE 

00001424 


c 


00001425 

000104 

HO 

READ Ul 

00001426 

000107 


IF lEOF.l) 3S.80 

00001427 


c 


00001428 


c 

INITIALIZE RECORD COUNT 

‘ OOOO 1429 


c 


0000 1430 

000112 

90 

IRI-C = 0 

00001431 

00011 J 


READ (11 LCHAR. ICHAR 

00001432 


c 


00001433 


c 

SCALE AND transfer 1N5TR. PARAMETERS 

00001434 


c 


00001435 

000122 


ICHAR(1)=RPAR( l)«l.0E-6 

00001436 

000125 


no 91 1=3.20 

00001437 

000127 

91 

ICHAR (I 1=RPAR( I)“l .0E4 

00001438 

000134 


ICHAR(21)=RPAft(2n«10.0 

00001439 

000137 


1CHAR(22)=RPAR(22)*10.0 

00001440 

000142 


no 92 1=23.75 

00001441 

000143 

92 

ICHAPtl)=RPARin“1.0E4 

OOOO 1442 

000150 


ICHAR (76) =spAR 1 1 ) "1 . OE-6 

00001443 

000153 


no 93 1=3.16 

00001444 

0001S5 

93 

ICHAR 1 1 *75) =SPAR ( I ) <> 1 . 0E4 

00001445 

000162 


no 96 1=17,31 

0000 1446 

000164 

96 

ICHAR (1*75) =SPAR ( I ) “SFACT (1-16) 

00001447 

000171 


no 97 1=32.85 

00001448 

000173 

97 

ICHAR (I*75)=SPAR[I)'»1.0£4 

OOOO 1449 

000200 


write (6,5000) (IHEAD ( I) , 1=1 ,7) 

OOOO 1450 

000212 

5000 

FORMAT (/IIX.OHAVE MISSION»I5»« FLT*I4.* DATE*>3I5* 

FLT LINE®A6, 00001451 


a 

“ RUN«I3) 

00001452 


C 


00001453 


C 

update FILE COUNTER. 

0000 1454 


C 


O00014S5 

000212 


NOF = NOF ♦ 1 

00001456 

000214 


WP1TE(2J IHEAD 

00001457 

000221 


WRITE(2) LCHAR, ICHAR 

00001453 

000230 

95 

call REAOl (DATA.EOFI ) 


000232 


IF (.not-. EOF 1 ) 60 TO no 


000234 

100 

WRITE (6.6000) IREC.NOF 

0000 1461 

000244 

6000 

format (//11X»I4» records have BEEN TRANSFERRED FROM 

FILE® I3///)00001462 

000244 


CALL ENOF2 


006245 


GO TO 35 

O0001464 


c 


■•.0001465 


c 

update RECORD COUNTER. 

OOOOl -66 


c 


00001467 

000246 

no 

IREC = IREC * 1 

00001468 


c 


00001469 


c 

scale AND TRANSFER RADSCAT )(EASUREHENTS. 

00001470 


OOOaSO IARRAY(1)=OATA{1)«100.0 

0002^3 IARRAY(2)=10ATA<2> 

00025A IARRAY(3)=0ATA (3)«100.0 

000257 I APRAY(A)=DArA (4) "100.0 

000261 1A0RAY(5I=IDATA(5) 

000263 IARRAY(6)=I0ATA(5) 

000264 I ARRAY (7) =DATA ( M “10000.0 

000267 IF {DATA{8> .LT. 1.02-03) GO TO 

000272 lARHAY(a) = DATA(8)"1.0 E06 

000275 TaRRAY(9)=1000000 

000276 no TO 130 

000277 120 IA“RAY(a)=r)ATA(d}«1.0C*09 

000302 TAPRAY(9)= 1000000000 

000304 130 IARRAY(10)=OATA<91“100.0 

000307 1ARRAY111)=OATAUO)«100.0 

000312 1ARRAY(12) = D AT A ( 1 1) « 1 0000 . 0 

000315 lARRAY(13>=lOArA(12) 

000316 lAPRAY(14)=OAfA(13)«10.0 

000321 1 ARRAY IIS) =l)At A i 14) <>100.0 

000324 TARHAY(16)=OArAn5)«100.0 

000327 CALL WR I TE2 ( ! ARRA Y ) 

000330 no TO 95 

000331 END 
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00001471 
00001472 
00001473 
00001474 
00001475 
000014 7,6 
O0001477 
00001478 
U0001479 
OOOO 14S0 
o'00014t)l 

<)0001<»rt2 

000014o3 

000014rt4 

00001485 

000014H6 

00001487 

00001483 

00001489 

00001490 

00001491 

00001-*93 

00001494 






Figure D“2. (Continued) 
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000005 


SUBROUTINE READl ( IDATA ,EOFl ) 
DIMENSION IDATA(20), I SLOCK ( 20 ♦ 20 ) 

000005 


LOGICAL FOFl 

000005 


DATA NDEX/0/ 

000005 


EOFl=. FALSE, 

000005 


IF(M(;EX,NE,0) GO TO 10 

000006 


fiOEx = l 

000007 


READ(l) IBLOCK 

000014 


IF(FnF,l) 100,10 

000020 

10 

IF(IRLOCK(l,NDEX) ,EOoO) GO TO 90 

000023 


DO 20 1=1,20 

000025 

20 

IDATA (I) =I BLOCK ( 1 ,NOEX> 

000035 


NDEX=NDEX*1 

000036 


IF (NDEX,GT.20) NDEX=0 

000041 


return 

000042 

90 

READ ( 1 ) 

000045 


IF(E0F,1) 100,90 

000051 

100 

■ EOFl=,TRUE. » 

000052 


NDEX = 0, 1 

000053 


return 

000054 


F'-ID 


000003 

000003 

000003 

00000b 

000015 

000016 

000021 

000026 

000027 

000030 

000036 

000042 

000043 

000044 

000054 

000061 

000062 

000064 

000065 


20 


25 


30 


SUBROUTINE WRITE2 (lARRAY) 
dimension IBLOCK (20?20) ♦ lARRAY(l) 
DATA INDEX /!/ 

DO 20 I = 1»20 

IBLOCK ( I » INDEX ) = lARRAY(I) 

INDEX = INDEX ♦ 1 

IF (INDEX oLEo 20) RETURN 

WRITE (2) I^LOCX 

I'^'OEX = 1 

RETURN 


ENTPY ENDF 2 

IF (INDEX 1) 

DO 25 I = INDEX, 20 

DO 25 J= 1,20 

IBLOCK<J,I.) = 0 ‘ 

WRITE <2> IBLOCK 

INDEX=1 

EN')F1LE 2 

RETURN 

E''>l) 


GO TO 30 




Figure D“2. (Conlinuod) 
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when f'here are no add^^ional ins true Hons. However, before ferminafing an end 
of reel s fare men ^ is wriffen on file unlf 2, 

The variables are defined in Table D“1 . 

Table D“1 Definifion of Variables 


IHEAD 

= 

File Label 

LCHAR 

= 

Insfrumenf Paramefers Label 

ICHAR 


Insfrumenf Paramefers 

DATA - 

- 

Inpuf Records 

IDATA 

= 

Infeger Equivalenf 

lARRAY 

= 

Oufpuf Records 

lEOR 

= 

IIHEND of Reel 

5 FACT 

— 

Scale Faefors 
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APPENDIX E 


SPECIAL PROGRAMS 


I. INTRODUCTION 

Two special engineering roufines WIDTH and GAIN prepare several anfenna 
paramelers essential to the reduction of the scatterometer data. WIDTH computes the 
equivalent beamwidth of the RADSCAT antenna; whereas GAIN computes the antenna 
directivit/. Other parameters of special interest to antenna engineers are also 
computed in GAIN. These programs are provided with the realization that the 
antenna pattern may eventually change through an alteration cf the antenna. In 
this case the equivalent beamwidth and gain should then be recomputed based on 
new pattern information and the instrument parameters changed accordingly. 

Each of the programs are described below. Operating instructions are given in 
Section IV F of the main text. 

II. PROGRAM WIDTH 


The theory for this program was described in Appendix A, Section III B. 
In that development it is apparent that the equivalent beamwidth is given by 


I ^ V n Tf ———————— 

Jo Jo ^ 


d4> 


(1) 


where P(6 <}>) is the normalized power pattern function described in the standard spherica 
coordinate system ( r ^ 0 , ). When the pattern Is elilptically symmetric, the 

integration may be limited to one quadrant in which case 


& 





n Ve o 

. P 


1 


slot 6 d& 


(2) 
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Here It is assumed that the main beam lies on the x“axls as shown in Figure £“1 . 

It Is convenient to work in terms of elevation angle c rather than polar angle 
in which case 

j C0S€ Ae d(p (3) 

Since quadrature techniques usually require rectangular domains of integration, it is 
necessary to apply the transformation 

1|/ = sin £ (4) 

so that 

* * > 

The full limits of integration are not necessary since the integral converges 
rapidly over the main beam. Program WIDTH whose listing is shown in Figure E-2 
use this property to advantage. As shown In the driver rjoutine, the upper limits of 
integration are progressively increased by 1/2 degree intervals, in e. The convergence of 
the equivalent beamwidth can then be observed as the domain increases to include 
more of the beam. 

A Gaussian- Legendre quadrature technique is employed to perform the double 

■ 

integration. This technique is embodied in RINTEG. The reader is referred to Klerer 
and Korn* or any other advanced text on numerically techniques for a discussion of 
this technique , 

The main beam of the RADSCAT antenna pattern is functionally represented 
in the FUNCTION subroutine FUN, RINTEG calls FUN repeatedly to evalute the 
pattern function squared. 


■k 

Klerer, M. and Korn, G. A., "Digital Computer Usei's Handbook," McGraw-Hill 
Book Company/ 1970. 
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V 

A 


FIGURE E-1 . ANTENNA PATTERN GEOMETRY FOR WIDTH 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


program width (OUTPUT . TAPE6=0UTPUT ) 

THIS ROUTINE determines THE EQUIVALENT (PENCIL) 
BEAHWlOTH 

BEAMWlDTH row AN ANTENNA WHOSE PATTERN IS SPECIEIEO 
IN SUBROUTINE PUN. SUUWOUTINE RINTEG INTEGRATES 
THE PATTERN BY A OUASSIAN (LEGENDRE) OUADRATURE METHOD. 
INTEGPAnON ARE REPETlVELY PERFORMED AT ONE DEGREE 
intervals OUT TO SIX DEGREES. THE EOUIVALENI 
BEAMWIOTH is CHOSEN FROM THE LISTING AT THE POINT AT 
' WHICH THE values BEGIN TO CONVERGE 

THIS program was PREPARED BY 

JOHN P. CLAaSSEN 

UNIVERSITY OF KANSAS 
CE''TER FOR RESEARCH. 


OOOOOOO I 
00000002 
OOOOOOOJ 
OOOOOOOA 
00000005 
00000006 
00000007 
UOOOOOOA 
ooooooov 
OOOOOOlO 
UOOOOOll 
00000012 
00000013 
OOOOOOlA 
00000015 
O0DO0016 
00000017 
OOOOOOIB 


C 00000019 


000003 


EXTERNAL FUN 


00000020 

000003 


data PI/3. 14159265/ 

. 

00000021 

000003 


WPITF (6.300> 


00000022 

000007 

300 

format (IHl//) 


00000023 

000007 


WOItE (6.100) 


00000024 

000013 

100 

F0RHAT(30X,5(1H») .’SEARCH FOR THE EQUIVALENT BEAMWIDTHoB ( IH’ 1 //) 

00000025 

000013 


DO 20 1=1.6 


00000026 

0000 15 


X2 = 0.017A532925»FlOaT(I)»0.50 

♦ 

00000027 

000017 


Y2 = X2 


00000026 

000021 


X2 = SIN(X2) 


00000029 

000023 


Xl = 0.0 


00000030 

C0O02A 


Yl = 0.0 


00000031 

00002S 


IX = lO’I 


00000032 

000030 


call RINTEG (FUN. Z, X2. Xli Y2. Yl. IX. IX. 4, 

4> 

00000033 

OOOOAl 


Bfar = SORT (16. O’Z/PI ) /O. 017*532925 


00000034 

OOOOA7 


w»ITE (6*200) I .BEAM 


00000035 

000056 

200 

FORMAT (lOX. ’ACTIVE BEAMWIOTH IS ’I3»’ DEGREES. 

WHEREAS’. 

00000036 



» ’ THE EOUIVALENT BEAMWlDTH IS’FS.2.’ DEGREES’) 


00000037 

000056 

20 

CONTINUE 


00000036 

000060 


STOP 


00000039 

000062 


END 


00000040 


000056 


00 30 I = 1 .HP 

00000083 

000057 


SY(I) = SAMPL£(MP.I)’HDELY 

OOOOOOB* 

000066 

30 

CONTINUE 

00000085 

000070 


DO SO I = 1 .NP 

00000086 

000071 


DO *0 J = 1 »MP 

00000087 

000072 


C(I.J) = C0EF(NP.I)’C0EF{HP,J> 

00000088 

000106 

*0 

CONTINUE 

00000089 

000110 

50 

CONTINUE 

00000090 

000! 12 


DO 90 N=1.NX 

00000091 

00011* 


X1=X2 , 

00000092 

000116 


XZsXl’OELX 

00000093 

000120 


XM = ,X1 ♦ hdelx 

0000009* 

000121 


Y2 = YIP 

00000095 

000122 


DO 60 M=1.MY 

00000096 

00012* 


Yl = Y2 

00000097 

000126 


Y2 = Yl ♦ DELY 

00000096 

000130 


YH = Yl ♦ HDELY 

00000099 

000131 


no 70 1=1. MR 

00000100 

000133 


X=SX(I).XH 

OOOOOIOI 

0001 36 


DO 60 J=1.HP 

O00001O2 

000137 


Y=5Y (J) *YM 

00000103 

0001*2 


SUM = SUH+C ( I . J) »FUN(X. Y) 

0000010* 

OOOIgO 

60 

CONTINUE 

.00000105 

000162 

70 

CONT INUE 

00000106 

00016* 

80 

CONTINUE 

00000107 

000167 

90 

CONTINUE 

00000106 

000171 


7.Z = SUH’D£LX’DELY’0,250 

00000109 

00017* 


RFTURN 

00000110 

OOOl 75 


END 

00000111 


FIGURE E“2. SOURCE LISTING FOR WIDTH 
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000015 

OOOOIS 

000015 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


000015 


5UHR0UT1ME RlUfEG (FUN.2Z. X2P, X H> . Y2P , Y 1 P , IJNX ,HHY,NNP, HHI>I 
THIS ROUTINE EHP|,0Y5 A OUASSIAN LECENDER OUAORAIUPE 

integration Procedure, iniecration over x-V segments 
ARE PERFORHFO AFTER TRANSLATION TO (-1 . 1 ) X (-1 » 1 I . 

F0N= INTEGRAND FUNCTION , 

T7. = EITHER A parameter FOR FUN OR 'THE INTEGRATED RESULT 
OK ROTH 1 

XEP.XIP = UPPER AND LOWER LIMITS ON X 
Y2p,Yip = Upper and uqwer limits '< 

NX. MY = SEGMENTS IN X AND Y. RESPECTIVELY 
NP.HP =5 DEGREE OF PRECISION IN X AND Y, RESPECTIVELY 
DIMENSION SAMPLE (S.Sl .COEF (5 . S3 » SX (5 ) . S Y (5) t C ( 5 . 5) .S V (25) .CV (25) 
EOUIVALENCE (SVni .SAMPLE(l.l) 1 • (CV ( I ) ♦ COEF ( 1 . 1 ) ) 

DATA SV/ 

« 1.0. -0.S773502GO, -0 . 77A596669 . -0.861136312. -0.906179846, 

•> 1.0. 0.S773S0269. 0.0, -0.339981044. -0.533669310. 

« 1.0. 0.0. .776596669. 0.339981066, 0.0, 

« 1.0, 0.0. 0.0. 0.861136312. 0.538469310. 

» 1.0. 0.0, 0.0. 0.0, 0.906179846/ 

DATA CV/ 

<* 1*0. 1.0. 0.55555S556. 0.36785484S. 0.236926885. 

« 1.0, 1.0. 0.088888889, 0.652165155, 0.478628670, 

« 1.0. 0.0. 0.555555556. 0.6S214S155, 0.568888888, 

« 1.0, 0.0, 0.0, 0.347854845, 0.478628670, 

■» 1.0. 0.0, 0.0. 0.0, 0.336926885/ 


000015 

SUM = 0.0 

00000066 

000016 

XI = XlP 

00000067 

000017 

X2 = X2P 

00000063 

000020 

Y1 = YIP 

00000069 

000021 

Y2 = y2P 

00000070 

000022 

NX = NNX 

00000071 

000023 

MY = HMY 

00000072 

000025 

NP = NNP 

00000073 

000026 

HP = HHP 

00000076 

000030 

DELX = (X2-X1 )/FL0AT (NX) 

00000075 

000033 

DELY = (Y2-Y1)/FL0AT(MY} 

00000076 

000037 

HOELX = OELX®0.50 

00000077 

000041 

HDELY = DELX«0.50 

00000078 

000042 

X2 = Xl 

00000079 

000063 

DO 20 I = l.NP 

00000080 

000065 

SX(T) = sample (NP, I )®HDELX 

’ 00000031 

000056 

20 CONTINUE 

00000082 




FUNCTION FUN CXX.YYJ 

00000112 


c' 

THIS ROUTINE DEFINES THE SQUARE OF THE ANJENNA 

U0000U3 


c 

pattern Over the main beam. 

uOOOO 1 14 

000005 


FSD = XX 

OOOOOl 15 

000005 


ESP = ATAM(ESP.SQRT(1.0-ESP"ESP) > 

00000116 

000015 


PKIP =YY 

.000001 17 

000020 


AHXEP = AHAXl (ESP.PHIPj 

00000118 

000024 


IF (AmXEP .lE. 0.0291 1 GO To 10 

OOOOOl 19 

000027 


IF (AMXEP .LE. 0.0349 ) GO TO 20 

00000120 

000031 


IF (AMXEP .LE. 0.0611 ) GO TO 30 

00000121 

000034 


IF (AMXEP .UE. 0.105 ) GO TO 40 

00000122 

000036 


IF (AMXEP .LE. 0.628) GO TO 50 

00000123 

000041 


FUN = J0.0F-3S 

00000124 

000042 

10 

FUN = EXP ( -4690.OESPOESP -3280 .-“PH IP^PHI P ) 

0000012S 

000050 


FUN = FUN»FUN 

0000012b 

OOOOS2 


RPTUP7 

00000127 

000052 

20 

FUN.= EXP (-121-1 "ESP-3280. «PHIPoPH!P> 

00000128 

000060 


FUN = FUN*fUN 

0000012V 

000062 


PFTURM 

-00000130 

000062 

30 

FUN = O,O39a<‘EXP(-021.»ESP»ESP-lO72."PHIP®PHlP) 

00000131 

000071 


FUN = FUN"FUN 

00000132 

000073 


RETURN 

00000133 

000073 

60 

FUN = 0.0159“EXP(-4R6.«ESP»ESP-53.6«PHIP) 

00000136 

000102 


FUN = FUN"FUN 

00000135 

000106 


RFTURN 

U0000l3b 

000106 

SO 

FUN =O.OOO159*EXP(-0.59"SORT IESP«ESP.PH IP»PHI P ) ) 

00000137 

000116 


FUN = FUN"FUN 

00000138 

000116 


RFTURN 

00000139 

000116 


END 

00000140 
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III. PROGMM GAIN 


A . Iniroducfion 

Program GAIN compules ^he anf'enna direcfivlly, beam solid angle and, antenna 
efficiency. The fheory and design of Ihe program is described below. 

B . Theory 

The an^enna beam solid angje is given by 


n TT 
0 


where 


( 6 ) 


P(0,4>) = normalized antenna pattern 
6 = polar angle 
~ azimuthal angle 

p(0,4>) represents the normalized power in the antenna pattern. It can be shown that 
the antenna directivity F is related to by the following expression 


r = 


4- 'Tt' 


( 7 ) 


The antenna efficiency, on the otherhand. Is the ratio of the power in the 
main beam to that in the entire pattern. Clearly then, the efficiency is given by 


pS'iT . 

“ J J P Sen & <10 


(8) 


where 0| is c polar angle which subdivides the main beam. In the above expression 
it Is tacitly assumed that the main beam lies on the positive z axis as Illustrated In 
Figure E"3. 9] is chosen somewhat arbitrarily. Commonly 0] is chosen at the ~3 
dB point or at the first null of the antenna pattern. For narrow beam antenna, the 
latter Is most appropriate for radiometric work. 
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C . Design and Opoi alien oT GAIN 

In ihe discussion which follows if will be helpful lo examine fhc FORTRAN 
lisfing for GAIN os shown in Figure E-4. Upon.cnfry into GAIN, an array BEAM 
in which various polarized beam solid angle faefors will be sfoiod, is cleared. The 
main beam solid angles are fhen compufod for various 0. spaced af 1/2 degree 
infervais. The infegrafion of fhe anfenna patfern is performed by RINTEG which in 
turn calls for polaiizod pattern infoi motion from FUNl . The results are stored In 
BEAM. The accumulated solid beam angle is stored in TOTAL. The remainder of 
the polarized beam solid angle is then compuled in the following statements. 

At this point the anienna gain is compuled and the gain and beam solid angle 
are printed . 

In the remainder of the program the power content in the cross-polarized 
pattern (FUN2) is computed and added to TOTAL, A DO loop terminating on 30 
then computes the antenna efficiency at one degree intervals out to 8 degrees. The 
result is listed. 

Finally, the total beam solid angle, including polarized and cross-polarized 
contributions, is printed. 

D . Subroutines FUNl and FUN2 

The ontenna patterns describee! in FUNl and FUN2 are actually represented 
in terms of another coordinate system,’ The relationship between the two coordinate 
systems is Illustrated in Figure E“3. In the primed system the pattern is described in 
terms of the elevation angle <r and azimuth as shown in Figure E“1 . The 
necessary transformotions between the two are performed in the opening stotements 
of FUNl and FUN2. The latter stalomenls simply evaluate the pattern over one .of 
six domains In FUNl and over five domains in FUN2. 
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PROGRAM GAIN tOU IPOT . T APE:6=0UTPUT I 


00000 l '*2 
00000 !•>) 

c THIS PROGRAM COMPUTES THF ANTENNA GAJN AND EFFECIENCY FROM 000001«A 
c ANTEHMA (POWER) PATTERNS OESCRIRED IN FUN! AND FUN2. I 00000 IsS 
c FUNl OESCHIGES THE P0I.ARI7ED PATTERN WHEREAS FUN3 OESCRIHES THE OOOOOlAh 
C CROSS-POLAR] 7 EO' pattern. > OOOOOIa? 
r ' 00000 lAA 



C 


THIS PMilGRAl* HAS PREPARED GY 


00000149 


C 




OOOOOISO 


c 


JOHN P, CLAASSEN 


oooooisi 


C 




00000152 


c 


UNIVERSITY OF KANSAS CENTER FOR RESEARCH, INC. 


00000153 


c 




00000154 

000003 



DIMENSION BEamUOI. EFFdOl 


OOOOOlbS 

OOOOOJ 



external FUNl 


00000156 

OOOOOJ 



external FUN2 


1*0000157 

000003' 



DATA Pl/3. 14159265/ 


000001S8 

000003 



DATA FACT'lR/0.0 174532925/ 

* 

00000159 

000003 



WRITF (6.5001 


00000160 

000007 


soo 

FORMAT (IHl//) 


OOOOOlol 

000007 



total = 0.0 


00000162 

OOOOlO 



0“T 10 I = 1.10 


0D000163 

000012 



beam (I) = 0.0 


00000164 

000013 


10 

roNTiHue 


00000165 

OOOOIS 



Y2 = IM/2.0 


U0000166 

OOOOlG 



VI = 0.0 


00000167 

000017 



WRITE (o.lOOl 


00000168 

000023 


100 

FORMAT (30X. '•antenna GAIN AND EFFICIENCY FACTORS*//) 


00000169 


c 




00000170 


c 


COMPUTE The main beam solid ANGLE 


00000171 


c 




00000172 

000023 



20 I = I.R 


OOO00173 

00002b 



XI = COS(FI OAT(l)»FACrOR“>O.SO) 


00000174 

000031 



X2 = COS(FLOAT(l-l)»FACTO'?vO.Sn) 


00000175 

0000<»0 



call RINTEG (FifNl. Z. X2. XI. Y2. Yl« 2* 90. 3. 3) 


00000176 

000051 



beam (Itp = MtAM (I) . 4.0®/ 


00000177 

00005S 



TOTAL = TOTAL ♦ 4.0*Z 


00000l7rt 

000057 


20 

CONTINUE 


00000179 


c 




00000180 


c 


COMPUTE The PEHAlNOER REAM SOLID ANGLE 


oooooiai 


c 




00000162 

000061 



X? = COS(4.0’'FACTOR) 


00000183 

00006A 



Xl = COSll?0.0”-FACT''P> 


00000184 

000071 



call PINTEG (FUNl. 7. X2. XI. Y2. Yl. 115. 90. 3. 3) 


00000185 

000102 



TOTAL = TOTAL . 4.0"/ 


00000186 


c 




00000187 


c 


COMPUTE GAIN 


00000188 


c 




00000189 

000105 



GAINl = 4.(1**PI/T0TAL 


00000190 

000107 



GA1h2 = 10,0«ALOG10IGAIN11 


00000191 

000113 



•/RITE ( 6.200I TOTAL. GAlNl. GAIN2 


OOOQOI92 

00012A 


200 

format ( 5X. •POLARIZED ‘beam SOLIO ANGLE=«E 1 4 ,6/ 


00000193 




Sx.opOLARIZED GAIN = *E14.6.» 0R*F6,1,<> DD»//1 


00000194 


c 




00000195 


c 


COMPUTE The C0NTRI8TI0N FROM THE c*<OSS PATTERN 


00000196 


c 




•T0000197 

00012“* 



72 = 1.0 


0000019S 

00012b 



XI = cos (FaCTOH"5.0) 


00000199 

000132 



call HINTEG (FUN2.Z.X2.X1 .Y2.Y1 .S.mO.G.A) 


00000200 

000 IA3 



TOTAL = total . 4.0*2 


U0000201 


c 




00000202 


c 


COMPUTE beam efficiency 


00000203 


c 




OOOOO 2 O 4 

000 1A6 



on 30 1 = 1 .8 


00000205 

000150 



EFF(I> = 8EAM(1«1) /TOTAL 


00000206 

000152 



WRITE (6.300) 1 . £FF( I ) 


00000207 

000162 


300 

format (Sx.*active beamwidth is *13. » degrees. *., 


00000208 



o 

* CORRESPONDING EFFICIENCY 1S*F6.3* ,=) 


OOOQ0209 

OOO 162 


30 

CONTINUE 


00000210 


c 




00000211 


c 


WHITE THE BEAM SOLIO ANGLE 


00000212 


c 




00000213 

OOO 16A 



write (6.400) total 


OOOOO 2 I 4 

000172 


AOO 

format (/5X.*rOTAL HI AM SOLID ANGL£ = “E1 A .6> 


U000021S 

000172 



STOP 


00000216 

00017a 



END 


00000217 


FIGURE E-4, SOURCE LISFING FOR GAIN 
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SUBHOUTIHE HIHTEG (FUN. Z^.X?P«X 1C< y2P. rlP.NNX.KMY.NNP. ««!'> OOOOOZIG 
C IKIS HOUTINE employs A CjUASSIAN LCiENDKE OUADWAfOWf; BOOOftZlV 
C integration PROCtUURE. IHTtGPAIiON OVER X-V SFGMENTS V0000^^0 
c are PENEOHHEU’AFrEH TRANSLATION TO (- I . I | X (- | . I) « 00000221 
C FUN= INTEGRAND FUNCTION 00000222 
C ZZ = EITHER A PARAMETER FOR FUN OR THE INTEGRATED RESULT 00000223 



c 



OR ROTH 

(100002,29 


c 



X2P.XIP => UPPER AND LOWER LIMITS ON X 

00000225 


c 



Y2P.Y1P = UPPER and LOWER LIMITS ON Y 

00000226 


c 



NX.HY = SEGH.NTS IN X AND Y. RESPECTIVELY 

00000227 


c 



NP.HP = DEGREE OF PRECISION IN X AND Y. RESPECTIVELY 

00000228 

000015 



DIMENSION SAMPLE (S. 5) .COEFIS.5) .SXIS) .SYI5I .C(5.5) .SV(2SJ .CV(25l 

00000229 

000015 



FOUTVALENCE (SVU) .SAMPLECl.D). (CV( 1 1 .COEF ( 1 . 1 )| 

00O002J0 

000015 



DATA 

SVZ 

000002J1 



« 

J.0« 

-0.577350269, -0.77A596669, -0.051136312. -0 . 9061 798A6, 

00000232 



o 

uo« 

0.577350269. 0.D« -0.33996I0A9. -0.53B'.69310, 

00000233 



o 

l.Ot 

0.0, .779596659, 0.339981099, 0.0, 

00000239 



« 

1.0« 

0.0. 0.0, 0.861136312, 0.538969310, , 

00000235 



« 

1*0« 

0.0, 0,0, 0.0, 0.906179896/ \ 

00000236 

000015 



DATA 

CV/ , 

00000237 



o 

1*0* 

1.0. 0.555555556. 0.397H5989S. 0.2359268BS. 

00000238 



e 

uo* 

1.0, 0.688888889. 0.652198155.' 0.978628670. 

00000239 



o 

KOt 

0.0, 0.555555556. 0.652195155. 0.568688888, 

00000290 



o 

1.0* 

0.0, 0.0. 0.397859895. 0.978628670. 

00000291 



o 

1*0« 

0.0. 0.0. 0.0. 0.236926885/ 

00000292 


c 




00000293 


c 



CLEAR SUMMING VARIABLE 

00000299 


c 




00000295 

000015 



SUM : 

, .0.0 

00000296 


c 




1)0000297 


c 



RE-ASSIGN INPUT ARGUMENTS 

00000298 


c 




00000299 

000016 



XI = 

XlP - 

00000250 

000017 



X2 » 

X2P 

00000251 

000020 



Y1 :5 

YIP 

00000252 

000021 



Y2 ? 

Y2P 

00000253 

000022 



NX 

NNX 

00000259 

000023 



HY s 

mhy 

00000255 

00002S 



NP :: 

NNP 

00000256 

000026 



HP ti 

HHP 

00000257 


c 




00000258 


c 



COMPUTE LENGTH OF CELL SIDES 

00000259 



C 




00000260 

000030 



oetx = CX2-X1) /FLOAT(NX) 


00000261 

000033 



OELY = <Y2-YI)/FL0AT<HY) 


00000262 

000037 



H0£LX s OELX^O.50 


00000263 

000041 



H0£UY s 0ELX«0.50 


00000269 


C 




00000265 


c 


FORM SAMPLE factor -FOR 

X 

00000266 


c 




00000267 

000042 



no 20 I =i l*NP “ 


00000268 

OOOOA3 



sx<i) = Sample <NP#n*»HOELX 


00000269 

000052 


20 

CONTINUE 


00000270 


c 




00000271 


c 


FORM SAMPLE FACTOR FOR 

Y 

00000272 


c 




00000273 

000054 



DO 30 I = 1*MP ^ 


00000279 

000055 



SY(I) ~ SAMPLE (HP in»HOELY 


00000275 

000064 


30 

continue 


00000276 


c 




00000277 


c 


FORM GAUSSIAN WEIGHTS 


00000278 


c 




00000279 

000066 



DO 50 I - l«NP 


00000260 

000067 



no 40 j s= I iHP 


000002B1 

000070 



C(l*j) « COEFCNP*I)«COEF(HP»J> 


00000262 

000104 


40 

continue . 


00000203 

000106 


50 

continue 


00000269 


c 




00000265 


c 


integrate in strips of 

DELX 

00000266 


c 




00000287 

0001 10 



X2 s XI 


00000208 

000)12 



DO RO Nfil»NX 


00000209 

000113 



XlaXZ 


U0000290 

OOOllS 



X2»X1^0ELX 


00000291 

000117 



XH a XI ♦ HOELX 


00000292 


c 




00000293 


c 


INTEGRATE ALONG Y 


00000299 


c 




00000295 

000120 



Y2 = YIP 


00000296 

000121 



no RO H«ltHY 


0000029; 

000123 



Yl =* V2 


00000298 

00012b 



y 2 Yl ♦ OELY 


00000299 

000127 



YH a Y1 ♦ HDELY 


UOOOOJOO 


c 




OOOOO 101 


c 


TRANSFORM TO CELL <-lt 

i)X(*i »n 

ooooo to? 




V 


Figur<3 E~4, (Continued) 
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c 


000130 



DO 70 1=1. NP 

000132 



X=SX( 1)*XH 

0001 JS 

c 


00 60 3=1. HP 


c 


FORM partial SUM 


c 



000136 



Y=SY(J>*YH 

000141 



SUM = SUH.C(I.J) “FUNlX.YI 

000157 


60 

continue 

000161 


70 

CONTINUE 

000163 


60 

continue 

000166 

c 

QO 

CONTINUE 


c 


APPLY JACOBIAN 

000170 

c 


ZZ = SUH=nELX<’DELY*0.250 

000173 



RETURN 

000174 



END 


UOOO0JO3 
000O0J04 
00000 SOS 
00000 S06 

00000 soo 
00000 soo 

00000310 
00000307 
0000031 1 
0000031^ 
00000313 
00000314 

00000315 

00000316 

00000317 

U000031B 

00000319 

000003^0 

00000321 


000005 

000013 

000016 

000022 

00002'* 

000031 

000034 

0000«5 

000046 

000051 

000052 


000054 

000057 

000062 

000064 

000067 

000071 

000074 

000077 

000106 

000110 

000123 

000124 

000124 

000133 

000133 

000142 

OOOl'*2 

000152 

000152 

000162 

000162 

000173 

000173 

000177 

000201 

000204 

000207 

000224 

000225 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


c 

c 

c 


FUNCTION FUNl IX. Y1 00000322 

00000323 

POLARIZED PATTERN U0000324 

THIS FUNCTION SUBHOUTIME EVALUATES THE PATTERN INTENSITY 00000325 
IN THE PRIMED coordinate SYSTEM. THE PATTERN IS 00000326 

EXPRESSED IM TERMS OF ELEVATION ANO A2IHTH IN THAT 00000327 

COORDINATE SYSTEM, THE INPUT ARGUMENTS. HOWEVER. ARE 00000328 

DESCRIBED IM AN UNPRIMLD COORDINATE SYSTEM RELATED 00000329 

TO THE PRIMED SYSTEM BY A ROTATION ABOUT THE 00000330 

Y AXIS. the BORE SIGHT IS ON THE X-PRIME AXIS IN 00000331 

THE PRIMED SYSTEM ANO ON TH£ Z AXIS IN THE UNPRIHEO 00000332 

SYSTEM. 00000333 

00000334 

00000335 

transform to pattern coordinates 00000336 

00000337 


SINX = SORTU.O - X«X) 

COSY = COb(Y) 

s iny=sqrt ( 1 . o-cosy«cosy) 

COSTP =-SINX«C0SY 

SINTP = SORT U.O-COSTP«COSTPl 

THETAP=ATAN2 (SIMTP.COSTP) 

PHIpr ATANaiSINXASINY.XI 
APHIP = ABS(PHIP) 

AESPP = ABSU.57079633-THETAP) 
ESO = AESPP»AESPP 
oSO = PhIP'*PHIP 


00000333 

00000339 

00000340 

00000341 

00000342 

00000343 

00000344 

00000345 

00000346 

00000347 

00000348 

00000349 


EVALUATE THE PATTERN FUNCTIONS 


00000350 

00000351 


AHXEP = AMAXI (AESPP, PHIPl 
IF (AMXEP .LE. 0.02S1 I GO TO 100 

IF (AMXEP ,l£. 0.0349 » 60 TO 200 * 

IF (AHXEP ,LE. O.06U ) 00 TO 300 

IF (AHXEP ,LE. '0,105 ).G0 TO 400 
IF (AHXEP .l£. 0.628 ) GO TO 500 i 
XI = -SIHTP*COS(PHlP) I 

XI = ATAN2(SORT(1.0-XI»XI) .Xll 
XlD = XI*S7. 2957796 

IF (XID ,LT. 100. .AND. XIO ,GT. 60. > GO TO 600 

FUNl =0.0 

RETURN 

100 FUNl = EXP < -4690. «eSO -3280. "PSO ) 

return 

200 FUNl =EXP(- 121.1«AESPP-3280,»PSQ) 

return 

300 FUNl = 0.039e<’£XP(-820.7**ES0-1072.'‘PSO) 

RETURN 

400 FUNl = 0.0159»6XP(-486.»ESO-53.6»APHIP» 

RFTuRN 

SOO FUNl = O.OOOIBO-EXPI-S.SROSORT (FSO.PSO) 1 

RETURN 

600 TANPSl = SINTP/COSTP»SIN(PHIP) 

PSl = ATANdANPSn 
COSPS = COS (PSI ) "’*2 
SINPS = SIN(PSI)*“2 

FUNl = (3.16E-0SOCOSPS ♦ 6 , 31E-06»SINPSI *• 

« EXP(-(Xin-8o.l»'*2/C10G. '•COSPS . 36.»S1NPS1) 
RFTURN 
E"0 


00000352 

00000353 

00000354 

00000355 

((0000356 

00000357 

00000356 

00000359 

00000360 

00000361 

00000362 

00000363 

00000364 

00000365 

00000366 

00000367 

00000368 

00000369 

00000370 

O0000371 

1)0000372 

OODOOJ73 

00000374 

00000375 

00000376 

00000377 

000003/6 

UOOOOJ'9 

OOOOOJAO 
IIOOOO 381 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION FUN? tX»Y) 0000038? 

00000333 

CROSS-POL AH I ZFn PATTERN 0000038A 

' THIS FUNCTION SUBROUTINE EVALUATES THE PATTERN INTENSITY 0000038S 
IN THE PRIMED COORDINATE SYSTEM. THE PATTERN IS 00000336 
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Y AXIS. IHE bore sight is on the X-PRIHE axis in 00000391 

THE PHIHED SYSTEM AND ON THE Z AXIS IN THE UNPRIMED 0000039? 

SYSTEM. 00000393 

0000039A 

00000395 

TRANSFORM TO PATTERN COORDINATES 00000396 

00000397 


000003 



SINX = SORTd.O - X»X) 

00000398 

000013 



COSY = COSlYl 

00000399 

000016 



SINY-SQHT( 1 .0-COSY*COSY) 

00000900 

0000?2 



rOSTP =“SIMX«COSY 

00000901 

OOOOZA 



SINTP = SORTn.O-COSTP«COSTP) 

00000902 

000031 



THETAP=ATAN2(51NTP.C0STP) 

00000903 

000039 



PHIP= ATAN2(SINX»SINY.X> 

00000909 

OOOO^S 



APHIP = ABSIPHIP) 

00000905 

000096 



AESPP = A8SI1.S7079633-TH£TAP> 

00000906 

OOOOSl 



ESO = AESPPOAESPP 

00000907 

000052 



PSQ = PHIP»PHIP 

OOOOO 9 O 8 


C 



00000909 


C 


EVALUATE THE PATTERN FUNCTIONS 

00000910 


C 


■ - 

00000911 

000059' 



AMXEP = AHAXl (AESPP, PHIP) 

00000912 

000057 



IF (AMXEP ,LE. 0.0291 1 GO TO 100 

00000913 

000062 



TF (AMXEP .LE. 0.0399 ) GO TO 200 

00000919 

000069 



IF (AMXEP .LE. 0.0611 ) GO TO 300 

00000915 

000067 



IF (AMXEP ,LE. 0.105 ) GO TO 900 

O00O0916 

000071 



FIIN2 = 0.0 

00000^17 

000072 



pftupn 

00000915 

000073 


100 

FUN2 = 0.0063'EXP < -9690 ."ESQ-3280 .*PSO) 

0000091 9 

000103 



return 

00000920 

000103 


200 

FUN2 = 0.0063»EXP ( -9690. *ESO -3280. «PSO ) 

00000921 

OOOl 13 



rfturn 

OOOOOh22 

0001 13 


300 

FUN2 = 0.00025»EXP(-821.»PSO-B21.*ESQ) 

00000923 


000123 

RETURN 

00000929 

000123 

900 FUN2 = 0.00025"EXP (-621 ."PSQ-821 ,»ESO) 

00000925 

000133 

RETURN 

00000926 

000133 

END 

00000927 
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